Human vs Sloth LastZ results

Human (Homo sapiens, GRCh38) and Sloth (Choloepus hoffmanni, choHof1) 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
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

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

Statistics over 1,484,560 alignment blocks

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

Uncovered: 2,075,082,329 out of 3,099,750,718
Covered: 1,024,668,389 out of 3,099,750,718

Uncovered: 12,023,347 out of 38,571,754
Matches: 22,136,645 out of 38,571,754
Mismatches: 4,082,479 out of 38,571,754
Insertions: 329,283 out of 38,571,754
Identity over aligned base-pairs: 83.4%

Sloth

Uncovered: 1,513,167,168 out of 2,467,493,193
Covered: 954,326,025 out of 2,467,493,193

Uncovered: 3,992,807 out of 19,729,254
Matches: 13,949,308 out of 19,729,254
Mismatches: 1,750,854 out of 19,729,254
Insertions: 36,285 out of 19,729,254
Identity over aligned base-pairs: 88.6%

Block size distribution

Size range All 1,484,560 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
10,751
68.3 kb
10 bp - 100 bp
114,410
7.2 Mb
23,210
1.6 Mb
100 bp - 1 kb
996,630
463.0 Mb
160,774
75.6 Mb
1 kb - 10 kb
362,729
605.3 Mb
148,966
474.1 Mb
10 kb - 100 kb
40
458.6 kb
23,206
475.3 Mb
100 kb - 1 Mb
345
49.3 Mb