Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,970,077	G→A	*349* (TAG→TAA)	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,077	0	G	A	100.0%	198.4 / NA	62	*349* (TAG→TAA)	wzzE	Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/35);  total (27/35)

TCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTG‑TC‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAA  >  NC_000913/3969945‑3970205                                                                                                                                       |                                                                                                                                 tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTTCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                                                                             >  1:1168937/1‑139 (MQ=255) tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                                                                             >  1:670900/1‑139 (MQ=255) tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                                                                             <  2:601751/139‑1 (MQ=255)           cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGt                                                                                                                   >  2:741723/1‑139 (MQ=255)                     aagaaCCGGTAAAACGCGATAGCCCACGCCGTGCCTTCCTGATGATTATGCGGGGCATTGTTG‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAg                                                                                                         >  2:436546/1‑139 (MQ=255)                        aaCCGGGAAAACGCGATAGCCGACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTCCTGCGGTGAGCGCAAAGgcg                                                                                                      <  2:771448/139‑1 (MQ=255)                          ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                                                                               <  2:1040399/112‑1 (MQ=255)                          ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                                                                               >  1:1040399/1‑112 (MQ=255)                          ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACActg                                                                                                                          >  2:137566/1‑117 (MQ=255)                          ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACActg                                                                                                                          <  1:137566/117‑1 (MQ=255)                          ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt                                                                                                    >  1:954555/1‑139 (MQ=255)                           cGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTc                                                                                                   <  1:436546/139‑1 (MQ=255)                               aaaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                                                                               <  2:643918/139‑1 (MQ=255)                                    gcgATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                                                                          <  2:1296764/139‑1 (MQ=255)                                       aTAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCg                                                                                       >  2:1110304/1‑139 (MQ=255)                                           cccACGTCGTGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAaga                                                                                   >  2:407989/1‑139 (MQ=255)                                                  cgtGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGTTCGGGGTTGGTGTCGCATTAACCCGCCGTTGCTCGAAAGAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                                                                            <  2:308496/139‑1 (MQ=255)                                                  cgtGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                                                                            <  1:1138521/139‑1 (MQ=255)                                                  cgtGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                                                                            <  2:1069538/139‑1 (MQ=255)                                                  cgtGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                                                                            <  1:1110304/139‑1 (MQ=255)                                                  cgtGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                                                                            <  2:954555/139‑1 (MQ=255)                                                    tGCCTTCCTGATGATTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                                                                          >  2:542137/1‑139 (MQ=255)                                                              atgatTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTact                                                                <  2:670900/139‑1 (MQ=255)                                                                gatTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                                              <  2:1168937/139‑1 (MQ=255)                                                                 atTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                                             <  1:542137/139‑1 (MQ=255)                                                                 atTATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                                             <  2:187832/139‑1 (MQ=255)                                                                   tATGCGGGGCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgactg                                                           >  2:272788/1‑139 (MQ=255)                                                                     tgcggggCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTa                                                         >  1:44516/4‑139 (MQ=255)                                                                        ggggCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCg                                                                             <  1:298953/116‑1 (MQ=255)                                                                        ggggCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCg                                                                             >  2:298953/1‑116 (MQ=255)                                                                          ggCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATc                                                                              >  1:1261191/1‑113 (MQ=255)                                                                          ggCATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATc                                                                              <  2:1261191/113‑1 (MQ=255)                                                                            cATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                                                                          <  1:283076/99‑1 (MQ=255)                                                                            cATTG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                                                                          >  2:283076/1‑99 (MQ=255)                                                                              ttG‑TCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTAcg                                                <  1:768237/139‑1 (MQ=255)                                                                                tgtcgggggggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                                                                          <  1:930354/107‑1 (MQ=255)                                                                                tgtcgggggggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTAGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                                                                          >  2:930354/6‑112 (MQ=255)                                                                                         ggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                                             <  2:724342/116‑1 (MQ=255)                                                                                         ggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                                             >  1:724342/1‑116 (MQ=255)                                                                                         ggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAg                                      >  1:907148/1‑139 (MQ=255)                                                                                            tGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAtt                                                                                         <  1:303217/85‑1 (MQ=255)                                                                                            tGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAtt                                                                                         >  2:303217/1‑85 (MQ=255)                                                                                                      tGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGcc                                            >  2:563692/1‑120 (MQ=255)                                                                                                      tGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGcc                                            <  1:563692/120‑1 (MQ=255)                                                                                                       ggTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcg                        <  2:44516/139‑1 (MQ=255)                                                                                                        gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                                             >  2:727017/1‑101 (MQ=255)                                                                                                        gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                                             <  1:727017/101‑1 (MQ=255)                                                                                                        gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                                              <  2:146446/100‑1 (MQ=255)                                                                                                        gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                                              >  1:146446/1‑100 (MQ=255)                                                                                                        gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcgc                       >  1:34829/1‑139 (MQ=255)                                                                                                        gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcg                        >  1:708522/1‑138 (MQ=255)                                                                                                        gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcg                        <  2:708522/138‑1 (MQ=255)                                                                                                          gtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCg                     <  2:34829/139‑1 (MQ=255)                                                                                                           tCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGt                    <  2:552209/139‑1 (MQ=255)                                                                                                             gCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                                             <  2:1005935/96‑1 (MQ=255)                                                                                                             gCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                                             >  1:1005935/1‑96 (MQ=255)                                                                                                             gCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTact                                                                <  1:1008316/93‑1 (MQ=255)                                                                                                             gCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTact                                                                >  2:1008316/1‑93 (MQ=255)                                                                                                                      cgccgTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCg         <  2:907148/139‑1 (MQ=255)                                                                                                                       gccgTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGt        <  1:272788/139‑1 (MQ=255)                                                                                                                        ccgTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                                                                            >  1:970669/1‑70 (MQ=255)                                                                                                                        ccgTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                                                                            <  2:970669/70‑1 (MQ=255)                                                                                                                           ttGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGa                                                                       <  1:985518/72‑1 (MQ=255)                                                                                                                           ttGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGa                                                                       >  2:985518/1‑72 (MQ=255)                                                                                                                             gCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCaa  >  1:1066893/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                 TCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTG‑TC‑GGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAA  >  NC_000913/3969945‑3970205

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