Human vs Chicken LastZ results

Human (Homo sapiens, GRCh38) and Chicken (Gallus gallus, bGalGal1.mat.broiler.GRCg7b) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 109. 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

ParameterHumanChicken
Chunk size30,000,00010,100,000
Overlap0100,000
Group set size10,100,000
Masking options

Statistics over 281,940 alignment blocks

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

Uncovered: 3,021,988,993 out of 3,099,750,718
Covered: 77,761,725 out of 3,099,750,718

Uncovered: 12,676,172 out of 38,571,754
Matches: 19,113,951 out of 38,571,754
Mismatches: 6,381,214 out of 38,571,754
Insertions: 400,417 out of 38,571,754
Identity over aligned base-pairs: 73.8%

Chicken

Uncovered: 989,127,915 out of 1,053,332,251
Covered: 64,204,336 out of 1,053,332,251

Uncovered: 8,569,932 out of 31,304,036
Matches: 16,937,059 out of 31,304,036
Mismatches: 5,507,626 out of 31,304,036
Insertions: 289,419 out of 31,304,036
Identity over aligned base-pairs: 74.5%

Block size distribution

Size range All 281,940 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
809
4.3 kb
10 bp - 100 bp
34,917
2.6 Mb
8,676
654.7 kb
100 bp - 1 kb
236,683
66.2 Mb
81,676
28.6 Mb
1 kb - 10 kb
9,523
14.7 Mb
14,487
34.0 Mb
10 kb - 100 kb
7
86.3 kb
697
15.5 Mb
100 kb - 1 Mb
1
100.2 kb
29
4.9 Mb