Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 43,125 | A→G | 21.6% | 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 | 21.6% | 8.1 ‑3.2 0.0 | 45 | pseudogene (184/420 nt) | caiB | crotonobetainyl‑CoA:carnitineCoA‑transferase; b0038_2 |
Reads supporting (aligned to +/- strand): new base (10/0): ref base (22/13): total (32/13) | ||||||||||
Fisher's exact test strand distribution p-value = 4.24e-02 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
TATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATGGTG > REL606/43092‑43159
|
tatTTAACGACAGCGCGTGCAAATTGCGGCGGggga > 1:2231674/1‑36 (MQ=255)
tatTTAACGACAGCGCGTGCAAATTGCGGCGGgaga < 1:7464751/36‑1 (MQ=255)
tatTTAACGACAGCGCGTGCAAATTGCGGCGGgaga > 1:1554687/1‑36 (MQ=255)
atTTAACGACAGCGCGTGCAAATTGCGGCGGGggag > 1:7745014/1‑36 (MQ=255)
tttAACGACAGCGCGTGCAAATTGCGGCGGGAGAGc < 1:3108512/36‑1 (MQ=255)
ttAACGACAGCGCGTGCAAATTGCGGCGGGggaggt > 1:1293633/1‑34 (MQ=38)
ttAACGACAGCGCGTGCAAATTGCGGCGGGGGAGCt > 1:1222267/1‑36 (MQ=255)
tAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCTg < 1:2339511/36‑1 (MQ=255)
tAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCTg < 1:405394/36‑1 (MQ=255)
aaCGACAGCGCGTGCAAATTGCGGCGGGGGAGGTgg > 1:671058/1‑35 (MQ=37)
aaCGACAGCGCGTGCAAATTGCGGCGGGGGAGGTgg > 1:2881510/1‑35 (MQ=37)
aCGACAGCGCGCGCAAATTGCGGCGGGGGAGCTgag > 1:7292718/1‑34 (MQ=37)
cGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCgg > 1:1131079/1‑36 (MQ=255)
cGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCgg > 1:5343062/1‑36 (MQ=255)
gACAGCGCGTGCAAATTGCGGCGGGGGAGCTGCggg > 1:2983786/1‑36 (MQ=255)
gACAGCGCGTGCAAATTGCGGCGGGGGAGCTGCggg > 1:5820515/1‑36 (MQ=255)
gACAGCGCGTGCAAATTGCGGCGGGGGAGCTGCggg > 1:1023867/1‑36 (MQ=255)
gACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCggg > 1:2904187/1‑36 (MQ=255)
gACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCggg < 1:2758375/36‑1 (MQ=255)
aCAGCGCGTGCAAATTGCGGCGGGGGAGGTGCgggg > 1:4200001/1‑35 (MQ=37)
aCAGCGCGTGCAAATTGCGGCGGGAGAGCTGAgggg > 1:4160406/1‑35 (MQ=255)
cAGCGCGTGCAAATTGCGGCGGGGGAGCTGCggggc > 1:7973139/1‑34 (MQ=38)
cAGCGCGTGCAAATTGCGGCGGGAGAGCTGCgggga > 1:7887207/1‑34 (MQ=255)
aGCGCGTGCAAATTGCGGCGGGGGAGCTGCggggcg > 1:2401345/1‑33 (MQ=37)
gcgTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTc < 1:6169597/36‑1 (MQ=255)
gTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCgg > 1:1348562/1‑36 (MQ=255)
gCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGtt < 1:3461791/36‑1 (MQ=255)
gCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGtt < 1:260091/36‑1 (MQ=255)
cAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTg < 1:4168681/36‑1 (MQ=255)
aaaTTGCGGCGGGAGAGCTGCGGGTAGTTTGGTTGa > 1:5353847/1‑36 (MQ=255)
aaaTTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGa < 1:5498077/36‑1 (MQ=255)
aaTTGCGGCGGGAGAGGTGCGGGTTGTTTGGGTgta > 1:169946/1‑34 (MQ=16)
aaTTGCGGCGGGAGAGCTGCGGGGAGTTCGGTTTaa > 1:928164/1‑36 (MQ=37)
aTTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGAAc < 1:7107837/36‑1 (MQ=255)
ttGCGGCGGGAGAGCTGCGGGGAGTTTGGTTTaaag > 1:1223776/1‑34 (MQ=21)
tGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGAAcgc > 1:401051/1‑36 (MQ=255)
gcggcgGGAGAGCTGCGGGTAGTTCGGTTGAAcgcg < 1:7035939/36‑1 (MQ=255)
cggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGa > 1:6099490/1‑36 (MQ=255)
cggcggGAGAGCTGCGGGGAGTTCGGTTGAACGCGa > 1:5664213/1‑36 (MQ=255)
ggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGaa > 1:5479470/1‑36 (MQ=255)
ggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGaa < 1:4009081/36‑1 (MQ=255)
ggcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGaa > 1:2513915/1‑36 (MQ=255)
gcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAAt > 1:4955443/1‑36 (MQ=255)
gcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAAt > 1:5391675/1‑36 (MQ=255)
gcggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAAt > 1:2250532/1‑36 (MQ=255)
cggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATg > 1:8146594/1‑36 (MQ=255)
cggGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATg < 1:5339606/36‑1 (MQ=255)
gggAGAGCTGCGGGTAGTTCGGTTGAACGCGAATgg > 1:2588378/1‑36 (MQ=255)
gggAGAGCTGCGGGTAGTTCGGTTGAACGCGAATgg > 1:1186747/1‑36 (MQ=255)
ggAGAGCTGCGGGTAGTTCGGGTGAACGCGAATGgt > 1:2905131/1‑36 (MQ=255)
gagaGCTGCGGGTAGTTCGGTTGAACGCGAATGgtg < 1:4729108/36‑1 (MQ=255)
gagaGCTGCGGGTAGTTCGGTTGAACGCGAATGgtg < 1:5683064/36‑1 (MQ=255)
gagaGCTGCGGGTAGTTCGGTTGAACGCGAATGgtg < 1:4730919/36‑1 (MQ=255)
|
TATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGCTGCGGGTAGTTCGGTTGAACGCGAATGGTG > REL606/43092‑43159
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 6 ≤
ATCG
< 19 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|