Human vs Pike-perch LastZ results

Human (Homo sapiens, GRCh38) and Pike-perch (Sander lucioperca, SLUC_FBN_1) 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
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

ParameterHumanPike-perch
Chunk size30,000,00010,100,000
Overlap0100,000
Group set size010,100,000
Masking options

Statistics over 267,667 alignment blocks

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

Uncovered: 3,061,408,345 out of 3,096,649,726
Covered: 35,241,381 out of 3,096,649,726

Uncovered: 14,479,832 out of 35,634,800
Matches: 14,514,077 out of 35,634,800
Mismatches: 6,280,538 out of 35,634,800
Insertions: 360,353 out of 35,634,800
Identity over aligned base-pairs: 68.6%

Pike-perch

Uncovered: 864,897,986 out of 900,477,767
Covered: 35,579,781 out of 900,477,767

Uncovered: 15,901,147 out of 40,902,852
Matches: 17,155,788 out of 40,902,852
Mismatches: 7,378,757 out of 40,902,852
Insertions: 467,160 out of 40,902,852
Identity over aligned base-pairs: 68.6%

Block size distribution

Size range All 267,667 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
1,900
8.2 kb
10 bp - 100 bp
56,298
4.2 Mb
19,637
1.5 Mb
100 bp - 1 kb
207,468
41.4 Mb
91,075
26.9 Mb
1 kb - 10 kb
1,992
3.2 Mb
9,639
18.8 Mb
10 kb - 100 kb
8
228.2 kb
94
1.7 Mb
100 kb - 1 Mb
1
100.2 kb
1
190.0 kb