Zebrafish vs Lamprey LastZ results

Zebrafish (Danio rerio, GRCz11) and Lamprey (Petromyzon marinus, Pmarinus_7.0) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 92. 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

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

Statistics over 118,596 alignment blocks

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

Uncovered: 1,353,521,715 out of 1,373,471,384
Covered: 19,949,669 out of 1,373,471,384

Uncovered: 28,000,241 out of 42,031,845
Matches: 9,435,959 out of 42,031,845
Mismatches: 4,437,387 out of 42,031,845
Insertions: 158,258 out of 42,031,845
Identity over aligned base-pairs: 67.2%

Lamprey

Uncovered: 872,565,139 out of 885,550,958
Covered: 12,985,819 out of 885,550,958

Uncovered: 4,838,883 out of 12,541,363
Matches: 5,182,668 out of 12,541,363
Mismatches: 2,388,654 out of 12,541,363
Insertions: 131,158 out of 12,541,363
Identity over aligned base-pairs: 67.3%

Block size distribution

Size range All 118,596 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
552
2.2 kb
10 bp - 100 bp
24,199
1.9 Mb
11,505
886.9 kb
100 bp - 1 kb
93,168
17.5 Mb
51,864
13.5 Mb
1 kb - 10 kb
629
1.1 Mb
3,103
5.8 Mb
10 kb - 100 kb
39
1.7 Mb
44
1.8 Mb
100 kb - 1 Mb
9
1.2 Mb
10
1.4 Mb