Human vs American black bear LastZ results

Human (Homo sapiens, GRCh38) and American black bear (Ursus americanus, ASM334442v1) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 95. 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

ParameterHumanAmerican black bear
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,200,350 alignment blocks

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

Uncovered: 1,446,760,370 out of 3,096,649,726
Covered: 1,649,889,356 out of 3,096,649,726

Uncovered: 671,955 out of 35,570,071
Matches: 30,230,231 out of 35,570,071
Mismatches: 4,410,398 out of 35,570,071
Insertions: 257,487 out of 35,570,071
Identity over aligned base-pairs: 86.6%

American black bear

Uncovered: 1,049,639,750 out of 2,588,393,916
Covered: 1,538,754,166 out of 2,588,393,916

Uncovered: 1,900,240 out of 30,870,688
Matches: 25,431,909 out of 30,870,688
Mismatches: 3,354,303 out of 30,870,688
Insertions: 184,236 out of 30,870,688
Identity over aligned base-pairs: 87.8%

Block size distribution

Size range All 2,200,350 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
78,865
525.7 kb
10 bp - 100 bp
272,098
12.9 Mb
18,291
1.2 Mb
100 bp - 1 kb
1,272,153
576.3 Mb
68,925
26.5 Mb
1 kb - 10 kb
576,541
1.1 Gb
35,244
104.8 Mb
10 kb - 100 kb
691
9.6 Mb
12,984
513.9 Mb
100 kb - 1 Mb
2
255.1 kb
5,190
1.1 Gb
1 Mb - 10 Mb
2
2.9 Mb