Human vs Brown trout LastZ results

Human (Homo sapiens, GRCh38) and Brown trout (Salmo trutta, fSalTru1.1) 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

ParameterHumanBrown trout
Chunk size30,000,00010,100,000
Overlap0100,000
Group set size010,100,000
Masking options

Statistics over 440,439 alignment blocks

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

Uncovered: 3,058,186,365 out of 3,096,649,726
Covered: 38,463,361 out of 3,096,649,726

Uncovered: 13,529,518 out of 35,619,351
Matches: 15,100,906 out of 35,619,351
Mismatches: 6,571,748 out of 35,619,351
Insertions: 417,179 out of 35,619,351
Identity over aligned base-pairs: 68.4%

Brown trout

Uncovered: 2,304,799,695 out of 2,371,880,186
Covered: 67,080,491 out of 2,371,880,186

Uncovered: 33,838,746 out of 79,685,618
Matches: 31,607,231 out of 79,685,618
Mismatches: 13,261,253 out of 79,685,618
Insertions: 978,388 out of 79,685,618
Identity over aligned base-pairs: 68.9%

Block size distribution

Size range All 440,439 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
3,277
14.6 kb
10 bp - 100 bp
86,326
6.4 Mb
28,775
2.2 Mb
100 bp - 1 kb
347,286
70.9 Mb
151,777
45.0 Mb
1 kb - 10 kb
3,517
5.5 Mb
16,778
32.9 Mb
10 kb - 100 kb
27
1.2 Mb
171
3.3 Mb
100 kb - 1 Mb
6
717.0 kb
9
1.3 Mb