Human vs Shrew mouse LastZ results

Human (Homo sapiens, GRCh38) and Shrew mouse (Mus pahari, PAHARI_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

ParameterHumanShrew 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,495,457 alignment blocks

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

Uncovered: 2,104,471,814 out of 3,096,649,726
Covered: 992,177,912 out of 3,096,649,726

Uncovered: 1,426,356 out of 35,437,039
Matches: 28,111,384 out of 35,437,039
Mismatches: 5,493,993 out of 35,437,039
Insertions: 405,306 out of 35,437,039
Identity over aligned base-pairs: 82.7%

Shrew mouse

Uncovered: 1,571,439,642 out of 2,475,012,951
Covered: 903,573,309 out of 2,475,012,951

Uncovered: 1,446,055 out of 32,578,444
Matches: 26,177,650 out of 32,578,444
Mismatches: 4,793,256 out of 32,578,444
Insertions: 161,483 out of 32,578,444
Identity over aligned base-pairs: 84.1%

Block size distribution

Size range All 1,495,457 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
20,514
139.0 kb
10 bp - 100 bp
143,468
7.7 Mb
11,093
752.6 kb
100 bp - 1 kb
995,301
461.2 Mb
41,774
15.3 Mb
1 kb - 10 kb
336,011
574.9 Mb
13,786
32.4 Mb
10 kb - 100 kb
163
2.1 Mb
887
19.8 Mb
100 kb - 1 Mb
160
60.6 Mb
1 Mb - 10 Mb
128
452.7 Mb
10 Mb - 100 Mb
27
464.6 Mb