Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,979,769	C→T	A390V (GCG→GTG)	wzyE →	enterobacterial common antigen polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,979,769	0	C	T	100.0%	129.1 / NA	33	A390V (GCG→GTG)	wzyE	enterobacterial common antigen polymerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/15);  total (18/15)

CCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTT  >  NC_000913/3979723‑3979815                                               |                                               ccTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGAt                                             >  1:512506/1‑50 (MQ=255)  cTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGTTTTGGGGTGATc                                            <  1:925700/50‑1 (MQ=255)  cTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATc                                            >  1:2385757/1‑50 (MQ=255)  cTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATc                                            >  1:1755571/1‑50 (MQ=255)     aTCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTc                                         <  1:7450/50‑1 (MQ=255)     aTCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTc                                         <  1:2086328/50‑1 (MQ=255)        gCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCaa                                       >  1:2296837/1‑49 (MQ=255)        gCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCaa                                       >  1:1733264/1‑49 (MQ=255)         cTATCAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAata                                     <  1:1974869/50‑1 (MQ=255)          tataAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAatat                                    >  1:1047807/1‑50 (MQ=255)            taAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAatat                                    <  1:1158979/48‑1 (MQ=255)            taAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGa                                  >  1:1060050/1‑50 (MQ=255)              aaGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGAt                                 <  1:2734368/49‑1 (MQ=255)              aaGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGAt                                 <  1:2209112/49‑1 (MQ=255)                gCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGt                              <  1:1237974/50‑1 (MQ=255)                 cTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTg                             >  1:1877361/1‑50 (MQ=255)                 cTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTg                             >  1:687984/1‑50 (MQ=255)                    cGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTg                          >  1:1925647/1‑50 (MQ=255)                       tatTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGgcg                       >  1:2423287/1‑50 (MQ=255)                         ttGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGgcgc                      <  1:1517387/49‑1 (MQ=255)                              caGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAg                >  1:16109/1‑50 (MQ=255)                               aGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAgg               >  1:899426/1‑50 (MQ=255)                                 tttCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGc             >  1:643179/1‑50 (MQ=255)                                   tCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTg           >  1:994661/1‑50 (MQ=255)                                   tCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTg           >  1:1759798/1‑50 (MQ=255)                                    cTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTg           <  1:2547802/49‑1 (MQ=255)                                    cTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTg           <  1:1058694/49‑1 (MQ=255)                                       cTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGAtt       >  1:256439/1‑50 (MQ=255)                                       cTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGAtt       >  1:93471/1‑50 (MQ=255)                                       cTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGAt        <  1:2331563/49‑1 (MQ=255)                                            gggTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGttt  <  1:1690523/50‑1 (MQ=255)                                            gggTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGttt  <  1:1150456/50‑1 (MQ=255)                                            gggTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGtt   <  1:864528/49‑1 (MQ=255)                                               |                                               CCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTT  >  NC_000913/3979723‑3979815

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