Cow vs Vaquita LastZ results

Cow (Bos taurus, ARS-UCD1.3) and Vaquita (Phocoena sinus, mPhoSin1.pri) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 111. Cow 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

ParameterCowVaquita
Chunk size10,000,00010,100,000
Overlap0100,000
Group set size10,100,00010,100,000
Masking options

Statistics over 2,062,320 alignment blocks

Genome coverage (bp) Coding exon coverage (bp)
Cow

Uncovered: 1,064,359,625 out of 2,711,209,831
Covered: 1,646,850,206 out of 2,711,209,831

Uncovered: 3,268,776 out of 37,648,360
Matches: 31,263,571 out of 37,648,360
Mismatches: 2,945,201 out of 37,648,360
Insertions: 170,812 out of 37,648,360
Identity over aligned base-pairs: 90.9%

Vaquita

Uncovered: 760,066,117 out of 2,371,524,154
Covered: 1,611,458,037 out of 2,371,524,154

Uncovered: 899,287 out of 31,056,201
Matches: 27,750,000 out of 31,056,201
Mismatches: 2,296,252 out of 31,056,201
Insertions: 110,662 out of 31,056,201
Identity over aligned base-pairs: 92.0%

Block size distribution

Size range All 2,062,320 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
72,199
461.1 kb
10 bp - 100 bp
186,815
9.7 Mb
20,960
1.4 Mb
100 bp - 1 kb
1,227,899
557.0 Mb
52,091
16.5 Mb
1 kb - 10 kb
573,970
1.1 Gb
9,338
23.6 Mb
10 kb - 100 kb
1,433
20.1 Mb
824
18.9 Mb
100 kb - 1 Mb
4
458.6 kb
55
20.7 Mb
1 Mb - 10 Mb
44
188.1 Mb
10 Mb - 100 Mb
55
1.4 Gb