Japanese medaka HdrR vs Golden-line barbel LastZ results

Japanese medaka HdrR (Oryzias latipes, ASM223467v1) and Golden-line barbel (Sinocyclocheilus grahami, SAMN03320097.WGS_v1.1) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 99. 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 HdrRGolden-line barbel
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 239,819 alignment blocks

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

Uncovered: 673,364,872 out of 734,057,086
Covered: 60,692,214 out of 734,057,086

Uncovered: 10,429,783 out of 40,334,262
Matches: 21,139,188 out of 40,334,262
Mismatches: 8,152,208 out of 40,334,262
Insertions: 613,083 out of 40,334,262
Identity over aligned base-pairs: 70.7%

Golden-line barbel

Uncovered: 1,696,344,923 out of 1,750,287,761
Covered: 53,942,838 out of 1,750,287,761

Uncovered: 43,556,178 out of 69,134,296
Matches: 18,279,943 out of 69,134,296
Mismatches: 6,859,198 out of 69,134,296
Insertions: 438,977 out of 69,134,296
Identity over aligned base-pairs: 71.5%

Block size distribution

Size range All 239,819 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
2,680
12.0 kb
10 bp - 100 bp
39,144
2.7 Mb
8,653
631.8 kb
100 bp - 1 kb
192,013
49.8 Mb
44,005
14.2 Mb
1 kb - 10 kb
5,934
10.6 Mb
11,910
31.8 Mb
10 kb - 100 kb
48
676.6 kb
779
15.8 Mb
100 kb - 1 Mb
10
1.3 Mb