Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,768,203 | 0 | G→T | 9.5% | 2.5 ‑10.4 0.0 | 85 | intergenic (‑5/+45) | ygbE/cysC | conserved inner membrane protein/adenylylsulfate kinase |
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands. | ||||||||||
Reads supporting (aligned to +/- strand): new base (10/0): ref base (23/52): total (33/52) | ||||||||||
Fisher's exact test strand distribution p-value = 2.96e-05 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
TTGTTAGTGTAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCC > REL606/2768168‑2768238
|
ttGTTAGTGTAATGTTATGGCTATTACGCATATTTg < 1:42803/36‑1 (MQ=255)
ttGTTAGTGTAATGTTATGGCTATTACGCATATTTg < 1:3601768/36‑1 (MQ=255)
gTTAGTGTAATGTTATGGCTATTACGCCAATTTTAt > 1:2556261/1‑36 (MQ=25)
gTTAGTGTAATGTTATGGCTATTACGCATATTTGAt < 1:3583497/36‑1 (MQ=255)
gTTAGTGTAATGTTATGGCTATTACGCATATTTGAt < 1:4630183/36‑1 (MQ=255)
gTTAGTGTAATGTTATGGCTATTACGCATATTTGAt < 1:5118570/36‑1 (MQ=255)
ttAGTGTAATGTTATGGCTATTACGCATATTTTATc > 1:3651341/1‑36 (MQ=255)
ttAGTGTAATGTTATGGCTAATACCCATATTTTATc > 1:875373/1‑36 (MQ=25)
tAGTGTAATGTTATGGCTATTTCGCATATTTGATcc < 1:3617659/36‑1 (MQ=255)
tAGTGTAATGTTATGGCTATTACGCATATTTTAttc > 1:804270/1‑34 (MQ=37)
tAGTGTAATGTTATGGCTATTACGCATATTTGATcc < 1:385517/36‑1 (MQ=255)
tAGTGTAATGTTATGGCTATTACGCATATTTGATcc < 1:5801844/36‑1 (MQ=255)
tAGTGTAATGTTATGGCTATTACGCATATTTGATcc < 1:2956918/36‑1 (MQ=255)
tAGTGTAATGTTATGGCTATTAAGCATATTTTATca > 1:5356116/1‑35 (MQ=37)
aGTGTAATGTTATGGCTATTACGCATATTTGATCCt < 1:5437385/36‑1 (MQ=255)
aGTGTAATGTTATGGCTATTACGCATATTTGATCCt < 1:1018564/36‑1 (MQ=255)
aGTGTAATGTTATGGCTATTAAGCATATTTTATCCt > 1:530074/1‑36 (MQ=37)
aGTGTAATGTGATGGCTATTACGCATATTTGATCCt < 1:5507213/36‑1 (MQ=255)
gtgtAATGTTATGGCTATTAAGCATATTTTATTCtt > 1:2570012/1‑35 (MQ=21)
tgtAATGTTATGGCTATTACGCATATTTGATCCTGt < 1:5487230/36‑1 (MQ=255)
tgtAATGTTATGGCTATTACGCATATTTGATCCTGt < 1:4545347/36‑1 (MQ=255)
tgtAATGTTATGGCTATTACGCATATTTGATCCTGt < 1:1085631/36‑1 (MQ=255)
tgtAATGTTATGGCTATTACGCATATTTGATCCTGt < 1:4758978/36‑1 (MQ=255)
tgtAATGTTATGGCTAATACTCATATTTTATTCttt > 1:2863667/1‑34 (MQ=16)
tgtAATGTTATGGCTAATAAGCATATTTGATCCttt > 1:5465812/1‑34 (MQ=37)
gtAATGTTATGGCTATTACGCATATTTGATCCTGtt < 1:1694459/36‑1 (MQ=255)
gtAATGTTATGGCTATTACGCATATTTGATCCTGtt < 1:4704975/36‑1 (MQ=255)
gtAATGTTATGGCTAATTCGCATATATGATCCTTtt > 1:1434244/1‑36 (MQ=19)
tAATGTTATGGCTATTACGCATATTTGATCCTGTTa < 1:2474385/36‑1 (MQ=255)
tAATGTTATGGCTATTACGCATATTTGATCCTGTTa < 1:2087476/36‑1 (MQ=255)
aaTGTTATGGCTATTACGCATATTTTATCCTTtttc > 1:224363/1‑34 (MQ=25)
aaTGTTATGGCTATTACGCATATTTGATCCTGTTAc < 1:60961/36‑1 (MQ=255)
aTGTTATGGCTATTACGCATATTTGATCCTGTTACt < 1:4298412/36‑1 (MQ=255)
aTGTTATGGCTATTACGCATATTTGATCCTGTTACt < 1:1890381/36‑1 (MQ=255)
aTGTTATGGCTATTACGCATATTTGATCCTGTTACt < 1:2262206/36‑1 (MQ=255)
tGTTATGGCTATTACGCATATTTGCTCCTGTTTCtt < 1:5767720/36‑1 (MQ=37)
gTTATGGCTATTACGCATATTTGATCCTGTTACttt < 1:383783/36‑1 (MQ=255)
gTTATGGCTATTACGCATATTTGATCCTGTTACttt < 1:4472580/36‑1 (MQ=255)
gTTATGGCTATTACGCATATTTGATCCTGTTACttt < 1:4711959/36‑1 (MQ=255)
gTTATGGCTATTACGCATATTTGATCCTGTTACttt < 1:2217971/36‑1 (MQ=255)
ttATGGCTATTACGCATATTTGATCCTGTTACTTTg < 1:645285/36‑1 (MQ=255)
tATGGCTATTACGCATATTTTATCCTGTTACtttta > 1:753336/1‑34 (MQ=38)
tATGGCTATTACGCATATTTGATCCTGTTACTTTGa < 1:421654/36‑1 (MQ=255)
aTGGCTATTACGCATATTTGATCCTGTTACTTTGAc < 1:3249875/36‑1 (MQ=255)
aTGGCTATTACGCATATTTGATCCTGTTACTTTGAc < 1:5643338/36‑1 (MQ=255)
aTGGCTATTACGCATATTTGATCCTGTTACTTTGAc > 1:4164678/1‑36 (MQ=255)
aTGGCTAATGCGCATATTTGATCCTGTTACTTTGAc < 1:546172/36‑1 (MQ=37)
tGGCTATTACGCATATTTGATCCTGTTACTTTGACg > 1:4019160/1‑36 (MQ=255)
gCTATTACGCATATTTGATCCTGTTACTTTGACGaa > 1:3977870/1‑36 (MQ=255)
gCTATTACGCATATTTGATCCTGTTACTTTGACGaa < 1:1333918/36‑1 (MQ=255)
gCTATTACGCATATTTGATCCTGTTACTTTGACGaa < 1:3378148/36‑1 (MQ=255)
gCTATTACGCATATTTGATCCTGTTACTTTGACGaa > 1:2759003/1‑36 (MQ=255)
gCTATTACGCATATTTGATCCTGTTAATTTGACGaa > 1:4023568/1‑36 (MQ=255)
cTATTACGCATATTTGATCCTGTTACTTTGACGAAc > 1:824952/1‑36 (MQ=255)
cTATTACGCATATTTGATCCTGTTACTTTGACGAAc < 1:1218515/36‑1 (MQ=255)
cTATTACGCATATTTGATCCTGTTACTTTGACGAAc < 1:214090/36‑1 (MQ=255)
tATTACGCATATTTGATCCTGTTACTTTGACGAAcc < 1:964447/36‑1 (MQ=255)
tATTACGCATATTTGATCCTGTTACTTTGACGAAcc < 1:1032148/36‑1 (MQ=255)
tATTACGCATATTTGATCCTGTTACTTTGACGAAcc < 1:66345/36‑1 (MQ=255)
aTTACGCATATTTGCTCCTGTTACTTTGACGAACCg < 1:5256241/36‑1 (MQ=255)
aTTACGCATATTTGATCCTGTTACTTTGACGAACCg > 1:4177002/1‑36 (MQ=255)
aTTACGCATATTTGATCCTGTTACTTTGACGAACCg > 1:1885279/1‑36 (MQ=255)
ttACGCATATTTGATCCTGTTACTTTGACGAACCgg > 1:2803240/1‑36 (MQ=255)
ttACGCATATTTGATCCTGTTACTTTGACGAACCgg > 1:2269533/1‑36 (MQ=255)
tACGCATATTTGATCCTGTTACTTTGACGAACCggg > 1:870501/1‑36 (MQ=255)
tACGCATATTTGATCCTGTTACTTTGACGAAACggg > 1:4853575/1‑36 (MQ=255)
aCGCATATTTGATCCTGTTACTTTGACGAACCGGGc > 1:2509976/1‑36 (MQ=255)
gCATATTTGATCCTGTTACTTTTACGAACCGGGCAt < 1:2951062/36‑1 (MQ=255)
gCATATTTGATCCTGTTACTTTGACGAACCGGGCAt < 1:3619391/36‑1 (MQ=255)
gCATATTTGATCCTGTTACTTTGACGAACCGGGCAt > 1:5559265/1‑36 (MQ=255)
gCATATTTGATCCTGTTACTTTGACGAACCGGGCAt < 1:4602067/36‑1 (MQ=255)
gCATATTTGATCCTGTTACTGTGACGAAACGGGGAt > 1:2377032/1‑36 (MQ=37)
cATATTTGATCCTGTTACTTTGACGAACCGGGCATg < 1:885511/36‑1 (MQ=255)
cATATTTGATCCTGTTACTTTGACGAACCGGGCATg < 1:3002551/36‑1 (MQ=255)
cATATTTGATCCTGTTACTTTGACGAACCGGGCATg > 1:5689302/1‑36 (MQ=255)
cATATTTGATCCTGTTACTTTGACGAACCGGGCATg < 1:4332479/36‑1 (MQ=255)
atatTTGATCCTGTTCCTTTGACGAACCGGGCATgg < 1:4571546/36‑1 (MQ=255)
atatTTGATCCTGTTACTTTGACGAACCGGGCATgg > 1:3689177/1‑36 (MQ=255)
atatTTGATCCTGTTACTTTGACGAACCGGGCATgg > 1:284445/1‑36 (MQ=255)
tatTTGATCCTGTTACTTTGACTAACCGGGCATGGa < 1:1968238/36‑1 (MQ=255)
atTTGATCCTGTTACTTTGACGAACCGGGCATGGaa < 1:5818526/36‑1 (MQ=255)
atTTGATCCTGTTACTTTGACGAACCGGGCATGGaa > 1:2340161/1‑36 (MQ=255)
atTTGATCCTGTTACTTTGACGAACCGGGCATGGaa > 1:275846/1‑36 (MQ=255)
tttGATCCTGTTACTTTGACGAACCGGGCATGGaaa < 1:1471334/36‑1 (MQ=255)
tttGATCCTGTTACTTTGACGAACCGGGCATGGaaa > 1:5191111/1‑36 (MQ=255)
ttGATTCTGTTACTTTGACTACCCGGGCATGGAAAc < 1:1498496/36‑1 (MQ=25)
ttGATCCTGTTACTTTGACGAACCGGGCATGGAAAc < 1:3420353/36‑1 (MQ=255)
gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc < 1:1564766/36‑1 (MQ=255)
gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc > 1:604849/1‑36 (MQ=255)
gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc > 1:4348973/1‑36 (MQ=255)
gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc < 1:2841623/36‑1 (MQ=255)
gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc > 1:5055288/1‑36 (MQ=255)
gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc < 1:3858514/36‑1 (MQ=255)
gATCCTGTTACTTTGACGAACCGGGCATGGAAAccc > 1:896148/1‑36 (MQ=255)
|
TTGTTAGTGTAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCC > REL606/2768168‑2768238
Base quality scores: 
ATCG
< 3 ≤
ATCG
< 10 ≤
ATCG
< 21 ≤
ATCG
< 30 ≤
ATCG
< 37 ≤
ATCG
|