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

Statistics over 745,885 alignment blocks

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

Uncovered: 554,831,132 out of 3,096,649,726
Covered: 2,541,818,594 out of 3,096,649,726

Uncovered: 397,210 out of 35,628,222
Matches: 34,045,742 out of 35,628,222
Mismatches: 1,119,136 out of 35,628,222
Insertions: 66,134 out of 35,628,222
Identity over aligned base-pairs: 96.6%

Macaque

Uncovered: 520,982,006 out of 2,971,314,966
Covered: 2,450,332,960 out of 2,971,314,966

Uncovered: 2,773,685 out of 37,297,175
Matches: 33,223,024 out of 37,297,175
Mismatches: 1,227,332 out of 37,297,175
Insertions: 73,134 out of 37,297,175
Identity over aligned base-pairs: 96.2%

Block size distribution

Size range All 745,885 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
16,098
90.3 kb
10 bp - 100 bp
54,131
2.7 Mb
13,172
776.5 kb
100 bp - 1 kb
189,049
87.0 Mb
29,908
9.3 Mb
1 kb - 10 kb
423,009
1.6 Gb
8,044
25.8 Mb
10 kb - 100 kb
63,552
990.8 Mb
2,501
79.3 Mb
100 kb - 1 Mb
46
5.8 Mb
467
105.4 Mb
1 Mb - 10 Mb
33
82.4 Mb
10 Mb - 100 Mb
33
1.4 Gb
100 Mb - 1 Gb
8
988.8 Mb