Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,969,685	A→G	I219V (ATC→GTC)	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,969,685	0	A	G	100.0%	175.6 / NA	50	I219V (ATC→GTC)	wzzE	Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/26);  total (24/26)

ACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCG  >  NC_000913/3969539‑3969830                                                                                                                                                   |                                                                                                                                                  aCGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTa                                                                                                                                               <  1:470783/151‑1 (MQ=255)   gCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACg                                                                                                                                             >  1:67977/1‑151 (MQ=255)      aataaCCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGAcc                                                                                                                                          >  1:195682/1‑151 (MQ=255)      aataaCCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGAcc                                                                                                                                          <  2:704896/151‑1 (MQ=255)        taaCCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGAccgc                                                                                                                                        >  1:237717/1‑151 (MQ=255)              gTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATg                                                                                                                                  >  2:481523/1‑151 (MQ=255)                  cGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACa                                                                                                                              <  2:66122/151‑1 (MQ=255)                   gTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAAc                                                                                                                               >  1:742807/1‑149 (MQ=255)                       gTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCAtt                                                                                                                         <  2:766559/151‑1 (MQ=255)                        tATGGTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCAtt                                                                                                                         <  1:761701/150‑1 (MQ=255)                         aTGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTg                                                                                                                        >  1:512817/1‑150 (MQ=255)                          tGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACAGCCGGATGAACAGCATTg                                                                                                                        >  1:634833/1‑149 (MQ=255)                          tGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGGGAAGCGTCAGGAAGCGGTGGCGAAAGCCGTCTACGACCGCAGGATGAACAGCATTg                                                                                                                        >  2:492675/1‑149 (MQ=255)                           gTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAg                                                                                                                      <  2:590300/150‑1 (MQ=255)                           gTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAg                                                                                                                      <  1:631620/150‑1 (MQ=255)                             tGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGAGTACCATCCAGATGAAAGCTCAGGAGAAGAGTCAGGAAGAGGTGGCGAAAGCCGTCTTCGACCGCCGGATGAACAGCATTGAg                                                                                                                      >  1:711605/1‑148 (MQ=255)                                  ttGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCt                                                                                                              >  2:188899/1‑151 (MQ=255)                                   tGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCt                                                                                                              <  2:298539/150‑1 (MQ=255)                                    gccagccaGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTg                                                                                                             >  2:64889/1‑150 (MQ=255)                                       agccagcGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGaaaa                                                                                                         <  2:541098/151‑1 (MQ=255)                                       agccagcGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGaaa                                                                                                          >  2:679612/1‑150 (MQ=255)                                        gccagcGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAAt                                                                                                        <  1:447691/151‑1 (MQ=255)                                                  agccagccATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGATGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGagcag                                                                                              <  1:44486/151‑1 (MQ=255)                                                   gccagccaTCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGagcagc                                                                                             <  2:190876/151‑1 (MQ=255)                                                    ccagccaTCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGagcagc                                                                                             >  1:388692/1‑150 (MQ=255)                                                     cagccaTCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCat                                                                                           <  1:64858/151‑1 (MQ=255)                                                         caTCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCataat                                                                                        >  2:743262/1‑150 (MQ=255)                                                              gAATGATGAGCCGAAAGGCGCAGGGGTGGCGAGTAACATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGACGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAatat                                                                                      <  1:310703/147‑1 (MQ=255)                                                                                   aTGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTTAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCt                                                             >  2:515736/1‑151 (MQ=255)                                                                                            gcgTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGaa                                                    <  2:726998/151‑1 (MQ=255)                                                                                            gcgTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGa                                                     <  2:278459/150‑1 (MQ=255)                                                                                               tACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAAtt                                                  >  2:265622/1‑150 (MQ=255)                                                                                                           gAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGaa                                     <  2:698517/151‑1 (MQ=255)                                                                                                           gAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGa                                      <  2:430926/150‑1 (MQ=255)                                                                                                           gAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGa                                      >  2:64067/1‑150 (MQ=255)                                                                                                            aaaGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGaa                                     <  1:87335/150‑1 (MQ=255)                                                                                                               gCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGt                                 <  2:770649/151‑1 (MQ=255)                                                                                                                  cAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGtt                                <  2:383579/149‑1 (MQ=255)                                                                                                                      tGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCt                          >  2:129073/1‑151 (MQ=255)                                                                                                                       gAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCt                          <  2:128970/150‑1 (MQ=255)                                                                                                                         aGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTgg                       >  2:46333/1‑151 (MQ=255)                                                                                                                            gTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCg                    >  1:164871/1‑151 (MQ=255)                                                                                                                              cAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTc                  <  2:103175/151‑1 (MQ=255)                                                                                                                                     agGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCt           <  2:732947/151‑1 (MQ=255)                                                                                                                                      ggTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTCCCTGCCGAGGAATTACCTGATTCAGAAGTGTTCCTGCTTGGGCGTCCAATGCt           <  1:286224/150‑1 (MQ=255)                                                                                                                                          gCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAgg      <  2:122011/151‑1 (MQ=255)                                                                                                                                            gAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCt    >  2:75809/1‑151 (MQ=255)                                                                                                                                             aaaGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTc   >  2:600478/1‑151 (MQ=255)                                                                                                                                              aaGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCg  >  2:561368/1‑151 (MQ=255)                                                                                                                                              aaGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCg  >  1:2846/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                  ACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCG  >  NC_000913/3969539‑3969830

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