Human vs Bushbaby LastZ results

Human (Homo sapiens, GRCh38) and Bushbaby (Otolemur garnettii, OtoGar3) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 76. 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

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

Statistics over 2,077,909 alignment blocks

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

Uncovered: 1,532,950,295 out of 3,096,649,726
Covered: 1,563,699,431 out of 3,096,649,726

Uncovered: 868,633 out of 35,228,305
Matches: 30,241,264 out of 35,228,305
Mismatches: 3,866,916 out of 35,228,305
Insertions: 251,492 out of 35,228,305
Identity over aligned base-pairs: 88.0%

Bushbaby

Uncovered: 1,060,491,436 out of 2,519,404,493
Covered: 1,458,913,057 out of 2,519,404,493

Uncovered: 1,607,429 out of 31,659,618
Matches: 26,849,703 out of 31,659,618
Mismatches: 3,104,789 out of 31,659,618
Insertions: 97,697 out of 31,659,618
Identity over aligned base-pairs: 89.3%

Block size distribution

Size range All 2,077,909 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
50,653
332.7 kb
10 bp - 100 bp
219,315
12.1 Mb
21,752
1.5 Mb
100 bp - 1 kb
1,271,682
561.9 Mb
63,172
22.3 Mb
1 kb - 10 kb
535,407
1.0 Gb
18,218
44.6 Mb
10 kb - 100 kb
852
11.1 Mb
1,569
38.4 Mb
100 kb - 1 Mb
223
87.0 Mb
1 Mb - 10 Mb
241
798.8 Mb
10 Mb - 100 Mb
38
637.6 Mb