Human vs Ferret LastZ results

Human (Homo sapiens, GRCh38) and Ferret (Mustela putorius furo, MusPutFur1.0) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 76. 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

ParameterHumanFerret
Chunk size30,000,00010,100,000
Overlap0100,000
Group set size010,100,000
Masking options{default_soft_masking => 1}

Statistics over 2,003,302 alignment blocks

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

Uncovered: 1,582,596,328 out of 3,096,649,726
Covered: 1,514,053,398 out of 3,096,649,726

Uncovered: 1,185,601 out of 35,228,305
Matches: 29,420,335 out of 35,228,305
Mismatches: 4,357,762 out of 35,228,305
Insertions: 264,607 out of 35,228,305
Identity over aligned base-pairs: 86.4%

Ferret

Uncovered: 979,992,275 out of 2,410,758,013
Covered: 1,430,765,738 out of 2,410,758,013

Uncovered: 1,370,312 out of 31,900,581
Matches: 26,690,001 out of 31,900,581
Mismatches: 3,675,633 out of 31,900,581
Insertions: 164,635 out of 31,900,581
Identity over aligned base-pairs: 87.4%

Block size distribution

Size range All 2,003,302 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
49,260
319.5 kb
10 bp - 100 bp
183,782
10.1 Mb
16,851
1.2 Mb
100 bp - 1 kb
1,239,330
560.8 Mb
57,886
21.3 Mb
1 kb - 10 kb
530,241
1.0 Gb
18,109
44.0 Mb
10 kb - 100 kb
689
8.8 Mb
1,566
38.9 Mb
100 kb - 1 Mb
285
108.5 Mb
1 Mb - 10 Mb
313
1.1 Gb
10 Mb - 100 Mb
21
316.3 Mb