Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,066,123 T→C Q161R (CAG→CGG)  cobT ← nicotinate‑nucleotide dimethylbenzimidazole‑P phophoribosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,066,1230TC90.5% 234.1 / 11.0 85Q161R (CAG→CGG) cobTnicotinate‑nucleotide dimethylbenzimidazole‑P phophoribosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (0/8);  new base C (76/0);  total (77/8)
Fisher's exact test for biased strand distribution p-value = 2.29e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TTTGCCATCCCCAGTTCACCTACACCAAACAGCGTGACACCGTTTTTTGCCAGCTCCTGCGTAT  >  W3110S.gb/2066066‑2066129
                                                         |      
tttGCCATCCCCAGTTCACCTACACCAAACAGCGTGACACCGTTTTTTGCCAGCTCCTGCGt    <  1:1493386/62‑1 (MQ=255)
tttGCCATCCCCAGTTCACCTACACCAAACAGCGTGACACCGTTTTTTGCCAGCTCCTGCGt    <  1:356063/62‑1 (MQ=255)
tttGCCATCCCCAGTTCACCTACACCAAACAGCGTGACACCGTTTTTTGCCAGCTCCTGCGt    <  1:2043560/62‑1 (MQ=255)
tttGCCATCCCCAGTTCACCTACACCAAACAGCGTGACACCGTTTTTTGCCAGCTCCTGCGt    <  1:2700160/62‑1 (MQ=255)
tttGCCATCCCCAGTTCACCTACACCAAACAGCGTGACACCGTTTTTTGCCAGCTCCTGCGt    <  1:1424854/62‑1 (MQ=255)
tttGCCATCCCCAGTTCACCTACACCAAACAGCGGGACACCGTTTTTTGCCAGCTCCTGCGt    <  1:2735413/62‑1 (MQ=255)
 ttGCCATCCCCAGTTCACCTACACCAAACAGCGTGACACCGTTTTTTGCCAGCTCCTGCGt    <  1:135978/61‑1 (MQ=255)
 ttGCCATCCCCAGTTCACCTACACCAAACAGCGTGACACCGTTTTTTGCCAGCTCCTGCGt    <  1:1481553/61‑1 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2491921/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2490938/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:333483/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:317495/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2876933/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2851619/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2851025/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2748318/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2736240/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2670099/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2616627/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2717693/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2561215/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:490880/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2675470/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2623759/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2575207/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:634888/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:961039/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:881997/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:796818/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:758441/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:735047/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:727396/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:72227/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:719411/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:637674/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:37825/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:628966/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:625307/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:604252/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:596199/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:584555/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:574116/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:521995/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2310922/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:463631/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1602985/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1755908/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2388318/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1749053/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1744406/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1732899/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1724385/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1650646/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:161596/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1603605/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1752952/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1598590/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:150839/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:141634/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1408176/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1215083/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1094202/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1091027/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1087969/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2084980/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2340050/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2317011/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1026078/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2302448/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2295774/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2266650/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2214773/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1818521/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2074901/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2048949/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2025625/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1965940/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1962917/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1877108/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1870377/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:1845359/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGCCAGCTCCAGCGtat  >  1:1741235/1‑36 (MQ=255)
                            aCAGCGTGACACCGTTTTTTGACAGCTCCCGCGtat  >  1:120643/1‑36 (MQ=37)
                            aCAGCGTGACACCGTTTTATGCCAGCTCCCGCGtat  >  1:1486192/1‑36 (MQ=37)
                            aCAGCGTCACACCGTTTTTTGCCAGCTCCCGCGtat  >  1:2307615/1‑36 (MQ=37)
                                                         |      
TTTGCCATCCCCAGTTCACCTACACCAAACAGCGTGACACCGTTTTTTGCCAGCTCCTGCGTAT  >  W3110S.gb/2066066‑2066129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: