Japanese medaka HdrR vs Huchen LastZ results

Japanese medaka HdrR (Oryzias latipes, ASM223467v1) and Huchen (Hucho hucho, ASM331708v1) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 97. Japanese medaka HdrR 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

ParameterJapanese medaka HdrRHuchen
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 314,800 alignment blocks

Genome coverage (bp) Coding exon coverage (bp)
Japanese medaka HdrR

Uncovered: 650,951,881 out of 734,057,086
Covered: 83,105,205 out of 734,057,086

Uncovered: 8,524,789 out of 40,334,262
Matches: 23,339,626 out of 40,334,262
Mismatches: 7,989,829 out of 40,334,262
Insertions: 480,018 out of 40,334,262
Identity over aligned base-pairs: 73.4%

Huchen

Uncovered: 2,406,956,094 out of 2,487,549,814
Covered: 80,593,720 out of 2,487,549,814

Uncovered: 42,184,910 out of 70,924,295
Matches: 21,049,696 out of 70,924,295
Mismatches: 6,982,524 out of 70,924,295
Insertions: 707,165 out of 70,924,295
Identity over aligned base-pairs: 73.2%

Block size distribution

Size range All 314,800 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
4,053
18.9 kb
10 bp - 100 bp
48,564
3.2 Mb
8,300
604.7 kb
100 bp - 1 kb
252,556
69.1 Mb
46,230
15.1 Mb
1 kb - 10 kb
9,580
16.9 Mb
11,611
32.6 Mb
10 kb - 100 kb
47
521.0 kb
1,284
31.3 Mb
100 kb - 1 Mb
64
10.1 Mb