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%	180.1 / NA	56	*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 (25/31);  total (25/31)

ATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑G‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGC  >  NC_000913/3969944‑3970193                                                                                                                                        |                                                                                                                     aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATaacaac                                                                                                                  >  2:1133442/1‑139 (MQ=255)                     gaagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCaaa                                                                                              >  1:1126110/1‑139 (MQ=255)                                  aaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                                 <  1:1133442/139‑1 (MQ=255)                                       gATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTc                                                                            <  2:100274/139‑1 (MQ=255)                                               aCGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGaa                                                                    >  1:363347/1‑139 (MQ=255)                                                cgtcgtGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAAGAACAACACTGCTGCg                                                                                                         <  2:956820/101‑1 (MQ=255)                                                cgtcgtGCCTTCCTGATGATTATGCGGGGCATTGTC‑GCGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTCCg                                                                                                         >  1:956820/1‑101 (MQ=255)                                                     tGCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                                                              <  2:519292/139‑1 (MQ=255)                                                      gCCTTCCTGATGATTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgt                                                             >  1:1010364/1‑139 (MQ=255)                                                             tgatgaTTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                                 <  1:450849/112‑1 (MQ=255)                                                             tgatgaTTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                                 >  2:450849/1‑112 (MQ=255)                                                               atgatTATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTact                                                    <  1:78545/139‑1 (MQ=255)                                                                   ttATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgact                                                <  2:1126110/139‑1 (MQ=255)                                                                    tATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                                 <  1:249723/105‑1 (MQ=255)                                                                    tATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                                 >  2:249723/1‑105 (MQ=255)                                                                    tATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                                 >  2:1032646/1‑137 (MQ=255)                                                                    tATGCGGGGCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                                 <  1:1032646/137‑1 (MQ=255)                                                                      tgcggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAtt                                                                             >  2:369240/4‑107 (MQ=255)                                                                      tgcggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAtt                                                                             <  1:369240/104‑1 (MQ=255)                                                                      tgcggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                                  >  1:654244/4‑134 (MQ=255)                                                                      tgcggggCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                                  <  2:654244/131‑1 (MQ=255)                                                                            gCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGaa                                                          <  2:1358606/120‑1 (MQ=255)                                                                            gCATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGaa                                                          >  1:1358606/1‑120 (MQ=255)                                                                             cATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt                                                                                        <  1:414046/89‑1 (MQ=255)                                                                             cATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt                                                                                        >  2:414046/1‑89 (MQ=255)                                                                             cATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                                                                   <  2:1429048/110‑1 (MQ=255)                                                                             cATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                                                                   >  1:1429048/1‑110 (MQ=255)                                                                             cATTGTCGG‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTa                                      <  1:1416976/139‑1 (MQ=255)                                                                                        ggggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                                 <  1:1295337/86‑1 (MQ=255)                                                                                        ggggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                                 >  2:1295337/1‑86 (MQ=255)                                                                                         gggCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGCGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGaa                           <  2:1010364/139‑1 (MQ=255)                                                                                         gggCTGATCGGGGCTAGTGTCGCATTAACCCGCCGTTGCTCGAAATAACCACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGaa                           <  2:1199468/139‑1 (MQ=255)                                                                                                   gggCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGt                                                                                                       <  1:1194081/53‑1 (MQ=255)                                                                                                   gggCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGt                                                                                                       >  2:1194081/1‑53 (MQ=255)                                                                                                    ggCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                                                                   >  1:244899/1‑88 (MQ=255)                                                                                                    ggCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                                                                   <  2:244899/88‑1 (MQ=255)                                                                                                      cTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATgg              <  2:363347/139‑1 (MQ=255)                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                                                                   <  2:865648/35‑1 (MQ=255)                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                                                                   >  1:865648/1‑35 (MQ=255)                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAgg                                                                                            <  2:1104931/58‑1 (MQ=255)                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAgg                                                                                            >  1:1104931/1‑58 (MQ=255)                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                                 <  2:510882/69‑1 (MQ=255)                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                                 >  1:510882/1‑69 (MQ=255)                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                                                                   <  2:29841/83‑1 (MQ=255)                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                                                                   >  1:29841/1‑83 (MQ=255)                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGAt                <  2:1400975/134‑1 (MQ=255)                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGAt                >  1:1400975/1‑134 (MQ=255)                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcg            >  1:666996/1‑138 (MQ=255)                                                                                                         gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcg            <  2:666996/138‑1 (MQ=255)                                                                                                          tgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgtg                                                            >  1:560903/1‑89 (MQ=255)                                                                                                          tgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgtg                                                            <  2:560903/89‑1 (MQ=255)                                                                                                                    aCCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                                                                   <  1:1101914/56‑1 (MQ=255)                                                                                                                    aCCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                                                                   >  2:1101914/1‑56 (MQ=255)                                                                                                                            ttGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTc                                                                                       <  1:56299/44‑1 (MQ=255)                                                                                                                            ttGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTc                                                                                       >  2:56299/1‑44 (MQ=255)                                                                                                                            ttGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGc  <  1:499076/129‑1 (MQ=255)                                                                                                                            ttGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGc  >  2:499076/1‑129 (MQ=255)                                                                                                                                        |                                                                                                                     ATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑G‑GGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGC  >  NC_000913/3969944‑3970193

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