Human vs C.intestinalis LastZ results

Human (Homo sapiens, GRCh38) and C.intestinalis (Ciona intestinalis, KH) 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

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

Statistics over 64,090 alignment blocks

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

Uncovered: 3,087,903,268 out of 3,096,649,726
Covered: 8,746,458 out of 3,096,649,726

Uncovered: 30,013,881 out of 35,228,305
Matches: 3,446,027 out of 35,228,305
Mismatches: 1,724,260 out of 35,228,305
Insertions: 44,137 out of 35,228,305
Identity over aligned base-pairs: 66.1%

C.intestinalis

Uncovered: 110,976,375 out of 115,226,814
Covered: 4,250,439 out of 115,226,814

Uncovered: 12,950,169 out of 16,105,584
Matches: 2,058,593 out of 16,105,584
Mismatches: 1,057,686 out of 16,105,584
Insertions: 39,136 out of 16,105,584
Identity over aligned base-pairs: 65.2%

Block size distribution

Size range All 64,090 alignment blocks Blocks grouped in nets
# blocks Total size (incl. gaps) # nets Total size (incl. gaps)
1 bp - 10 bp
265
1.3 kb
10 bp - 100 bp
18,488
1.4 Mb
2,378
180.9 kb
100 bp - 1 kb
45,271
7.9 Mb
15,742
5.7 Mb
1 kb - 10 kb
198
247.9 kb
2,171
3.5 Mb
10 kb - 100 kb
4
79.2 kb