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%	125.9 / NA	36	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 (17/19);  total (17/19)

GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑GGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATG  >  NC_000913/3969888‑3970133                                                                                                                                |                                                                                                                         ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCAt                                                                                                                   <  2:116250/135‑1 (MQ=255) ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCAt                                                                                                                   >  1:116250/1‑135 (MQ=255) ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                               >  1:376119/1‑139 (MQ=255)        ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggggg‑gg                                                                                                      <  1:601821/139‑2 (MQ=255)        ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggggg‑gg                                                                                                      <  2:247528/139‑2 (MQ=255)                  tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCgggg                                                                                             <  1:379252/139‑1 (MQ=255)                          ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTc                                                                                     >  2:311350/1‑139 (MQ=255)                              cAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCAt                                                                                 <  1:672270/139‑1 (MQ=255)                                            cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGTCATTGTCGGGGGG‑GCTGATCGGGGc                                                                                            <  1:646646/114‑1 (MQ=255)                                            cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGc                                                                                            >  2:646646/1‑114 (MQ=255)                                                        tATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                       <  1:311350/139‑1 (MQ=255)                                                          tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgc                                                                         >  2:261495/1‑119 (MQ=255)                                                          tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgc                                                                         <  1:261495/119‑1 (MQ=255)                                                          tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                     >  1:50505/1‑139 (MQ=255)                                                          tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                     <  2:376119/139‑1 (MQ=255)                                                                    cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTg                                                                                          <  2:550649/92‑1 (MQ=255)                                                                    cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTg                                                                                          >  1:550649/1‑92 (MQ=255)                                                                    cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                       >  1:1052938/1‑111 (MQ=255)                                                                    cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGt                                           >  1:564446/1‑139 (MQ=255)                                                                    cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCGGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                       <  2:1052938/111‑1 (MQ=255)                                                                        cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCAt                                                                                 >  2:720165/1‑97 (MQ=255)                                                                        cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCAt                                                                                 <  1:720165/97‑1 (MQ=255)                                                                        cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                               <  2:523168/115‑1 (MQ=255)                                                                        cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                               >  1:523168/1‑115 (MQ=255)                                                                           cGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCa                                    <  1:55420/139‑1 (MQ=255)                                                                                                             tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                       >  1:835136/1‑70 (MQ=255)                                                                                                             tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                       <  2:413601/70‑1 (MQ=255)                                                                                                             tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                       >  1:413601/1‑70 (MQ=255)                                                                                                             tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                       <  2:835136/70‑1 (MQ=255)                                                                                                             tGCCTTCCTGATGATTATGCGGGGCATTGTC‑GGGGGAGCTGATCGGGACTGGTGTCGCATTAAGCCGCCGTTGCTCGAAATAACAACGCTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg  >  1:722650/1‑139 (MQ=255)                                                                                                               ccTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                 >  2:270852/1‑74 (MQ=255)                                                                                                               ccTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                 <  1:270852/74‑1 (MQ=255)                                                                                                                   ccTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                       <  1:692804/64‑1 (MQ=255)                                                                                                                   ccTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                       >  2:692804/1‑64 (MQ=255)                                                                                                                            tATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCg                                             >  2:1024154/1‑81 (MQ=255)                                                                                                                            tATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCg                                             <  1:1024154/81‑1 (MQ=255)                                                                                                                                |                                                                                                                         GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑GGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATG  >  NC_000913/3969888‑3970133

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: ‑