Human vs Greater horseshoe bat LastZ results

Human (Homo sapiens, GRCh38) and Greater horseshoe bat (Rhinolophus ferrumequinum, mRhiFer1_v1.p) 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
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

ParameterHumanGreater horseshoe bat
Chunk size30,000,00010,100,000
Overlap0100,000
Group set size10,100,000
Masking options{default_soft_masking => 1}{default_soft_masking => 1}

Statistics over 1,850,961 alignment blocks

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

Uncovered: 1,569,182,869 out of 3,099,750,718
Covered: 1,530,567,849 out of 3,099,750,718

Uncovered: 1,635,669 out of 38,571,754
Matches: 31,836,426 out of 38,571,754
Mismatches: 4,672,511 out of 38,571,754
Insertions: 427,148 out of 38,571,754
Identity over aligned base-pairs: 86.2%

Greater horseshoe bat

Uncovered: 640,908,923 out of 2,075,768,562
Covered: 1,434,859,639 out of 2,075,768,562

Uncovered: 950,907 out of 32,593,495
Matches: 27,928,785 out of 32,593,495
Mismatches: 3,569,440 out of 32,593,495
Insertions: 144,363 out of 32,593,495
Identity over aligned base-pairs: 88.3%

Block size distribution

Size range All 1,850,961 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
61,886
421.8 kb
10 bp - 100 bp
233,545
10.9 Mb
18,989
1.3 Mb
100 bp - 1 kb
1,002,172
460.9 Mb
57,977
20.9 Mb
1 kb - 10 kb
551,951
1.1 Gb
18,367
45.7 Mb
10 kb - 100 kb
1,396
19.6 Mb
1,750
46.4 Mb
100 kb - 1 Mb
11
1.3 Mb
233
60.3 Mb
1 Mb - 10 Mb
72
273.0 Mb
10 Mb - 100 Mb
41
1.2 Gb