Human vs Lumpfish LastZ results

Human (Homo sapiens, GRCh38) and Lumpfish (Cyclopterus lumpus, fCycLum1.pri) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 101. 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

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

Statistics over 254,034 alignment blocks

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

Uncovered: 3,065,976,316 out of 3,099,750,718
Covered: 33,774,402 out of 3,099,750,718

Uncovered: 17,620,997 out of 38,571,754
Matches: 14,351,792 out of 38,571,754
Mismatches: 6,223,695 out of 38,571,754
Insertions: 375,270 out of 38,571,754
Identity over aligned base-pairs: 68.5%

Lumpfish

Uncovered: 540,256,093 out of 572,885,595
Covered: 32,629,502 out of 572,885,595

Uncovered: 13,112,326 out of 36,193,391
Matches: 15,877,180 out of 36,193,391
Mismatches: 6,766,506 out of 36,193,391
Insertions: 437,379 out of 36,193,391
Identity over aligned base-pairs: 68.8%

Block size distribution

Size range All 254,034 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
1,915
8.5 kb
10 bp - 100 bp
52,446
3.9 Mb
14,459
1.1 Mb
100 bp - 1 kb
197,968
39.3 Mb
77,983
23.9 Mb
1 kb - 10 kb
1,695
2.6 Mb
9,869
19.3 Mb
10 kb - 100 kb
10
439.3 kb
96
1.9 Mb
100 kb - 1 Mb
1
103.9 kb