Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,970,062	C→T	R344R (CGC→CGT)	wzzE →	Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,970,062	0	C	T	100.0%	173.8 / NA	46	R344R (CGC→CGT)	wzzE	Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/24);  total (22/24)

CCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGA  >  NC_000913/3969916‑3970209                                                                                                                                                   |                                                                                                                                                    ccACCCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGtt                                                                                                                                                 <  2:577602/151‑1 (MQ=255) ccAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCcgtcgt                                                                                                                                                  <  2:342578/150‑1 (MQ=255)  cAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTg                                                                                                                                                <  1:93427/151‑1 (MQ=255)         ggATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATTATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGaaa                                                                                                                                         <  2:314134/151‑1 (MQ=255)         ggATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGa                                                                                                                                           <  2:75134/149‑1 (MQ=255)             ccGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGc                                                                                                                                     <  1:865327/151‑1 (MQ=255)                cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCaa                                                                                                                                   >  1:197771/1‑150 (MQ=255)                cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCaa                                                                                                                                   >  2:32165/1‑150 (MQ=255)                cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAac                                                                                                                                  <  1:669007/151‑1 (MQ=255)                           cTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGt                                                                                                                       >  2:243406/1‑151 (MQ=255)                                 cTATTTGCGGACGCCGGAAGAACCGGGAAAACGCGATAGCCCACGTCGTGACTTCCTGATGATGATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTg                                                                                                                      <  1:769421/146‑1 (MQ=255)                                  tATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAgcgc                                                                                                                 >  1:102119/1‑150 (MQ=255)                                     ttGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCaa                                                                                                               >  2:257784/1‑149 (MQ=255)                                      tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAg                                                                                                             <  1:547807/150‑1 (MQ=255)                                       gCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGgc                                                                                                           <  1:669377/151‑1 (MQ=255)                                               cGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                                                                                   >  1:341640/1‑151 (MQ=255)                                               cGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                                                                                   <  2:432712/151‑1 (MQ=255)                                                ggAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccgc                                                                                                  >  2:601961/1‑151 (MQ=255)                                                  aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCtt                                                                                                <  1:445204/151‑1 (MQ=255)                                                       ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCg                                                                                           >  1:708623/1‑151 (MQ=255)                                                             aaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGTCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAag                                                                                        >  2:861624/1‑148 (MQ=255)                                                              aaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAgaga                                                                                     >  1:229711/1‑150 (MQ=255)                                                               aCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                                                                                   <  2:908972/151‑1 (MQ=255)                                                                  cgATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCg                                                                                 >  2:86567/1‑150 (MQ=255)                                                                    aTAGCCCACGTCGTGCTTTCCTGATGATTATGTGGGGCATTGTCGGGGCGCCGATCGGGGCTGGTGTCGAATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGTGGGGAATGCTAAAGGGCTCGCCGCTTATTCGAAGAGAATCg                                                                                 >  2:446063/1‑148 (MQ=255)                                                                           aCGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGt                                                                       <  2:197746/151‑1 (MQ=255)                                                                             gtcgtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTac                                                                     <  2:551525/151‑1 (MQ=255)                                                                               cgtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg                                                                   >  1:695094/1‑151 (MQ=255)                                                                                  gCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgact                                                                <  2:199535/151‑1 (MQ=255)                                                                                     ttCCTGATGATTATGTGGTGCATTGGCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgactg                                                               <  2:185440/149‑1 (MQ=255)                                                                                       ccTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTAtt                                                           >  2:99401/1‑151 (MQ=255)                                                                                         tgatgaTTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTg                                                         >  1:254566/1‑151 (MQ=255)                                                                                           atgatTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGt                                                       <  1:604412/151‑1 (MQ=255)                                                                                             gatTAGGTGGGGCTTTGTCGGGGGGCTGATCGTTGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAAGAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTAc                                                     <  2:248678/151‑1 (MQ=255)                                                                                                          aTTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGcc                                        >  2:855921/1‑151 (MQ=255)                                                                                                            tGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCAt                                      >  2:345762/1‑151 (MQ=255)                                                                                                              tCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATc                                     >  2:65889/1‑150 (MQ=255)                                                                                                              tCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAACCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATc                                     >  1:87504/1‑150 (MQ=255)                                                                                                               cGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCaa                                   <  1:858232/151‑1 (MQ=255)                                                                                                               cGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCa                                    <  2:557334/150‑1 (MQ=255)                                                                                                                     gCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgc                             <  1:637136/151‑1 (MQ=255)                                                                                                                        gATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGt                                                                       >  1:412465/1‑106 (MQ=255)                                                                                                                        gATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGt                                                                       <  2:412436/106‑1 (MQ=255)                                                                                                                        gATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGcc                          >  1:355760/1‑151 (MQ=255)                                                                                                                                  ggTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCa                 <  2:447188/150‑1 (MQ=255)                                                                                                                                       cGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGtt           <  1:2274/151‑1 (MQ=255)                                                                                                                                               cccGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCaaa     >  2:899733/1‑149 (MQ=255)                                                                                                                                                ccgtcgTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCaaa     >  1:261421/1‑148 (MQ=255)                                                                                                                                                ccgtcgTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGa  <  2:708574/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAAGA  >  NC_000913/3969916‑3970209

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