Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 3,304,275 | 0 | C→A | 13.2% | 4.2 ‑9.5 0.0 | 86 | 400 | sspB | ClpXP protease specificity‑enhancing factor |
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands. | ||||||||||
Reads supporting (aligned to +/- strand): new base (0/15): ref base (38/33): total (38/48) | ||||||||||
Fisher's exact test strand distribution p-value = 7.22e-05 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
ATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATG > REL606/3304240‑3304310
|
aTGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc < 1:7400594/36‑1 (MQ=255)
tGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa < 1:6660308/36‑1 (MQ=255)
tGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa < 1:4598952/36‑1 (MQ=255)
tGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa < 1:4330236/36‑1 (MQ=255)
tGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa < 1:4321748/36‑1 (MQ=255)
tGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa < 1:1808831/36‑1 (MQ=255)
gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa > 1:7561913/1‑36 (MQ=255)
gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa < 1:2591158/36‑1 (MQ=255)
gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa < 1:3564865/36‑1 (MQ=255)
gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa < 1:209344/36‑1 (MQ=255)
gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa < 1:2353235/36‑1 (MQ=255)
aGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt < 1:4216631/36‑1 (MQ=255)
aGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt < 1:3253037/36‑1 (MQ=255)
gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa < 1:3587971/36‑1 (MQ=255)
gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa > 1:3808376/1‑36 (MQ=255)
gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa < 1:1943920/36‑1 (MQ=255)
tgtCATCATCGTGATCTGGCTTGTCGCCATCAATaa > 1:2967477/1‑36 (MQ=255)
tgtCATCATCGTGATCTGGCTTGTCGCCATCAATaa < 1:862914/36‑1 (MQ=255)
gtCATCATCGTGATCTGGCTTGTCGCCATCAATAAc > 1:5637189/1‑36 (MQ=255)
gtCATCATCGTGATCTGGCTTGTCGCCATCAATAAc < 1:4066736/36‑1 (MQ=255)
tcatcaTCGTGATCTGGCTTGTCGCCATCCATAAcc < 1:6643762/36‑1 (MQ=255)
tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc < 1:6809978/36‑1 (MQ=255)
tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc < 1:3415714/36‑1 (MQ=255)
tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc < 1:3463209/36‑1 (MQ=255)
tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc < 1:2484326/36‑1 (MQ=255)
tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc < 1:3651558/36‑1 (MQ=255)
tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc < 1:4890253/36‑1 (MQ=255)
catcatCGTGATCTGGCTTGTCGCCATCAATAACCg > 1:1732722/1‑36 (MQ=255)
catcatCGTGATCTGGCTTGTCGCCATCAATAACCg > 1:7382900/1‑36 (MQ=255)
tcatcGTGATCTGGCTTGTCGCCATCAATAACCGAc < 1:761193/36‑1 (MQ=255)
tcatcGTGATCTGGCTTGTCGCCATCAATAACCGAc < 1:1773492/36‑1 (MQ=255)
tcatcGTGATCTGGCTTGTCGCCATCAATAACCGAc < 1:5054900/36‑1 (MQ=255)
catcGTGATCTGGCTTGTCGCCATCAATAACCGACa > 1:4344950/1‑36 (MQ=255)
atcGTGATCTGGCTTTTCGCCATCAATAACCGACAt < 1:6856034/36‑1 (MQ=255)
atcGTGATCTGGCTTGTCGCCATCAATAACCGACAt > 1:4866153/1‑36 (MQ=255)
atcGTGATCTGGCTTGTCGCCATCAATAACCGACAt > 1:4639302/1‑36 (MQ=255)
atcGTGATCTGGCTTGTCGCCATCAATAACCGACAt < 1:1897014/36‑1 (MQ=255)
atcGTGATCTGGCTTGTCGCCATCAATAACCGACAt > 1:5299033/1‑36 (MQ=255)
atcGTGATCTGGCTTGTCGCCATCAATAACCGACAt > 1:5301662/1‑36 (MQ=255)
agcGTGATCTGGCTTGTCGCCATCAATAACCGACAt < 1:2390317/34‑1 (MQ=255)
tcGTGATCTGGCTTGTCGCCATCAATAACCGACATa > 1:4422839/1‑36 (MQ=255)
cGTGATCTGGCTTGTCGCCATCAATAACCGACATaa > 1:7000089/1‑36 (MQ=255)
gTGATCTGGCTTGTCGCCATCAATAACCGACATAAc > 1:3827992/1‑36 (MQ=255)
gTGATCTGGCTTGTCGCCAAAAAAAACCGACATAAc < 1:7124110/36‑1 (MQ=21)
tGATCTGGTTTGTCACCATAAATAACCGACATAACg < 1:4773768/36‑1 (MQ=25)
tGATCTGGCTTGTCGCCATCAATAACCGACATAACg > 1:3678126/1‑36 (MQ=255)
gATCTGTCTTCTCGCAAAAAATAACCGACATAACgg < 1:1336613/36‑1 (MQ=12)
gATCTGGCTTGTCGCCATCAATAACCGACATAACgg > 1:842225/1‑36 (MQ=255)
gATCTGGCTTGTCGCCATCAATAACCGACATAACgg > 1:7506715/1‑36 (MQ=255)
gATCTGGCTTGTCGCCATAAATAACCGACATAACgg < 1:59899/36‑1 (MQ=255)
aTCTGGCTTGTCGCCATAAAAAACCGACATAACGGt < 1:1822065/36‑1 (MQ=37)
tCTGGCTTGTCGCCATCAATAACCGACATAACGGtt > 1:1957554/1‑36 (MQ=255)
cgggCTTGTCGCAATAAATAACCGACATAACGGttt < 1:5702652/34‑1 (MQ=37)
cTGGCTTGTCGCCATCAATAACCGACATAACGGttt > 1:485364/1‑36 (MQ=255)
cTGGCTTGTCGCCATCAATAACCGACATAACGGttt > 1:2981015/1‑36 (MQ=255)
cTGGCTTGTCGCCATCAATAACCGACATAACGGttt > 1:6005296/1‑36 (MQ=255)
cTGGCTTGTCGCCATCAATAACCGACATAACGGttt > 1:7551454/1‑36 (MQ=255)
cTGGCTTGTCGCCATAAATAACCGACATAACGGttt < 1:6117806/36‑1 (MQ=255)
cTGGCTTGTCGCAATAAATAACCGACATAACGGttt < 1:587298/36‑1 (MQ=37)
tGGCTTGTCGCCATCAATAACCGACATAACGGTTTc < 1:7587774/36‑1 (MQ=255)
tGGCTTGTCGCCATCAATAACCGACATAACGGTTTc < 1:6628959/36‑1 (MQ=255)
tGGCTTGTCGCAATCAATAACCGACATAACGGTTTc < 1:4883891/36‑1 (MQ=255)
ggCTTGTCGCCATCAATAACCGACATAACGGTTTCg > 1:525499/1‑36 (MQ=255)
gCTTGTCGCCATCAAAAACCGACATAACGGTTTCGt < 1:2551649/36‑1 (MQ=255)
cTTGTCGCCATCAATAACCGACATAACGGTTTCGtt > 1:4439245/1‑36 (MQ=255)
cTTGTCGCCATAAATAACCGACATAACGGTTTCGtt < 1:3313446/36‑1 (MQ=255)
cTTGTCGCCATAAATAACCGACATAACGGTTTCGtt < 1:7812143/36‑1 (MQ=255)
cTTGTCACCATAAATAACCGACATAACGGTTTCGtt < 1:3714796/36‑1 (MQ=37)
cTTGCCGCCATAAATAACCGACATAACGGTTTCGtt < 1:2543404/36‑1 (MQ=37)
cTTGCCGCCATAAAAACCCGACATAACGGTTTCGtt < 1:6368225/36‑1 (MQ=19)
ttGTCGCCATCAATAACCGACATAACGGTTTCGTTg > 1:2651840/1‑36 (MQ=255)
ttGTCGCCATCAATAACCGACATAACGGTTTCGTTg > 1:6995327/1‑36 (MQ=255)
tGTCGCCATCAATAACCGACATAACGGTTTCGTTGt > 1:2896478/1‑36 (MQ=255)
gTCGCCATCAATAACCGACATAACGGTTTCGTTGTc > 1:6637869/1‑36 (MQ=255)
tCGCCATCAATAACCGACATAACGGTTTCGTTGTCt > 1:1549311/1‑36 (MQ=255)
tCGCCATCAATAACCGACATAACGGTTTCGTTGTCt > 1:425605/1‑36 (MQ=255)
tCGCCATAAATAACCGACATAACGGTTTCGTTGTCt < 1:5470358/36‑1 (MQ=255)
gCCATCAATAACCGACATAACGGTTTCGTTGTCTGc > 1:6963054/1‑36 (MQ=255)
gCCATCAATAACCGACATAACGGTTTCGTTGTCTGc > 1:6562222/1‑36 (MQ=255)
aaaaTCAATCACCGACATAACGGTTTCGTTGTCTGc < 1:403527/33‑1 (MQ=37)
cATCAATAACCGACATAACGGTTTCGTTGTCTGCCg > 1:6021335/1‑36 (MQ=255)
cATCAATAACCGACATAACGGTTTCGTTGTCTGCCg > 1:944769/1‑36 (MQ=255)
cATAAATAACCGACATAACGGTTTCGTTGTCTGCCg < 1:5430382/36‑1 (MQ=255)
aacAATAACCGACATAACGGTTTCGTTGTCTGCCGa < 1:1953447/34‑1 (MQ=255)
aTCAATAACCGACATAACGGTTTCGTTGTCTGCCGa > 1:739814/1‑36 (MQ=255)
aTCAATAACCGACATAACGGTTTCGTTGTCTGCCGa > 1:5710992/1‑36 (MQ=255)
aTCAATAACCGACATAACGGTTTCGTTGTCTGCCGa > 1:1746781/1‑36 (MQ=255)
tCAATAACCGACATAACGGTTTCGTTGTCTGCCGAt > 1:6816129/1‑36 (MQ=255)
cAATAACCGACATAACGGTTTCGTTGTCTGCCGATg > 1:5438916/1‑36 (MQ=255)
|
ATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATG > REL606/3304240‑3304310
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 10 ≤
ATCG
< 19 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|