Japanese medaka HdrR vs Burton's mouthbrooder LastZ results

Japanese medaka HdrR (Oryzias latipes, ASM223467v1) and Burton's mouthbrooder (Haplochromis burtoni, AstBur1.0) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 94. 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 HdrRBurton's mouthbrooder
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 353,508 alignment blocks

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

Uncovered: 565,804,233 out of 734,057,086
Covered: 168,252,853 out of 734,057,086

Uncovered: 5,743,889 out of 40,334,262
Matches: 27,216,542 out of 40,334,262
Mismatches: 7,101,223 out of 40,334,262
Insertions: 272,608 out of 40,334,262
Identity over aligned base-pairs: 78.7%

Burton's mouthbrooder

Uncovered: 658,604,644 out of 831,411,547
Covered: 172,806,903 out of 831,411,547

Uncovered: 4,914,739 out of 37,764,289
Matches: 26,008,309 out of 37,764,289
Mismatches: 6,326,800 out of 37,764,289
Insertions: 514,441 out of 37,764,289
Identity over aligned base-pairs: 79.2%

Block size distribution

Size range All 353,508 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
1,210
7.4 kb
10 bp - 100 bp
26,719
1.7 Mb
2,585
188.3 kb
100 bp - 1 kb
281,818
110.8 Mb
16,388
6.0 Mb
1 kb - 10 kb
43,627
71.4 Mb
4,809
12.4 Mb
10 kb - 100 kb
134
2.4 Mb
903
29.5 Mb
100 kb - 1 Mb
347
105.1 Mb
1 Mb - 10 Mb
22
33.1 Mb