Zebrafish vs Guppy LastZ results

Zebrafish (Danio rerio, GRCz11) and Guppy (Poecilia reticulata, Guppy_female_1.0_MT) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 94. 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

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

Statistics over 283,415 alignment blocks

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

Uncovered: 1,313,482,546 out of 1,373,471,384
Covered: 59,988,838 out of 1,373,471,384

Uncovered: 10,308,735 out of 42,031,845
Matches: 22,300,919 out of 42,031,845
Mismatches: 8,938,818 out of 42,031,845
Insertions: 483,373 out of 42,031,845
Identity over aligned base-pairs: 70.3%

Guppy

Uncovered: 676,515,753 out of 731,622,281
Covered: 55,106,528 out of 731,622,281

Uncovered: 7,714,381 out of 36,750,983
Matches: 20,591,508 out of 36,750,983
Mismatches: 7,848,993 out of 36,750,983
Insertions: 596,101 out of 36,750,983
Identity over aligned base-pairs: 70.9%

Block size distribution

Size range All 283,415 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
1,852
8.2 kb
10 bp - 100 bp
42,710
3.1 Mb
13,472
1.0 Mb
100 bp - 1 kb
232,885
56.7 Mb
72,222
22.4 Mb
1 kb - 10 kb
5,895
11.5 Mb
14,280
35.8 Mb
10 kb - 100 kb
73
997.7 kb
637
11.7 Mb
100 kb - 1 Mb
11
1.4 Mb