Human vs Macaque LastZ results

Human (Homo sapiens, GRCh38) and Macaque (Macaca mulatta, Mmul_10) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 98. 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

ParameterHumanMacaque
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 714,896 alignment blocks

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

Uncovered: 557,932,124 out of 3,099,750,718
Covered: 2,541,818,594 out of 3,099,750,718

Uncovered: 567,749 out of 38,571,754
Matches: 36,625,147 out of 38,571,754
Mismatches: 1,279,430 out of 38,571,754
Insertions: 99,428 out of 38,571,754
Identity over aligned base-pairs: 96.4%

Macaque

Uncovered: 522,663,791 out of 2,971,331,530
Covered: 2,448,667,739 out of 2,971,331,530

Uncovered: 2,789,803 out of 37,308,545
Matches: 33,343,167 out of 37,308,545
Mismatches: 1,101,593 out of 37,308,545
Insertions: 73,982 out of 37,308,545
Identity over aligned base-pairs: 96.6%

Block size distribution

Size range All 714,896 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
15,380
86.4 kb
10 bp - 100 bp
50,539
2.5 Mb
12,021
707.5 kb
100 bp - 1 kb
177,347
81.8 Mb
27,381
8.5 Mb
1 kb - 10 kb
408,931
1.6 Gb
7,260
23.2 Mb
10 kb - 100 kb
62,670
976.4 Mb
2,168
65.9 Mb
100 kb - 1 Mb
29
3.6 Mb
367
85.8 Mb
1 Mb - 10 Mb
21
51.6 Mb
10 Mb - 100 Mb
33
1.4 Gb
100 Mb - 1 Gb
8
988.8 Mb