Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,979,723	C→T	P375S (CCT→TCT)	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,723	0	C	T	100.0%	86.3 / NA	26	P375S (CCT→TCT)	wzyE	enterobacterial common antigen polymerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/12);  total (14/12)

TCCGCGATTTCTATGTCTTTATCCCTTCCTGGCTGTGGCCGGGTCGCCCGAGTATGGTGCTGAACTCAGCCAACTACTTTACCTGGGAAGTGCTGAATAACCACTCCGGACTGGCGATCTCGCCTACGCTTATAGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGCCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATAACAATGAATAACAACA  >  NC_000913/3979481‑3979968                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      tCCGCGATTTCTATGTCTTTATCCCTTCCTGGCTGTGGCCGGGTCGCCCGAGTATGGTGCTGAACTCAGCCAACTACTTTACCTGGGAAGTGCTGAATAACCACTCCGGACTGGCGATCTCGCCTACGCTTATAGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGc                                                                                                                                                                                                                                               >  2:180088/1‑251 (MQ=255)   cgcgCTTTCTATGTCTTTATCCCTTCCTGGCTGTGGCCGGGTCGCCCGAGTATGGTGCTGAACTCAGCCAACTACTTTACCTGGGAAGTGCTGAATAACCACTCCGGACTGGCGATCTCGCCTACGCTTATAGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGCTCGACTGGCTGGATGAGCTGGGCAACCGCGAGTCTAATCGCt                                                                                                                                                                                                                                              >  2:33342/1‑250 (MQ=255)                             cTGGCTGTGGCCGGGTCGCCCGAGTATGGTGCTGAACTCAGCCAACTACTTTACCTGGGAAGTGCTGAATAACCACTCCGGACTGGCGATCTCGCCTACGCTTATAGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCt                                                                                                                                                                                                                   >  1:38098/1‑251 (MQ=255)                                gCTGTGGCCGGGTCGCCCGAGTATGGTGCTGAACTCAGCCAACTACTTTACCTGGGACGTGCTGAATAACCACTCCGGACTGGCGATCTCGCCTACGCTTATAGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCt                                                                                                                                                                                                                >  2:212593/1‑251 (MQ=255)                                                 cGAGTATGGTGCTGAACTCAGCCAACTACTTTACCTGGGAAGTGCTGAATAACCACTCCGGACTGGCGATCTCGCCTACGCTTATAGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAat                                                                                                                                                                                               <  1:102930/251‑1 (MQ=255)                                                        ggTGCTGAACTCAGCCAACTACTTTACCTGGGAAGTGCTGAATAACCACTCCGGACTGGCGATCTCGCCTACGCTTATAGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATc                                                                                                                                                                                         <  1:219671/250‑1 (MQ=255)                                                                   cAGCCCACTCCTTTACCTGGGACGTGCTGAATAACCACTCCGGACTGGCGATCTCGCCTACGCTTATAGGCTCACTGGCGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGCCTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATCTTGCACAGTTTCTGCTTTGGTGCGATCTTCAATATGATCTTGCTGGCGCGt                                                                                                                                                                             >  2:354707/1‑251 (MQ=255)                                                                      ccAACTACTTTACCTGGGAAGTGCTGAATAACCACTCCGGACTGGCGATCTCGCCTACGCTTATAGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGaa                                                                                                                                                                          <  2:168272/251‑1 (MQ=255)                                                                          ctactTTACCTGGGAAGTGCTGAATAACCACTCCGGACTGGCGATCTCGCCTACGCTTATAGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGc                                                                                                                                                                      >  2:317727/1‑251 (MQ=255)                                                                                                aaTAACCACTCCGGACTGGCGATCTCGCCTACGCTTATAGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGt                                                                                                                                                <  1:105299/251‑1 (MQ=255)                                                                                                                         cGCCTACGCTTATAGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCt                                                                                                                            <  1:180088/246‑1 (MQ=255)                                                                                                                                    ataGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCtgat                                                                                                            <  1:159002/251‑1 (MQ=255)                                                                                                                                      aGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCtgatg                                                                                                           >  1:330537/1‑250 (MQ=255)                                                                                                                                       ggCTCACTGGGGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGTTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTgatgat                                                                                                         <  2:184648/251‑1 (MQ=255)                                                                                                                                                                          cccGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGcc                                                                       <  2:105043/250‑1 (MQ=255)                                                                                                                                                                           ccGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCgg                                                                     <  2:281579/251‑1 (MQ=255)                                                                                                                                                                                 ggggCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCGGACtcat                                                               >  2:193310/1‑251 (MQ=255)                                                                                                                                                                                       aTCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATaa                                                         <  1:371392/251‑1 (MQ=255)                                                                                                                                                                                       aTCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATa                                                          >  1:82906/1‑250 (MQ=255)                                                                                                                                                                                                            aaaTGGTTCGACTGGCTGTATGAGCTGGGCAAACGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCg                                    >  2:19993/1‑251 (MQ=255)                                                                                                                                                                                                                 gTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGc                               >  1:39851/1‑251 (MQ=255)                                                                                                                                                                                                                      aCTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGtt                          >  1:55957/1‑251 (MQ=255)                                                                                                                                                                                                                      aCTGGCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGtt                          >  2:50391/1‑251 (MQ=255)                                                                                                                                                                                                                         ggCTGTATGAGCTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCCCCGTTTCTGCTTGGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGGGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGGACAAAAGCATCGCTCCGGACGCATGTTGaa                       <  2:231394/251‑1 (MQ=255)                                                                                                                                                                                                                                    cTGGGCAACCGCGAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATAACAATg            >  1:16154/1‑251 (MQ=255)                                                                                                                                                                                                                                              gcgAGTCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATAACAATGAATAacaaca  <  2:56097/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      TCCGCGATTTCTATGTCTTTATCCCTTCCTGGCTGTGGCCGGGTCGCCCGAGTATGGTGCTGAACTCAGCCAACTACTTTACCTGGGAAGTGCTGAATAACCACTCCGGACTGGCGATCTCGCCTACGCTTATAGGCTCACTGGTGGTGATGGGCGGCGCGTTGTTCATCCCGCTCGGGGCGATCGTGGTTGGTCTGATCATCAAATGGTTCGACTGGCTGTATGAGCTGGGCAACCGCGAGCCTAATCGCTATAAAGCTGCGATATTGCACAGTTTCTGCTTTGGGGCGATCTTCAATATGATCGTGCTGGCGCGTGAAGGGCTGGATTCGTTTGTCTCACGCGTGGTCTTTTTTATCGTGGTCTTCGGCGCATGTCTGATGATCGCAAAACTGTTGTACTGGCTTTTTGAAAGCGCCGGACTCATTCATAAACGTACAAAATCATCGCTCCGGACGCAGGTTGAAGGATAACAATGAATAACAACA  >  NC_000913/3979481‑3979968

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