Human vs Steppe mouse LastZ results

Human (Homo sapiens, GRCh38) and Steppe mouse (Mus spicilegus, MUSP714) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 96. 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)3000
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

ParameterHumanSteppe 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,469,949 alignment blocks

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

Uncovered: 2,095,017,308 out of 3,096,649,726
Covered: 1,001,632,418 out of 3,096,649,726

Uncovered: 1,503,406 out of 35,604,148
Matches: 28,056,612 out of 35,604,148
Mismatches: 5,633,467 out of 35,604,148
Insertions: 410,663 out of 35,604,148
Identity over aligned base-pairs: 82.3%

Steppe mouse

Uncovered: 1,586,052,209 out of 2,496,083,520
Covered: 910,031,311 out of 2,496,083,520

Uncovered: 3,656,119 out of 35,178,654
Matches: 26,342,686 out of 35,178,654
Mismatches: 4,988,448 out of 35,178,654
Insertions: 191,401 out of 35,178,654
Identity over aligned base-pairs: 83.6%

Block size distribution

Size range All 1,469,949 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
21,179
143.6 kb
10 bp - 100 bp
133,710
7.1 Mb
8,785
612.0 kb
100 bp - 1 kb
974,467
455.8 Mb
43,103
16.8 Mb
1 kb - 10 kb
340,229
587.1 Mb
22,118
57.4 Mb
10 kb - 100 kb
363
5.1 Mb
2,026
60.2 Mb
100 kb - 1 Mb
1
100.3 kb
1,151
444.5 Mb
1 Mb - 10 Mb
283
475.7 Mb