Human vs Golden Hamster LastZ results

Human (Homo sapiens, GRCh38) and Golden Hamster (Mesocricetus auratus, MesAur1.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

ParameterHumanGolden Hamster
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,458,823 alignment blocks

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

Uncovered: 2,132,401,956 out of 3,096,649,726
Covered: 964,247,770 out of 3,096,649,726

Uncovered: 2,458,873 out of 35,387,657
Matches: 27,261,250 out of 35,387,657
Mismatches: 5,264,132 out of 35,387,657
Insertions: 403,402 out of 35,387,657
Identity over aligned base-pairs: 82.8%

Golden Hamster

Uncovered: 1,634,584,393 out of 2,504,925,039
Covered: 870,340,646 out of 2,504,925,039

Uncovered: 0 out of 11,328
Matches: 7,884 out of 11,328
Mismatches: 3,376 out of 11,328
Insertions: 68 out of 11,328
Identity over aligned base-pairs: 69.6%

Block size distribution

Size range All 1,458,823 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
21,951
148.6 kb
10 bp - 100 bp
143,525
7.7 Mb
11,311
785.4 kb
100 bp - 1 kb
970,227
446.7 Mb
49,632
18.9 Mb
1 kb - 10 kb
322,917
555.5 Mb
16,958
40.4 Mb
10 kb - 100 kb
203
2.5 Mb
1,061
23.6 Mb
100 kb - 1 Mb
279
115.6 Mb
1 Mb - 10 Mb
240
716.2 Mb
10 Mb - 100 Mb
8
97.0 Mb