Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 1,080,000 | C→G | 31.1% | K22N (AAG→AAC) | cbpA ← | curved DNA‑binding protein, DnaJ homologue that functions as a co‑chaperone of DnaK |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 1,080,000 | 0 | C→G | 31.1% | 2.3 ‑4.5 0.0 | 20 | 66 | cbpA | curved DNA‑binding protein, DnaJ homologue that functions as a co‑chaperone of DnaK |
Reads supporting (aligned to +/- strand): new base (6/0): ref base (4/9): total (11/9) | ||||||||||
Fisher's exact test strand distribution p-value = 1.08e-02 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
GGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGC > REL606/1079965‑1080034
|
ggATGGTATTTGCGGGCAAGTCGACGATAGGCGGTc > 1:5871450/1‑36 (MQ=255)
aTGGTATTTGCGGGCAAGTCGACGATAGGCGGTGtt > 1:3321513/1‑36 (MQ=255)
gTATTTGCGGGCAAGTCGACGATAGGCGGTTttttt > 1:3581141/1‑33 (MQ=37)
gTATTTGCGGGCAAGTCGACGATAGGCGGGGTTGAt > 1:15222/1‑36 (MQ=37)
tATTTGCGGGCAAGTCGACGATAGGCGGTCTTGAtt < 1:861393/36‑1 (MQ=255)
tATTTGCGGGCAAGTCGACGATAGGCGGTCTTGAtt < 1:2913646/36‑1 (MQ=255)
tttGCGGGCAAGTCGACGATAGGCGGTGTTGCtttt > 1:4256885/1‑34 (MQ=25)
gCGGGCAAGTCGACGATCGGCGGTCTTGATTGTCtt < 1:4238418/36‑1 (MQ=255)
gCGGGCAAGTCGACGATAGGCGGGGTTGGTTGTTtt > 1:3547137/1‑36 (MQ=19)
cgcgCAAGTCGACGATAGGCGGTCTTGATTGTCTTg < 1:3602957/33‑1 (MQ=255)
cGGGCAAGTCGACGATAGTCGGTCTTGATTGTCTTg < 1:2168356/36‑1 (MQ=255)
cGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTg < 1:4736438/36‑1 (MQ=255)
ggCAAGTCGACGATAGGCGGTCTTGATTGTCTTgag < 1:4204998/36‑1 (MQ=255)
gCAAGTCGACGATAGGCGGTCTTGATTGTCTTgaga < 1:1912255/36‑1 (MQ=255)
gTCGACGATAAGCGGGGTTATTTTTCTTGAGATtgt > 1:1776614/1‑36 (MQ=2)
tCGACGATAGGCGGAGTTGATTGTTTTTTGATcgtg > 1:1424925/1‑35 (MQ=12)
gacgaTAGGCGGTCTTGTTTTTCTTTAGATcgtccg > 1:2305144/1‑36 (MQ=37)
gacgaTAGGCGGTCTTGATTTTCTTTTGATcgtccg > 1:5039936/1‑36 (MQ=25)
ggCGGTCTTGATTGTCTTGAGATCGTCCGTCGGttt < 1:4052650/36‑1 (MQ=255)
ggTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCAc > 1:2359877/1‑36 (MQ=255)
tCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGc > 1:5863834/1‑36 (MQ=255)
|
GGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTTGATTGTCTTGAGATCGTCCGTCGGTTTCACGC > REL606/1079965‑1080034
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 5 ≤
ATCG
< 12 ≤
ATCG
< 26 ≤
ATCG
< 37 ≤
ATCG
|