Human vs Horse LastZ results

Human (Homo sapiens, GRCh38) and Horse (Equus caballus, EquCab3.0) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 95. Human was used as the reference species. After running LastZ, the raw LastZ alignment blocks were chained according to their location in both genomes. During the final netting process, the best sub-chain was chosen in each region on the reference species.

Configuration parameters

ParameterValue
Gap open penalty (O)400
Gap extend penalty (E)30
HSP threshold (K)3000
Threshold for gapped extension (L)3000
Threshold for alignments between gapped alignment blocks (H)2200
Masking count (M)
Seed and Transition value (T)1
Scoring matrix (Q)
Default:
     A    C    G    T
    91 -114  -31 -123
  -114  100 -125  -31
   -31 -125  100 -114
  -123  -31 -114   91

Chunking parameters

ParameterHumanHorse
Chunk size30,000,00010,100,000
Overlap0100,000
Group set size010,100,000
Masking options{default_soft_masking => 1}{default_soft_masking => 1}

Statistics over 2,108,479 alignment blocks

Genome coverage (bp) Coding exon coverage (bp)
Human

Uncovered: 1,376,321,059 out of 3,096,649,726
Covered: 1,720,328,667 out of 3,096,649,726

Uncovered: 656,977 out of 35,570,071
Matches: 30,785,219 out of 35,570,071
Mismatches: 3,899,140 out of 35,570,071
Insertions: 228,735 out of 35,570,071
Identity over aligned base-pairs: 88.2%

Horse

Uncovered: 881,239,329 out of 2,506,966,135
Covered: 1,625,726,806 out of 2,506,966,135

Uncovered: 3,547,903 out of 37,218,147
Matches: 29,642,954 out of 37,218,147
Mismatches: 3,790,263 out of 37,218,147
Insertions: 237,027 out of 37,218,147
Identity over aligned base-pairs: 88.0%

Block size distribution

Size range All 2,108,479 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
86,533
583.1 kb
10 bp - 100 bp
276,934
12.6 Mb
20,734
1.4 Mb
100 bp - 1 kb
1,118,241
510.5 Mb
61,913
21.8 Mb
1 kb - 10 kb
624,860
1.3 Gb
20,187
53.1 Mb
10 kb - 100 kb
1,895
27.0 Mb
2,385
65.1 Mb
100 kb - 1 Mb
16
2.0 Mb
283
72.5 Mb
1 Mb - 10 Mb
63
209.1 Mb
10 Mb - 100 Mb
48
1.5 Gb