Human vs Sailfin molly LastZ results

Human (Homo sapiens, GRCh38) and Sailfin molly (Poecilia latipinna, P_latipinna-1.0) were aligned using the LastZ alignment algorithm (LastZ) in Ensembl release 94. 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

ParameterHumanSailfin molly
Chunk size30,000,00010,100,000
Overlap0100,000
Group set size010,100,000
Masking options

Statistics over 264,805 alignment blocks

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

Uncovered: 3,062,281,946 out of 3,096,649,726
Covered: 34,367,780 out of 3,096,649,726

Uncovered: 15,123,479 out of 35,619,351
Matches: 14,006,239 out of 35,619,351
Mismatches: 6,131,516 out of 35,619,351
Insertions: 358,117 out of 35,619,351
Identity over aligned base-pairs: 68.3%

Sailfin molly

Uncovered: 781,716,557 out of 815,144,743
Covered: 33,428,186 out of 815,144,743

Uncovered: 13,018,983 out of 35,799,401
Matches: 15,654,877 out of 35,799,401
Mismatches: 6,706,819 out of 35,799,401
Insertions: 418,722 out of 35,799,401
Identity over aligned base-pairs: 68.7%

Block size distribution

Size range All 264,805 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
1,570
6.8 kb
10 bp - 100 bp
53,754
4.1 Mb
16,267
1.2 Mb
100 bp - 1 kb
207,334
41.3 Mb
87,278
26.3 Mb
1 kb - 10 kb
2,140
4.2 Mb
10,410
20.8 Mb
10 kb - 100 kb
7
174.0 kb
75
1.1 Mb
100 kb - 1 Mb
2
212.7 kb