Human vs Tetraodon LastZ results

Human (Homo sapiens, GRCh38) and Tetraodon (Tetraodon nigroviridis, TETRAODON8) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 76. Human 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
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

ParameterHumanTetraodon
Chunk size30,000,00010,100,000
Overlap0100,000
Group set size10,100,000
Masking options

Statistics over 250,608 alignment blocks

Genome coverage (bp) Coding exon coverage (bp)
Human

Uncovered: 3,062,734,735 out of 3,099,750,718
Covered: 37,015,983 out of 3,099,750,718

Uncovered: 18,056,949 out of 38,571,754
Matches: 13,860,202 out of 38,571,754
Mismatches: 6,255,322 out of 38,571,754
Insertions: 399,281 out of 38,571,754
Identity over aligned base-pairs: 67.6%

Tetraodon

Uncovered: 325,014,075 out of 358,618,246
Covered: 33,604,171 out of 358,618,246

Uncovered: 8,382,245 out of 30,080,167
Matches: 14,914,836 out of 30,080,167
Mismatches: 6,386,702 out of 30,080,167
Insertions: 396,384 out of 30,080,167
Identity over aligned base-pairs: 68.7%

Block size distribution

Size range All 250,608 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
1,670
7.8 kb
10 bp - 100 bp
51,551
3.8 Mb
14,643
1.1 Mb
100 bp - 1 kb
194,101
41.0 Mb
74,873
23.9 Mb
1 kb - 10 kb
3,246
5.0 Mb
11,883
23.4 Mb
10 kb - 100 kb
37
1.0 Mb
116
2.2 Mb
100 kb - 1 Mb
3
374.7 kb
5
645.2 kb