Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,969,234	G→A	M68I (ATG→ATA)	wzzE →	enterobacterial common antigen polysaccharide co‑polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,969,234	0	G	A	100.0%	90.6 / NA	26	M68I (ATG→ATA)	wzzE	enterobacterial common antigen polysaccharide co‑polymerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/14);  total (12/14)

TCGCGCTGGCGTATACGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCAAACATGGCTTCTGCCGACCAACC  >  NC_000913/3969158‑3969311                                                                             |                                                                              tCGCGCTGGCGTATACGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTgg                                                                           >  1:778656/1‑81 (MQ=255)   gcgcTGGCGTATACGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTgggg                                                                         >  1:778644/1‑81 (MQ=255)      cTGGCGTATACGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGa                                                                      >  2:778659/1‑81 (MQ=255)      cTGGCGTATACGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGa                                                                      >  1:778649/1‑81 (MQ=255)         gCGTATACGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGata                                                                    <  2:778609/80‑1 (MQ=255)           gTATACGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTa                                                                 >  1:778645/1‑81 (MQ=255)            tataCGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTAc                                                                <  1:778610/81‑1 (MQ=255)            tataCGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTAc                                                                <  2:778611/81‑1 (MQ=255)              taCGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTc                                                              <  2:778612/81‑1 (MQ=255)               aCGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCg                                                             >  1:778654/1‑81 (MQ=255)                     tttGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCgcagcag                                                       >  1:778658/1‑81 (MQ=255)                      ttGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCgcagcag                                                       >  1:778655/1‑80 (MQ=255)                        gCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCGCAGCAGCaa                                                    >  1:778653/1‑81 (MQ=255)                         cTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCGCAGCAGCAAt                                                   <  2:778613/81‑1 (MQ=255)                           cGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCGCAGCAGCAAttt                                                 >  2:778647/1‑81 (MQ=255)                               aGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCGCAGCAGCAATTTTTg                                              <  1:778615/80‑1 (MQ=255)                               aGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCGCAGCAGCAATTTTTGc                                             <  1:778616/81‑1 (MQ=255)                                 gAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGt                                           <  2:778617/81‑1 (MQ=255)                                            gACGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGAt                                  <  1:778618/79‑1 (MQ=255)                                              cGGCGATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTcc                              <  1:778619/81‑1 (MQ=255)                                                  gATTACCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTc                          <  1:778620/81‑1 (MQ=255)                                                      aCCGATCGTCCAACGGTGAATATACTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCaaa                       <  2:778621/80‑1 (MQ=255)                                                                cAACGGTGAATATACTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCAAACATGGCTTCt             <  1:778622/80‑1 (MQ=255)                                                                   cGGTGAATATACTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCAAACATGGCTTCTGcc          >  2:778676/1‑80 (MQ=255)                                                                    ggTGAATATACTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCAAACATGGCTTCTGCCGa        >  1:778662/1‑81 (MQ=255)                                                                     gTGAATATACTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCAAACATGGCTTCTGCCGac       <  2:778623/81‑1 (MQ=255)                                                                          tataCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCAAACATGGCTTCTGCCGaccaacc  <  2:778624/81‑1 (MQ=255)                                                                             |                                                                              TCGCGCTGGCGTATACGTTTTTTGCTCGTCAGGAGTGGAGCTCGACGGCGATTACCGATCGTCCAACGGTGAATATGCTGGGGGGATATTACTCGCAGCAGCAATTTTTGCGTAACCTGGATGTCCGTTCAAACATGGCTTCTGCCGACCAACC  >  NC_000913/3969158‑3969311

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