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

AATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTT  >  NC_000913/3979726‑3979814                                            |                                              aaTCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCt                                          >  1:1444423/1‑49 (MQ=255) aaTCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCt                                          >  1:207892/1‑49 (MQ=255)  aTCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCtt                                         >  1:770400/1‑49 (MQ=255)  aTCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCtt                                         <  1:1677148/49‑1 (MQ=255)  aTCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCt                                          <  1:34902/48‑1 (MQ=255)  aTCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCt                                          <  1:2785524/48‑1 (MQ=255)   tCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCtt                                         >  1:207622/1‑48 (MQ=255)    cGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTc                                        >  1:1239982/1‑48 (MQ=255)     gCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCaa                                      >  1:1827314/1‑49 (MQ=255)     gCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCaa                                      >  1:836427/1‑49 (MQ=255)     gCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCaa                                      >  1:290929/1‑49 (MQ=255)     gCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCa                                       <  1:1708631/48‑1 (MQ=255)      cTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCaa                                      >  1:735290/1‑48 (MQ=255)        ataAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAatat                                   >  1:2354105/1‑49 (MQ=255)        ataAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAatat                                   >  1:1801407/1‑49 (MQ=255)         taAAGCTGCGATATTGCACATTTTCTGCTTTGGGGTGATCTTCAata                                    <  1:2554929/47‑1 (MQ=255)         taAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAatat                                   <  1:1408540/48‑1 (MQ=255)         taAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAatat                                   <  1:2091543/48‑1 (MQ=255)         taAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAatat                                   >  1:2229397/1‑48 (MQ=255)         taAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAata                                    <  1:939784/47‑1 (MQ=255)         taAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAata                                    <  1:2024405/47‑1 (MQ=255)         taAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATg                                  <  1:317512/49‑1 (MQ=255)           aaGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGa                                 <  1:2160804/48‑1 (MQ=255)           aaGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGAt                                >  1:1601212/1‑49 (MQ=255)           aaGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGAt                                >  1:1128122/1‑49 (MQ=255)           aaGCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGAt                                <  1:584810/49‑1 (MQ=255)             gCTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATct                              <  1:412914/49‑2 (MQ=255)              cTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGt                             >  1:685979/1‑49 (MQ=255)              cTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGt                             <  1:711938/49‑1 (MQ=255)              cTGCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGt                             >  1:928106/1‑49 (MQ=255)                gCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGc                           >  1:1653945/1‑49 (MQ=255)                gCGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGc                           >  1:2841566/1‑49 (MQ=255)                 cGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTg                            <  1:664297/47‑1 (MQ=255)                 cGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCt                          <  1:51264/49‑1 (MQ=255)                 cGATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCt                          <  1:244688/49‑1 (MQ=255)                  gATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTg                         >  1:2479057/1‑49 (MQ=255)                  gATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTg                         >  1:770043/1‑49 (MQ=255)                  gATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTg                         <  1:2802761/49‑1 (MQ=255)                  gATATTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTg                         <  1:1996665/49‑1 (MQ=255)                   atatTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTgg                        >  1:2319895/1‑49 (MQ=255)                    tatTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGgc                       >  1:1675212/1‑49 (MQ=255)                    tatTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGgc                       <  1:379996/49‑1 (MQ=255)                    tatTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGgc                       <  1:652985/49‑1 (MQ=255)                     atTGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGgcg                      >  1:231863/1‑49 (MQ=255)                       tGCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGcgcg                    >  1:1508941/1‑49 (MQ=255)                        gCACAGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGcgcg                    <  1:1191280/48‑1 (MQ=255)                         cacaGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTg                  <  1:837257/49‑1 (MQ=255)                          acaGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGa                 <  1:39754/49‑1 (MQ=255)                            aGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAg               >  1:2798019/1‑49 (MQ=255)                            aGTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAg               >  1:2503930/1‑49 (MQ=255)                             gTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAgg              >  1:507954/1‑49 (MQ=255)                             gTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAgg              <  1:736391/49‑1 (MQ=255)                             gTTTCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAgg              <  1:1006561/49‑1 (MQ=255)                               ttCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGc            >  1:935504/1‑49 (MQ=255)                                tCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGc            <  1:1183325/48‑1 (MQ=255)                                tCTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCt           >  1:1745946/1‑49 (MQ=255)                                 cTGCTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTg          >  1:2415710/1‑49 (MQ=255)                                    cTTTGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGAt       >  1:648847/1‑49 (MQ=255)                                     tttGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGAtt      >  1:1383388/1‑49 (MQ=255)                                       tGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTc     <  1:2274100/48‑1 (MQ=255)                                       tGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTc     <  1:2057857/48‑1 (MQ=255)                                       tGGGGTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCg    >  1:340785/1‑49 (MQ=255)                                        ggggTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGt   >  1:694116/1‑49 (MQ=255)                                         gggTGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGtt  >  1:1113782/1‑49 (MQ=255)                                            |                                              AATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTT  >  NC_000913/3979726‑3979814

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