Human vs Domestic yak LastZ results

Human (Homo sapiens, GRCh38) and Domestic yak (Bos grunniens, LU_Bosgru_v3.0) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 99. 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

ParameterHumanDomestic yak
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,336,359 alignment blocks

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

Uncovered: 1,696,244,926 out of 3,096,649,726
Covered: 1,400,404,800 out of 3,096,649,726

Uncovered: 1,270,322 out of 35,602,741
Matches: 29,458,677 out of 35,602,741
Mismatches: 4,546,046 out of 35,602,741
Insertions: 327,696 out of 35,602,741
Identity over aligned base-pairs: 85.8%

Domestic yak

Uncovered: 1,520,827,624 out of 2,832,776,395
Covered: 1,311,948,771 out of 2,832,776,395

Uncovered: 4,509,744 out of 34,628,786
Matches: 26,170,946 out of 34,628,786
Mismatches: 3,763,512 out of 34,628,786
Insertions: 184,584 out of 34,628,786
Identity over aligned base-pairs: 86.9%

Block size distribution

Size range All 2,336,359 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
75,996
507.7 kb
10 bp - 100 bp
283,034
14.1 Mb
17,035
1.2 Mb
100 bp - 1 kb
1,521,327
675.0 Mb
57,534
20.6 Mb
1 kb - 10 kb
455,787
780.9 Mb
16,942
40.9 Mb
10 kb - 100 kb
210
3.2 Mb
1,405
32.5 Mb
100 kb - 1 Mb
5
569.7 kb
137
46.1 Mb
1 Mb - 10 Mb
113
461.4 Mb
10 Mb - 100 Mb
43
871.6 Mb