Human vs Long-tailed chinchilla LastZ results

Human (Homo sapiens, GRCh38) and Long-tailed chinchilla (Chinchilla lanigera, ChiLan1.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

ParameterHumanLong-tailed chinchilla
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,122,104 alignment blocks

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

Uncovered: 1,640,040,176 out of 3,096,649,726
Covered: 1,456,609,550 out of 3,096,649,726

Uncovered: 1,471,522 out of 35,387,657
Matches: 28,805,955 out of 35,387,657
Mismatches: 4,743,218 out of 35,387,657
Insertions: 366,962 out of 35,387,657
Identity over aligned base-pairs: 84.9%

Long-tailed chinchilla

Uncovered: 1,051,219,465 out of 2,390,868,971
Covered: 1,339,649,506 out of 2,390,868,971

Uncovered: 0 out of 11,351
Matches: 8,071 out of 11,351
Mismatches: 3,198 out of 11,351
Insertions: 82 out of 11,351
Identity over aligned base-pairs: 71.1%

Block size distribution

Size range All 2,122,104 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
58,155
389.0 kb
10 bp - 100 bp
248,634
12.3 Mb
16,271
1.1 Mb
100 bp - 1 kb
1,320,486
586.1 Mb
60,568
21.1 Mb
1 kb - 10 kb
494,364
920.8 Mb
17,510
42.9 Mb
10 kb - 100 kb
465
5.9 Mb
1,441
33.9 Mb
100 kb - 1 Mb
183
69.4 Mb
1 Mb - 10 Mb
198
734.4 Mb
10 Mb - 100 Mb
40
622.6 Mb