Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,970,015	T→C	W329R (TGG→CGG)	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,015	0	T	C	100.0%	197.4 / NA	57	W329R (TGG→CGG)	wzzE	Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/28);  total (29/28)

CAGAATCGGGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑GGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTG  >  NC_000913/3969880‑3970148                                                                                                                                        |                                                                                                                                        cAGAATCGGGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCggg                                                                                                                                      <  2:907343/139‑1 (MQ=255)    aaTCGGGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCa                                                                                                                                   >  1:131111/1‑139 (MQ=255)         ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCa                                                                                                                                   <  2:25812/134‑1 (MQ=255)         ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCa                                                                                                                                   >  1:25812/1‑134 (MQ=255)         ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                                              >  2:382612/1‑139 (MQ=255)         ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                                              >  1:48366/1‑139 (MQ=255)         ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                                              >  2:251999/1‑139 (MQ=255)         ggCCATGTTAAACACCCTAAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                                              >  2:718698/1‑139 (MQ=255)           ccATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gg                                                                                                                           <  2:48366/139‑1 (MQ=255)                ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggggggg                                                                                                                      <  1:382612/139‑2 (MQ=255)                ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggggggg                                                                                                                      <  1:623920/139‑2 (MQ=255)                     cacCCTGAATGTTGGTCCAACCCTGGATCCGCGTTCTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGAt                                                                                                                 <  1:782117/139‑1 (MQ=255)                         cTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCggg                                                                                                             >  1:121477/1‑139 (MQ=255)                          tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCgggg                                                                                                            >  1:1210293/1‑139 (MQ=255)                                 tGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGgtgt                                                                                                     <  1:718698/139‑1 (MQ=255)                                  ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTc                                                                                                    >  1:901104/1‑139 (MQ=255)                                  ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gggggggctga                                                                                                                  <  2:40896/125‑6 (MQ=255)                                  ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gggggggctga                                                                                                                  >  1:40896/1‑120 (MQ=255)                                         cccTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggggggg                                                                                                                      <  1:1091272/114‑2 (MQ=255)                                         cccTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggggggg                                                                                                                      >  2:1091272/1‑113 (MQ=255)                                                   gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTg                                                                                   >  2:292376/1‑139 (MQ=255)                                                    cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                                                  >  2:673928/1‑139 (MQ=255)                                                     gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                 <  1:251999/139‑1 (MQ=255)                                                                tATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                      <  1:292376/139‑1 (MQ=255)                                                                 aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCa                                                                                                                                   >  1:847648/1‑78 (MQ=255)                                                                 aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCa                                                                                                                                   <  2:847648/78‑1 (MQ=255)                                                                 aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGc                                                                                                           >  2:311314/1‑101 (MQ=255)                                                                 aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGc                                                                                                           <  1:311314/101‑1 (MQ=255)                                                                 aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                              >  2:620748/1‑130 (MQ=255)                                                                 aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                              <  1:620748/130‑1 (MQ=255)                                                                  tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                    >  2:1006010/1‑139 (MQ=255)                                                                        tttGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCa                                                                                                                                   >  1:833143/1‑71 (MQ=255)                                                                        tttGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGAGCa                                                                                                                                   <  2:833143/71‑1 (MQ=255)                                                                            cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGt                                                          <  2:121477/139‑1 (MQ=255)                                                                                      aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                      >  2:1102690/1‑117 (MQ=255)                                                                                      aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                      <  1:1102690/117‑1 (MQ=255)                                                                                           ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gggggg                                                                                                                       <  2:257002/63‑1 (MQ=255)                                                                                           ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gggggg                                                                                                                       >  1:257002/1‑63 (MQ=255)                                                                                           ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gggggg                                                                                                                       <  2:714831/63‑1 (MQ=255)                                                                                           ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gggggg                                                                                                                       >  1:714831/1‑63 (MQ=255)                                                                                               tAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCAt                                                                                                                                  <  2:731492/49‑1 (MQ=255)                                                                                               tAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCAt                                                                                                                                  >  1:731492/1‑49 (MQ=255)                                                                                                aaaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                      <  2:1210293/139‑1 (MQ=255)                                                                                                                   cgtGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                   <  2:901104/139‑1 (MQ=255)                                                                                                                     tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACccgcc                                                                                       >  1:591630/1‑69 (MQ=255)                                                                                                                     tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACccgcc                                                                                       <  2:591630/69‑1 (MQ=255)                                                                                                                      gCCTTCCTGATGATTATGCGGGGCATTGTC‑GGGGGGACTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                <  2:861092/75‑1 (MQ=255)                                                                                                                      gCCTTCCTGATGATTATGCGGGGCATTGTC‑GGGGGGACTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                >  1:861092/1‑75 (MQ=255)                                                                                                                       ccTTCCTGTTGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                      <  2:1002169/68‑1 (MQ=255)                                                                                                                       ccTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                      >  1:1002169/1‑68 (MQ=255)                                                                                                                        cTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTTTTCGAAGAGAATCGATGTGa              <  1:673928/139‑1 (MQ=255)                                                                                                                        cTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGa              <  1:612089/139‑1 (MQ=255)                                                                                                                          tCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGaaa            >  2:424782/1‑139 (MQ=255)                                                                                                                                   ttATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgact   >  2:556106/1‑139 (MQ=255)                                                                                                                                    tATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                 <  2:643069/108‑1 (MQ=255)                                                                                                                                    tATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                 >  1:643069/1‑108 (MQ=255)                                                                                                                                    tATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgactg  >  2:295483/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        CAGAATCGGGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑GGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTG  >  NC_000913/3969880‑3970148

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