Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,659,834	G→A	A190A (GCC→GCT)	wzxE ←	O‑antigen translocase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,659,834	0	G	A	84.0%	117.9 / 14.1	50	A190A (GCC→GCT)	wzxE	O‑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 = 1.00e+00

GGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAACC  >  W3110S.gb/3659802‑3659893                                 |                                                            gggATGACACCACGTTTGATCAACATGATGGCGGCAGGa                                                       >  1:1749227/1‑39 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAGTTACCACCAGCGCGGGAATCAgcgc                             >  1:2182333/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2917134/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1974630/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2026686/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:2117113/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:2280043/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:183791/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1011386/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1178782/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1240901/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1267539/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1318952/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1335253/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1422134/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1493599/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:1697437/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1826684/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:178223/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1777830/1‑60 (MQ=255)      gACACCACGTTTGATCAACATGATAACAGCAGGAATTACCACCAGCGCGGGAATCAgcgc                             >  1:1767432/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:1207528/1‑54 (MQ=255)                             tgGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGccc            >  1:2944920/1‑54 (MQ=255)                             tgGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGccc            >  1:2580817/1‑54 (MQ=255)                               gcggcAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAAcc  >  1:1355745/1‑62 (MQ=255)                               gcggcAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAAcc  >  1:70303/1‑62 (MQ=255)                               gcggcAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAAcc  >  1:1006622/1‑62 (MQ=255)                               gcggcAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAAcc  >  1:864696/1‑62 (MQ=255)                               gcggcAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCATCATAAcc  >  1:1447934/1‑62 (MQ=255)                               gcggcAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCAAGCAACGCCCCt          >  1:550307/1‑54 (MQ=255)                                 |                                                            GGGATGACACCACGTTTGATCAACATGATGGCGGCAGGAATTACCACCAGCGCGGGAATCAGCGCCAGACCCAGCAACGCCCCTTCATAACC  >  W3110S.gb/3659802‑3659893

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