Human vs Coquerel's sifaka LastZ results

Human (Homo sapiens, GRCh38) and Coquerel's sifaka (Propithecus coquereli, Pcoq_1.0) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 91. 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)
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

ParameterHumanCoquerel's sifaka
Chunk size30,000,00010,100,000
Overlap0100,000
Group set size010,100,000
Masking options{default_soft_masking => 1}{default_soft_masking => 1}

Statistics over 2,025,030 alignment blocks

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

Uncovered: 1,398,806,143 out of 3,096,649,726
Covered: 1,697,843,583 out of 3,096,649,726

Uncovered: 990,624 out of 35,444,608
Matches: 30,918,429 out of 35,444,608
Mismatches: 3,330,321 out of 35,444,608
Insertions: 205,234 out of 35,444,608
Identity over aligned base-pairs: 89.7%

Coquerel's sifaka

Uncovered: 1,217,014,756 out of 2,798,152,141
Covered: 1,581,137,385 out of 2,798,152,141

Uncovered: 593,889 out of 29,815,600
Matches: 26,637,533 out of 29,815,600
Mismatches: 2,499,736 out of 29,815,600
Insertions: 84,442 out of 29,815,600
Identity over aligned base-pairs: 91.2%

Block size distribution

Size range All 2,025,030 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
70,667
479.3 kb
10 bp - 100 bp
286,765
13.8 Mb
33,018
2.2 Mb
100 bp - 1 kb
1,095,494
479.7 Mb
136,933
51.4 Mb
1 kb - 10 kb
569,914
1.2 Gb
47,705
95.4 Mb
10 kb - 100 kb
2,190
28.2 Mb
1,856
47.5 Mb
100 kb - 1 Mb
578
229.6 Mb
1 Mb - 10 Mb
417
1.2 Gb
10 Mb - 100 Mb
8
109.2 Mb