Human vs Ryukyu mouse LastZ results

Human (Homo sapiens, GRCh38) and Ryukyu mouse (Mus caroli, CAROLI_EIJ_v1.1) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 90. 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)
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

ParameterHumanRyukyu mouse
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 1,498,157 alignment blocks

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

Uncovered: 2,097,857,953 out of 3,096,649,726
Covered: 998,791,773 out of 3,096,649,726

Uncovered: 1,420,929 out of 35,437,039
Matches: 28,110,109 out of 35,437,039
Mismatches: 5,499,361 out of 35,437,039
Insertions: 406,640 out of 35,437,039
Identity over aligned base-pairs: 82.6%

Ryukyu mouse

Uncovered: 1,642,896,636 out of 2,553,112,587
Covered: 910,215,951 out of 2,553,112,587

Uncovered: 1,595,746 out of 32,770,133
Matches: 26,198,236 out of 32,770,133
Mismatches: 4,815,420 out of 32,770,133
Insertions: 160,731 out of 32,770,133
Identity over aligned base-pairs: 84.0%

Block size distribution

Size range All 1,498,157 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
21,117
142.6 kb
10 bp - 100 bp
143,896
7.7 Mb
11,591
789.4 kb
100 bp - 1 kb
994,543
461.4 Mb
42,286
15.6 Mb
1 kb - 10 kb
338,400
581.0 Mb
13,780
32.5 Mb
10 kb - 100 kb
201
2.5 Mb
859
19.7 Mb
100 kb - 1 Mb
156
59.2 Mb
1 Mb - 10 Mb
129
453.3 Mb
10 Mb - 100 Mb
27
471.6 Mb