Zebrafish vs Fugu LastZ results

Zebrafish (Danio rerio, GRCz11) and Fugu (Takifugu rubripes, fTakRub1.2) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 99. Zebrafish 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

ParameterZebrafishFugu
Chunk size10,000,00010,100,000
Overlap0100,000
Group set size010,100,000
Masking options

Statistics over 292,593 alignment blocks

Genome coverage (bp) Coding exon coverage (bp)
Zebrafish

Uncovered: 1,318,373,344 out of 1,373,471,384
Covered: 55,098,040 out of 1,373,471,384

Uncovered: 10,172,720 out of 42,031,845
Matches: 22,555,448 out of 42,031,845
Mismatches: 8,816,168 out of 42,031,845
Insertions: 487,509 out of 42,031,845
Identity over aligned base-pairs: 70.8%

Fugu

Uncovered: 327,598,218 out of 384,110,215
Covered: 56,511,997 out of 384,110,215

Uncovered: 8,774,941 out of 37,455,454
Matches: 20,276,388 out of 37,455,454
Mismatches: 7,784,100 out of 37,455,454
Insertions: 620,025 out of 37,455,454
Identity over aligned base-pairs: 70.7%

Block size distribution

Size range All 292,593 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
1,940
8.9 kb
10 bp - 100 bp
47,130
3.4 Mb
13,633
1.0 Mb
100 bp - 1 kb
237,850
59.0 Mb
64,166
20.4 Mb
1 kb - 10 kb
5,658
9.2 Mb
14,084
35.5 Mb
10 kb - 100 kb
14
261.0 kb
688
13.7 Mb
100 kb - 1 Mb
1
116.8 kb
11
1.5 Mb