Human vs Gibbon LastZ results

Human (Homo sapiens, GRCh38) and Gibbon (Nomascus leucogenys, Nleu_3.0) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 91. 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)5000
Threshold for gapped extension (L)5000
Threshold for alignments between gapped alignment blocks (H)3000
Masking count (M)10
Seed and Transition value (T)1
Scoring matrix (Q)
Primate:
    A    C    G    T
   90 -330 -236 -356
 -330  100 -318 -236
 -236 -318  100 -330
 -356 -236 -330   90

Chunking parameters

ParameterHumanGibbon
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 551,209 alignment blocks

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

Uncovered: 533,141,053 out of 3,096,649,726
Covered: 2,563,508,673 out of 3,096,649,726

Uncovered: 1,453,644 out of 35,444,608
Matches: 33,187,034 out of 35,444,608
Mismatches: 737,339 out of 35,444,608
Insertions: 66,591 out of 35,444,608
Identity over aligned base-pairs: 97.6%

Gibbon

Uncovered: 496,719,118 out of 2,962,060,179
Covered: 2,465,341,061 out of 2,962,060,179

Uncovered: 646,737 out of 33,058,749
Matches: 31,727,849 out of 33,058,749
Mismatches: 649,611 out of 33,058,749
Insertions: 34,552 out of 33,058,749
Identity over aligned base-pairs: 97.9%

Block size distribution

Size range All 551,209 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
8,391
46.3 kb
10 bp - 100 bp
39,162
2.0 Mb
12,136
702.6 kb
100 bp - 1 kb
123,338
55.3 Mb
27,233
8.7 Mb
1 kb - 10 kb
303,218
1.2 Gb
8,474
27.5 Mb
10 kb - 100 kb
77,100
1.3 Gb
2,630
75.2 Mb
100 kb - 1 Mb
1
107.1 kb
158
36.2 Mb
1 Mb - 10 Mb
45
201.3 Mb
10 Mb - 100 Mb
73
2.1 Gb
100 Mb - 1 Gb
1
100.2 Mb