Human vs Lamprey LastZ results

Human (Homo sapiens, GRCh38) and Lamprey (Petromyzon marinus, Pmarinus_7.0) 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
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

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

Statistics over 104,310 alignment blocks

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

Uncovered: 3,079,627,499 out of 3,096,649,726
Covered: 17,022,227 out of 3,096,649,726

Uncovered: 23,737,659 out of 35,228,305
Matches: 7,806,565 out of 35,228,305
Mismatches: 3,533,954 out of 35,228,305
Insertions: 150,127 out of 35,228,305
Identity over aligned base-pairs: 67.9%

Lamprey

Uncovered: 876,368,657 out of 885,550,958
Covered: 9,182,301 out of 885,550,958

Uncovered: 6,160,507 out of 12,541,363
Matches: 4,295,936 out of 12,541,363
Mismatches: 1,974,003 out of 12,541,363
Insertions: 110,917 out of 12,541,363
Identity over aligned base-pairs: 67.3%

Block size distribution

Size range All 104,310 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
356
1.5 kb
10 bp - 100 bp
23,464
1.8 Mb
10,563
814.6 kb
100 bp - 1 kb
79,876
14.7 Mb
46,596
12.3 Mb
1 kb - 10 kb
612
911.5 kb
2,520
4.2 Mb
10 kb - 100 kb
2
29.2 kb
6
134.1 kb