Human vs Hagfish LastZ results

Human (Homo sapiens, GRCh38) and Hagfish (Eptatretus burgeri, Eburgeri_3.2) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 93. 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
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

ParameterHumanHagfish
Chunk size30,000,00010,100,000
Overlap0100,000
Group set size010,100,000
Masking options{default_soft_masking => 1}{default_soft_masking => 1}

Statistics over 113,237 alignment blocks

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

Uncovered: 3,074,690,312 out of 3,096,649,726
Covered: 21,959,414 out of 3,096,649,726

Uncovered: 23,156,388 out of 35,522,272
Matches: 8,242,169 out of 35,522,272
Mismatches: 3,970,549 out of 35,522,272
Insertions: 153,166 out of 35,522,272
Identity over aligned base-pairs: 66.7%

Hagfish

Uncovered: 2,596,381,930 out of 2,608,383,542
Covered: 12,001,612 out of 2,608,383,542

Uncovered: 16,375,373 out of 23,977,300
Matches: 5,077,585 out of 23,977,300
Mismatches: 2,411,593 out of 23,977,300
Insertions: 112,749 out of 23,977,300
Identity over aligned base-pairs: 66.8%

Block size distribution

Size range All 113,237 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
583
2.6 kb
10 bp - 100 bp
23,794
1.8 Mb
12,524
980.4 kb
100 bp - 1 kb
87,372
16.1 Mb
61,235
14.9 Mb
1 kb - 10 kb
1,435
3.7 Mb
2,476
5.6 Mb
10 kb - 100 kb
52
628.1 kb
66
865.9 kb
100 kb - 1 Mb
1
105.4 kb
1
105.4 kb