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%	139.5 / NA	39	A390V (GCG→GTG)	wzyE	enterobacterial common antigen polymerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/17);  total (22/17)

GCCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTT  >  NC_000913/3979722‑3979814                                                |                                              gCCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGa                                             <  1:364547/50‑1 (MQ=255)  ccTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGAt                                            <  1:1242628/50‑1 (MQ=255)     aaTCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCtt                                         >  1:2583666/1‑50 (MQ=255)       tCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCa                                       <  1:446144/50‑1 (MQ=255)          cTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAata                                    <  1:1270290/50‑1 (MQ=255)            ataAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAatat                                   >  1:1103407/1‑49 (MQ=255)            ataAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAatat                                   >  1:1617326/1‑49 (MQ=255)             taAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAatat                                   <  1:1295196/48‑1 (MQ=255)             taAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAatat                                   <  1:282886/48‑1 (MQ=255)             taAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAatat                                   <  1:1814841/48‑1 (MQ=255)             taAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATg                                  <  1:1851410/49‑1 (MQ=255)               aaGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGAt                                <  1:2179612/49‑1 (MQ=255)                aGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCg                              >  1:958551/1‑50 (MQ=255)                aGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCg                              >  1:2418747/1‑50 (MQ=255)                  cTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTg                            >  1:266737/1‑50 (MQ=255)                  cTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTg                            <  1:728379/50‑1 (MQ=255)                   tGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGc                           >  1:956862/1‑50 (MQ=255)                    gCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCt                          >  1:1362390/1‑50 (MQ=255)                      gATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTgg                        >  1:942050/1‑50 (MQ=255)                       atatTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGgc                       >  1:516928/1‑50 (MQ=255)                        tatTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGgcg                      >  1:874415/1‑50 (MQ=255)                        tatTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGgcg                      >  1:1349110/1‑50 (MQ=255)                            gCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCgtgt                   <  1:249008/49‑1 (MQ=255)                             cacaGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTg                  >  1:1022829/1‑49 (MQ=255)                             cacaGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTg                  <  1:5247/49‑1 (MQ=255)                              acaGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGaa                <  1:223844/50‑1 (MQ=255)                                 gTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAggg             >  1:1180112/1‑50 (MQ=255)                                    tCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCt           <  1:1339205/49‑1 (MQ=255)                                    tCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTg          >  1:1830603/1‑50 (MQ=255)                                     cTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTg          >  1:553520/1‑49 (MQ=255)                                       gCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGAt       >  1:532674/1‑50 (MQ=255)                                       gCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGAt       >  1:2109310/1‑50 (MQ=255)                                        cTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGAt       <  1:2045575/49‑1 (MQ=255)                                          ttGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCg    >  1:1861865/1‑50 (MQ=255)                                           tGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGt   >  1:494453/1‑50 (MQ=255)                                            ggggTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGtt  >  1:1652168/1‑50 (MQ=255)                                            ggggTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGtt  >  1:1528358/1‑50 (MQ=255)                                            ggggTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGt   <  1:2036420/49‑1 (MQ=255)                                             gggTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGtt  <  1:2412366/49‑1 (MQ=255)                                                |                                              GCCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTT  >  NC_000913/3979722‑3979814

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