Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 43,125 | A→G | 17.5% | pseudogene (184/420 nt) | caiB ← | crotonobetainyl‑CoA:carnitineCoA‑transferase; b0038_2 |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 43,125 | 0 | A→G | 17.5% | 5.1 ‑2.5 ‑0.1 | 45 | pseudogene (184/420 nt) | caiB | crotonobetainyl‑CoA:carnitineCoA‑transferase; b0038_2 |
Reads supporting (aligned to +/- strand): new base (8/0): ref base (23/13): total (32/13) | ||||||||||
Fisher's exact test strand distribution p-value = 8.20e-02 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.77e-01 |
TATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATGGTGT > REL606/43092‑43160
|
tatTTAACGACAGCGCGTGCAAATTGCGGCGGgaga > 1:7840020/1‑36 (MQ=255)
tatTTAACGACAGCGCGTGCAAATTGCGGCGGgaga < 1:5442039/36‑1 (MQ=255)
tatTTAACGACAGCGCGTGCAAATTGCGGCGGgaga < 1:6978909/36‑1 (MQ=255)
atTTAACGACAGCGCGTGCAAATTGCGGCGGGggag > 1:3337430/1‑36 (MQ=255)
tttAACGACAGCGCGTGCAAATTGCGGCGGGGGAGc > 1:4029417/1‑36 (MQ=255)
ttAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCt > 1:6834904/1‑36 (MQ=255)
tAACGACAGCGCGTGCAAATTGCGGCGGGGGAGCTg > 1:5222489/1‑36 (MQ=255)
aaCGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGc < 1:4938438/36‑1 (MQ=255)
aCGACAGCGCGTGCAAATTGCGGCGGGGGGGCTGCg > 1:21011/1‑36 (MQ=38)
aCGACAGCGCGTGCAAATTGCGGCGGGGGAGCTGCg > 1:5638474/1‑36 (MQ=255)
aCGACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCg > 1:6070470/1‑36 (MQ=255)
aCGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGct < 1:4975068/36‑2 (MQ=255)
cGACAGCGCGTGCAAATTGCGGCGGGTGAGCCGCgg > 1:1813915/1‑36 (MQ=38)
cGACAGCGCGTGCAAATTGCGGCGGGGGAGGCGCgg > 1:3485333/1‑36 (MQ=21)
gACAGCGCGTGCAAATTGCGGCGGTGGAGCCTCggg > 1:2659819/1‑36 (MQ=19)
gACAGCGCGTGCAAATTGCGGCGGGGGAGGTGCggg > 1:2732905/1‑36 (MQ=37)
gACAGCGCGTGCAAATTGCGGCGGGGGAGGTGCggg > 1:3290011/1‑36 (MQ=37)
gACAGCGCGTGCAAATTGCGGCGGGAGAGGTGAggg > 1:5432402/1‑36 (MQ=37)
gACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCggg > 1:3374722/1‑36 (MQ=255)
cAGCGCGTGCAAATTGCGGCGGGAGAGCTGCggggg > 1:3871197/1‑34 (MQ=255)
aGCGCGTGCAAATTGCGGCGGGAGAGCTGCGGGTAg > 1:6231787/1‑36 (MQ=255)
cgcgTGCAAATTGCGGCGGGAGAGCTGCGGGTAGtt > 1:577408/1‑36 (MQ=255)
gcgTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTc < 1:4344905/36‑1 (MQ=255)
gcgTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTc < 1:2814275/36‑1 (MQ=255)
gcgTGCAAATTGCGGCGGGAGAGCTGCGGGGGGTTc > 1:164532/1‑36 (MQ=37)
gTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTTgg > 1:6337435/1‑36 (MQ=255)
tcAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGtt < 1:3821515/35‑1 (MQ=255)
gCAAATTGCGGCGGGAGAGCTGGGGGTAGTTTgggt > 1:3862627/1‑34 (MQ=25)
gCAAATTGCGGCGGGAGAGCTGGCGGGAGGTCGGtt > 1:5264637/1‑36 (MQ=16)
gCAAATTGCGGCGGGAGAGCTGCGGGGAGTTCGGtt > 1:3366802/1‑36 (MQ=255)
aaaTTGCGGCGGGAGAGCTGCGGATAGTTCGGTTGa < 1:204931/36‑1 (MQ=255)
aTTGCGGCGGGAGAGCTGCGGGTAGTTTGGTTGaag > 1:6229322/1‑35 (MQ=255)
aTTGCGGCGGGACAGCTGCGGGTGGTTCGGGTGatc > 1:1134957/1‑34 (MQ=25)
ttGCGGCGGGAGAGCTGCCGGTAGGTCGGGTGAAcg > 1:7383682/1‑36 (MQ=25)
gcggcgGGAGAGCTGCGGGTAGGTTGGTTGAAcgcg > 1:5444135/1‑36 (MQ=37)
cggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGa > 1:7776457/1‑36 (MQ=255)
cggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGa < 1:3623427/36‑1 (MQ=255)
cggcggGAGAGCTGCAGGTAGTTCGGTTGAACGCGa < 1:3151600/36‑1 (MQ=255)
ggggggAGTGCTGCGGGTAGTTCGGTTGAACGCGaa < 1:256805/33‑1 (MQ=38)
ggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGaa > 1:2959747/1‑36 (MQ=255)
ggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGaa > 1:2487032/1‑36 (MQ=255)
gcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAAt < 1:5155652/36‑1 (MQ=255)
gcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAAt < 1:2314663/36‑1 (MQ=255)
gcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAAt > 1:1079708/1‑36 (MQ=255)
cggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATg < 1:5956666/36‑1 (MQ=255)
gggAGAGCTGCGGGTAGTTCGGTTGAACGCGAATgg > 1:5262359/1‑36 (MQ=255)
gggAGAGCTGCGGGTAGTTCGGTTGAACGCGAATgg > 1:4607032/1‑36 (MQ=255)
ggAGAGCTGCGGGTAGTTCGGTTTAACGCGAATGgt > 1:1232837/1‑36 (MQ=255)
ggAGAGCTGCGGGTAGTTCGGTTGAACGCGAATGgt < 1:7872151/36‑1 (MQ=255)
agagCTGCGGGTAGTTCGGTTGAACGCGAATGgtgt < 1:2298587/36‑1 (MQ=255)
agagCTGCGGGTAGTTCGGTTGAACGCGAATGgtgt > 1:813008/1‑36 (MQ=255)
|
TATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATGGTGT > REL606/43092‑43160
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 6 ≤
ATCG
< 16 ≤
ATCG
< 29 ≤
ATCG
< 37 ≤
ATCG
|