Human vs Degu LastZ results

Human (Homo sapiens, GRCh38) and Degu (Octodon degus, OctDeg1.0) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 90. 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

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

Statistics over 2,024,428 alignment blocks

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

Uncovered: 1,817,075,946 out of 3,096,649,726
Covered: 1,279,573,780 out of 3,096,649,726

Uncovered: 1,357,156 out of 35,444,608
Matches: 28,620,848 out of 35,444,608
Mismatches: 5,064,016 out of 35,444,608
Insertions: 402,588 out of 35,444,608
Identity over aligned base-pairs: 84.0%

Degu

Uncovered: 1,820,280,886 out of 2,995,872,505
Covered: 1,175,591,619 out of 2,995,872,505

Uncovered: 2,236,226 out of 29,121,407
Matches: 23,240,790 out of 29,121,407
Mismatches: 3,535,620 out of 29,121,407
Insertions: 108,771 out of 29,121,407
Identity over aligned base-pairs: 86.4%

Block size distribution

Size range All 2,024,428 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
48,827
327.4 kb
10 bp - 100 bp
228,009
11.7 Mb
16,653
1.1 Mb
100 bp - 1 kb
1,328,142
590.5 Mb
56,555
20.8 Mb
1 kb - 10 kb
419,128
737.0 Mb
19,707
46.6 Mb
10 kb - 100 kb
322
4.1 Mb
1,348
32.0 Mb
100 kb - 1 Mb
309
126.5 Mb
1 Mb - 10 Mb
319
1.0 Gb
10 Mb - 100 Mb
7
93.8 Mb