Human vs Common canary LastZ results

Human (Homo sapiens, GRCh38) and Common canary (Serinus canaria, SCA1) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 96. 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

ParameterHumanCommon canary
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 292,978 alignment blocks

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

Uncovered: 3,015,766,328 out of 3,096,649,726
Covered: 80,883,398 out of 3,096,649,726

Uncovered: 9,729,110 out of 35,604,148
Matches: 18,899,822 out of 35,604,148
Mismatches: 6,635,622 out of 35,604,148
Insertions: 339,594 out of 35,604,148
Identity over aligned base-pairs: 73.0%

Common canary

Uncovered: 1,085,983,697 out of 1,152,083,301
Covered: 66,099,604 out of 1,152,083,301

Uncovered: 4,592,690 out of 25,036,602
Matches: 15,267,080 out of 25,036,602
Mismatches: 4,941,299 out of 25,036,602
Insertions: 235,533 out of 25,036,602
Identity over aligned base-pairs: 74.7%

Block size distribution

Size range All 292,978 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
653
3.4 kb
10 bp - 100 bp
41,548
3.1 Mb
14,000
1.0 Mb
100 bp - 1 kb
241,570
66.7 Mb
90,000
30.7 Mb
1 kb - 10 kb
9,189
14.3 Mb
14,072
33.2 Mb
10 kb - 100 kb
17
318.9 kb
661
14.8 Mb
100 kb - 1 Mb
1
100.2 kb
26
4.8 Mb