Zebrafish vs Burton's mouthbrooder LastZ results

Zebrafish (Danio rerio, GRCz11) and Burton's mouthbrooder (Haplochromis burtoni, AstBur1.0) 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

ParameterZebrafishBurton's mouthbrooder
Chunk size10,000,00010,100,000
Overlap0100,000
Group set size010,100,000
Masking options

Statistics over 299,231 alignment blocks

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

Uncovered: 1,309,144,489 out of 1,373,471,384
Covered: 64,326,895 out of 1,373,471,384

Uncovered: 9,850,544 out of 42,031,845
Matches: 22,672,340 out of 42,031,845
Mismatches: 9,050,934 out of 42,031,845
Insertions: 458,027 out of 42,031,845
Identity over aligned base-pairs: 70.5%

Burton's mouthbrooder

Uncovered: 770,495,598 out of 831,411,547
Covered: 60,915,949 out of 831,411,547

Uncovered: 7,349,568 out of 37,764,289
Matches: 21,638,697 out of 37,764,289
Mismatches: 8,114,233 out of 37,764,289
Insertions: 661,791 out of 37,764,289
Identity over aligned base-pairs: 71.1%

Block size distribution

Size range All 299,231 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
1,866
8.3 kb
10 bp - 100 bp
44,762
3.2 Mb
13,942
1.0 Mb
100 bp - 1 kb
246,214
61.3 Mb
73,983
23.1 Mb
1 kb - 10 kb
6,224
10.7 Mb
14,551
35.9 Mb
10 kb - 100 kb
165
3.2 Mb
811
16.7 Mb
100 kb - 1 Mb
11
1.7 Mb