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%	136.6 / NA	43	*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 (23/20);  total (23/20)

ATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCA‑TTGTC‑GGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGC  >  NC_000913/3969944‑3970184                                                                                                                                         |                                                                                                            aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCA‑TTGTCGGGGGGG‑GCGATCGGGGCTGGTGTCGCCTTAACCCCCCCTTGGTCGGAAAaacaac                                                                                                         >  1:962267/1‑139 (MQ=255) aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATaacaac                                                                                                         >  2:489795/1‑139 (MQ=255)  tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                                                        <  2:539996/139‑1 (MQ=255)             gTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTg                                                                                             >  2:352114/1‑139 (MQ=255)               aCGGCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCA‑TTGTC‑GGGGGGACTGATCGGGGCTGGTGTAGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAg                                                                                           >  2:1065479/1‑139 (MQ=255)               aCGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAg                                                                                           >  1:1174937/1‑139 (MQ=255)                                  aaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                        <  2:1258696/139‑1 (MQ=255)                                     gcgATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCA‑TTGTCGGGGGGG‑CTGATTGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                                                     <  2:1003616/139‑1 (MQ=255)                                               aCGTCGTGCCTTCCTGATGATTATGCGGGGCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGaa                                                           >  2:840865/1‑139 (MQ=255)                                                   cgtGCCTTCCTGATGATTATGCGGGGCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCCGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                                                       <  2:1174937/139‑1 (MQ=255)                                                     tGCCTTCCTGATGATTATGCGGGGCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                                                     >  1:104588/1‑139 (MQ=255)                                                         ttCCTGATGATTATGTGGGGCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGaa                                                 <  1:352114/139‑1 (MQ=255)                                                          tCCTGATGATTATGCGGGGCA‑TTGGCGGGGGGG‑CTGATCGGGGCGGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGaaa                                                <  1:994763/139‑1 (MQ=255)                                                             tgatgaTTATGCGGGGCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTa                                             <  1:1065479/139‑1 (MQ=255)                                                                  atTATGCGGGGCATTTGCG‑GGGGGG‑CTGATCGGGGCTGGGGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                        <  2:170220/139‑1 (MQ=255)                                                                    tATGCGGGGCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                                                     >  2:488141/1‑124 (MQ=255)                                                                    tATGCGGGGCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                                                     <  1:488141/124‑1 (MQ=255)                                                                      tgcggggCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTa                                    >  2:12034/4‑139 (MQ=255)                                                                           ggCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                                                          <  1:72193/64‑1 (MQ=255)                                                                           ggCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                                                          >  2:72193/1‑64 (MQ=255)                                                                           ggCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                                                          >  1:1118588/1‑96 (MQ=255)                                                                           ggCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                                                          <  2:1118588/96‑1 (MQ=255)                                                                           ggCA‑TTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTgg                               >  1:1310090/1‑139 (MQ=255)                                                                                 tgtcggggggg‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                        <  2:1279110/88‑1 (MQ=255)                                                                                 tgtcggggggg‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                        >  1:1279110/6‑93 (MQ=255)                                                                                     ggggggg‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGccc                      >  2:523223/2‑139 (MQ=255)                                                                                      gggggg‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg                                          >  2:1634774/1‑118 (MQ=255)                                                                                      gggggg‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg                                          <  1:1046444/118‑1 (MQ=255)                                                                                      gggggg‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg                                          <  1:1634774/118‑1 (MQ=255)                                                                                      gggggg‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGAGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg                                          >  2:1046444/1‑118 (MQ=255)                                                                                       ggggg‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCgg                    <  1:489795/139‑1 (MQ=255)                                                                                           g‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGc                >  1:1440329/1‑139 (MQ=255)                                                                                              tGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                                                       <  1:1307595/98‑1 (MQ=255)                                                                                              tGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                                                       >  2:1307595/1‑98 (MQ=255)                                                                                                     ggCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                                                     <  2:1614426/93‑1 (MQ=255)                                                                                                     ggCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                                                     >  1:1614426/1‑93 (MQ=255)                                                                                                          gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                                                          >  1:194376/1‑83 (MQ=255)                                                                                                          gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                                                          <  2:194376/83‑1 (MQ=255)                                                                                                          gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                        >  1:1251362/1‑101 (MQ=255)                                                                                                          gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgac                                        <  2:1251362/101‑1 (MQ=255)                                                                                                          gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcgc  >  1:1090525/1‑139 (MQ=255)                                                                                                          gtgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcgc  >  1:294478/1‑139 (MQ=255)                                                                                                           tgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGgcg                                                                                 <  1:968765/59‑1 (MQ=255)                                                                                                           tgtCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGgcg                                                                                 >  2:968765/1‑59 (MQ=255)                                                                                                                     aCCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                         >  2:1028152/1‑89 (MQ=255)                                                                                                                     aCCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                         <  1:1028152/89‑1 (MQ=255)                                                                                                                                         |                                                                                                            ATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCA‑TTGTC‑GGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGC  >  NC_000913/3969944‑3970184

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