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 size10,100,000
Masking options{default_soft_masking => 1}{default_soft_masking => 1}

Statistics over 575,265 alignment blocks

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

Uncovered: 536,242,045 out of 3,099,750,718
Covered: 2,563,508,673 out of 3,099,750,718

Uncovered: 1,663,660 out of 38,571,754
Matches: 35,972,812 out of 38,571,754
Mismatches: 847,600 out of 38,571,754
Insertions: 87,682 out of 38,571,754
Identity over aligned base-pairs: 97.5%

Gibbon

Uncovered: 493,043,003 out of 2,962,060,179
Covered: 2,469,017,176 out of 2,962,060,179

Uncovered: 637,349 out of 33,058,749
Matches: 31,732,120 out of 33,058,749
Mismatches: 654,751 out of 33,058,749
Insertions: 34,529 out of 33,058,749
Identity over aligned base-pairs: 97.9%

Block size distribution

Size range All 575,265 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
8,727
48.1 kb
10 bp - 100 bp
41,884
2.2 Mb
12,794
742.1 kb
100 bp - 1 kb
131,016
58.7 Mb
28,984
9.3 Mb
1 kb - 10 kb
314,910
1.2 Gb
9,601
31.3 Mb
10 kb - 100 kb
78,725
1.3 Gb
3,080
89.1 Mb
100 kb - 1 Mb
3
350.5 kb
351
82.9 Mb
1 Mb - 10 Mb
48
207.3 Mb
10 Mb - 100 Mb
73
2.1 Gb
100 Mb - 1 Gb
1
100.2 Mb