Human vs Prairie vole LastZ results

Human (Homo sapiens, GRCh38) and Prairie vole (Microtus ochrogaster, MicOch1.0) 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)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

ParameterHumanPrairie vole
Chunk size10,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,386 alignment blocks

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

Uncovered: 2,096,567,724 out of 3,096,649,726
Covered: 1,000,082,002 out of 3,096,649,726

Uncovered: 1,461,299 out of 35,387,657
Matches: 28,108,599 out of 35,387,657
Mismatches: 5,411,483 out of 35,387,657
Insertions: 406,276 out of 35,387,657
Identity over aligned base-pairs: 82.9%

Prairie vole

Uncovered: 1,379,865,667 out of 2,287,357,235
Covered: 907,491,568 out of 2,287,357,235

Uncovered: 0 out of 11,336
Matches: 8,042 out of 11,336
Mismatches: 3,235 out of 11,336
Insertions: 59 out of 11,336
Identity over aligned base-pairs: 70.9%

Block size distribution

Size range All 1,495,386 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
22,172
150.5 kb
10 bp - 100 bp
144,530
7.7 Mb
10,852
741.7 kb
100 bp - 1 kb
990,103
457.9 Mb
41,621
15.7 Mb
1 kb - 10 kb
338,362
582.8 Mb
14,706
35.0 Mb
10 kb - 100 kb
219
2.8 Mb
923
20.9 Mb
100 kb - 1 Mb
267
109.8 Mb
1 Mb - 10 Mb
211
685.4 Mb
10 Mb - 100 Mb
13
183.9 Mb