Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 393,086 | A→G | 15.0% | pseudogene (19/345 nt) | yajB ← | putative glycoprotein; b0404_1 |
Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 393,086 | 0 | A→G | 15.0% | 4.1 ‑3.9 ‑0.0 | 37 | pseudogene (19/345 nt) | yajB | putative glycoprotein; b0404_1 |
Reads supporting (aligned to +/- strand): new base (6/0): ref base (13/18): total (19/18) | ||||||||||
Fisher's exact test strand distribution p-value = 1.97e-02 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.59e-01 |
CGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAG > REL606/393051‑393121
|
cgtaaaGTCCGGCTACAGCTGCGACCAGTGGCGGga < 1:7862558/33‑1 (MQ=37)
cGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGga < 1:1725142/36‑1 (MQ=255)
cGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGga < 1:3491921/36‑1 (MQ=255)
gAAAGTCCGGCGACAGCTGCGACCAGTGGCGGggaa > 1:4986795/1‑36 (MQ=255)
aaaGTCCGGCGACAGCTGCGACCAGTGGCGGgaaag < 1:6022126/36‑1 (MQ=255)
aaGTCCGGCGACAGCTGCGACCAGTGGCGGgtaaga > 1:513164/1‑36 (MQ=255)
aGTCCGGCGACAGCTGCGACCAGTGGCGGgaaagaa < 1:4164052/36‑1 (MQ=255)
aGTCCGGCGACAGCTGCGACCAGTGGCGGgaaagaa > 1:6290984/1‑36 (MQ=255)
gTCCGGCGACAGCTGCGACCAGTGGCGGggaagaaa > 1:7366278/1‑36 (MQ=255)
gTCCGGCGACAGCTGCGACCAGTGGCGGggaagaaa > 1:7218655/1‑36 (MQ=255)
gTCCGGCGACAGCTGCGACCAGTGGCGGgaaagaaa > 1:1390412/1‑36 (MQ=255)
gTCCGGCGACAGCTGCGACCAGTGGCGGgaaagaaa < 1:5332403/36‑1 (MQ=255)
tCCGGCGACAGCTGCGACCAGTGGCGGgaaagaaaa > 1:5511123/1‑35 (MQ=255)
tCCGGCGACAGCTGCGACCAGTGGCGGGaaagaaag < 1:4604926/36‑1 (MQ=255)
tCCGGCGACAGCTGCGACCAGTGGCGGGaaagaaag < 1:5108603/36‑1 (MQ=255)
ccGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGt < 1:1041500/36‑1 (MQ=255)
cGGCGACAGCTGCGACCAGTGGCGGGgaagaaagag > 1:3376985/1‑34 (MQ=37)
cGGCGACAGCTGCGACCAGTGGCGGGGAAGAAAATg > 1:2831257/1‑36 (MQ=37)
ggCGACAGCTGCGACCAGTGGCGGGGAAGAAAGGGa > 1:3343233/1‑36 (MQ=37)
ggCGACAGCTGCGACCAGTGGCGGGGAAGAAAGAGa > 1:5311355/1‑36 (MQ=37)
gCGACAGCTGCGACCAGTGGCGGGGAAGAAAGGGAt > 1:3260400/1‑36 (MQ=37)
gCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGAt < 1:2680947/36‑1 (MQ=255)
gCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGAt < 1:7432576/36‑1 (MQ=255)
cGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATc > 1:6829108/1‑36 (MQ=255)
cGACAGCTGCGACCAGTGGCGGGAAAGAAAATGGTc > 1:5098546/1‑36 (MQ=25)
gACAGCTGCGACCAGTGGCGGGGAAGAAAGGGATcc > 1:4798762/1‑36 (MQ=37)
aCAGCTGCGACCAGTGGCGGGAAAGAAAGTGATccc < 1:3859293/36‑1 (MQ=255)
aCAGCTGCGACCAGTGGCGGGAAAGAAAGTGATccc > 1:296220/1‑36 (MQ=255)
cAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCca < 1:933156/36‑1 (MQ=255)
cAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCca < 1:566375/36‑1 (MQ=255)
aGCTGCGACCAGTGGGGGGAAAGAAAGTGATCCcac < 1:3029580/36‑1 (MQ=255)
aGCTCCGACCAGTGGCGGGAAAGAAAGTGATCCcac < 1:1066417/36‑1 (MQ=255)
gCTGCGACCAGTGGCGGGAAAGAAAGTGATCCcaca < 1:1557267/36‑1 (MQ=255)
tGCGACCAGTGGCGGGGAAGAAAATGATCCCCCATg > 1:2106821/1‑36 (MQ=25)
cGACCAGTGGCGGGAAAGAAAGTGATCCCACATGAc < 1:4231336/36‑1 (MQ=255)
cGACCAGTGGCGGGAAAGAAAGTGATCCCACATGAc < 1:7758527/36‑1 (MQ=255)
cGACCAGTGGCGGGAAAGAAAGTGATCCCACATGAc < 1:7050790/36‑1 (MQ=255)
gACCAGTGGCGGGAAAGAAAGTGATCCCACATGACa > 1:6730802/1‑36 (MQ=255)
aGTGGCGGGAAAGAAAGTGATCCCACATGACATCCa < 1:339929/36‑1 (MQ=255)
aGTGGCGGGAAAGAAAGTGATCCCACATGACATCCa < 1:1318505/36‑1 (MQ=255)
tGGCGGGAAAGAAAGTGATCCCACATGACATCCAGc > 1:7772886/1‑36 (MQ=255)
tGGCGGGAAAGAAAGTGATCCCACATGACATCCAGc > 1:6432073/1‑36 (MQ=255)
ggCGGGAAAGAAAGTGATCCCACATGACATCCAGCg > 1:3946057/1‑36 (MQ=255)
ggCGGGAAAGAAAGTGATCCCACATGACATCCAGCg < 1:1727205/36‑1 (MQ=255)
gCGGGAAAGAAAGTGATCCCACATGACATCCAGCGt < 1:358905/36‑1 (MQ=255)
cGGGAAAGAAAGTGATCCCACATGACATCCAGCGTa > 1:5834509/1‑36 (MQ=255)
cGGGAAAGAAAGTGATCCCACATGACATCCAGCGTa > 1:643453/1‑36 (MQ=255)
gggAAAGAAAGTGATCCCACATGACATCCAGCGTaa > 1:6710595/1‑36 (MQ=255)
gggAAAGAAAGTGATCCCACATGACATCCAGCGTaa > 1:6461996/1‑36 (MQ=255)
gggAAAGAAAGTGATCCCACATGACATCCAGCGTaa < 1:5002028/36‑1 (MQ=255)
gggAAAGAAAGTGATCCCACATGACATCCAGCGTaa > 1:5977256/1‑36 (MQ=255)
ggAAAGAAAGTGATCCCACATGACATCCAGCGTAAt > 1:1625624/1‑36 (MQ=255)
ggAAAGAAAGTGATCCCACATGACATCCAGCGTAAt < 1:2778696/36‑1 (MQ=255)
aaagaaagTGATCCCACATGACATCCAGCGTAATAg > 1:4205106/1‑36 (MQ=255)
aaagaaagTGATCCCACATGACATCCAGCGTAATAg < 1:2980368/36‑1 (MQ=255)
|
CGGAAAGTCCGGCGACAGCTGCGACCAGTGGCGGGAAAGAAAGTGATCCCACATGACATCCAGCGTAATAG > REL606/393051‑393121
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 10 ≤
ATCG
< 20 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|