Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,659,831 2 bp→AA coding (572‑573/1251 nt) wzxE ← O‑antigen translocase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,659,8310GA95.5% 133.8 / ‑1.7 44A191V (GCC→GTTwzxEO‑antigen translocase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (42/0);  total (43/1)
Fisher's exact test for biased strand distribution p-value = 4.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
*W3110S.gb3,659,8320GA84.0% 115.3 / 15.0 50A191V (GCC→GTT) wzxEO‑antigen translocase
Reads supporting (aligned to +/- strand):  ref base G (7/1);  new base A (42/0);  total (49/1)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAACC  >  W3110S.gb/3659773‑3659893
                                                          ||                                                             
tATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATggc                                                              >  1:435429/1‑61 (MQ=255)
tATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATggc                                                              >  1:2643187/1‑61 (MQ=255)
                             gggATGACACCACGTTTGATCAACATGATGGCGGCAGGa                                                       >  1:1749227/1‑39 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAGTTACCACCAGCGCGGGAATCAgcgc                             >  1:2182333/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2917134/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2026686/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2117113/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2280043/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2348655/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2367994/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2436253/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2506348/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2589145/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2593645/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1011386/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:354342/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:362496/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:513960/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:733462/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:833855/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:861890/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:901173/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:912172/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:948078/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:981523/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1974630/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1422134/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1335253/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1493599/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1318952/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1267539/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1532401/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1589848/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:166235/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1674158/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1240901/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1697437/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1767432/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1777830/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:178223/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1178782/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1826684/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:183791/1‑60 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcg                              >  1:449877/1‑59 (MQ=255)
                                  gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGAGCGGGAATCAgcg                              >  1:1864688/1‑59 (MQ=255)
                                         cGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCgg                                       <  1:2337825/43‑1 (MQ=255)
                                                         tgGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGAccc                      >  1:74753/1‑44 (MQ=255)
                                                         tgGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGccc            >  1:799071/1‑54 (MQ=255)
                                                         tgGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGccc            >  1:2960512/1‑54 (MQ=255)
                                                         tgGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGccc            >  1:2944920/1‑54 (MQ=255)
                                                         tgGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGccc            >  1:2580817/1‑54 (MQ=255)
                                                         tgGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGccc            >  1:1207528/1‑54 (MQ=255)
                                                           gcggcAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAAcc  >  1:1355745/1‑62 (MQ=255)
                                                           gcggcAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAAcc  >  1:70303/1‑62 (MQ=255)
                                                           gcggcAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAAcc  >  1:864696/1‑62 (MQ=255)
                                                           gcggcAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAAcc  >  1:1006622/1‑62 (MQ=255)
                                                           gcggcAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCATCATAAcc  >  1:1447934/1‑62 (MQ=255)
                                                           gcggcAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCAAGCAACGCCCCt          >  1:550307/1‑54 (MQ=255)
                                                          ||                                                             
TATCCCAGCTGGGTTTCAGATAGCTTAACGGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAACC  >  W3110S.gb/3659773‑3659893

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