Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,970,062	C→T	R344R (CGC→CGT)	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,062	0	C	T	100.0%	161.3 / NA	44	R344R (CGC→CGT)	wzzE	Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/23);  total (21/23)

TGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAA  >  NC_000913/3969923‑3970207                                                                                                                                            |                                                                                                                                                  tGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGa                                                                                                                                         >  2:528122/1‑150 (MQ=255)  ggATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGaaa                                                                                                                                       <  1:84668/151‑1 (MQ=255)     tCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAAt                                                                                                                                      >  1:705569/1‑149 (MQ=255)      ccGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGc                                                                                                                                   >  1:334376/1‑151 (MQ=255)        gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCa                                                                                                                                  >  1:538954/1‑150 (MQ=255)             ttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACAct                                                                                                                             <  2:425621/150‑1 (MQ=255)              tCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACActg                                                                                                                            <  1:265589/150‑1 (MQ=255)              tCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGT‑GGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACActgc                                                                                                                           <  2:128838/150‑1 (MQ=255)                   ccTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCgg                                                                                                                      >  1:763377/1‑151 (MQ=255)                      aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTg                                                                                                                    <  1:537434/150‑1 (MQ=255)                      aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGa                                                                                                                   <  1:711830/151‑1 (MQ=255)                          cTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAgcgc                                                                                                               >  2:137818/1‑151 (MQ=255)                            aTTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAgcgc                                                                                                               <  2:2846/149‑1 (MQ=255)                               tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAgg                                                                                                          >  1:575671/1‑151 (MQ=255)                               tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCAATGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAgg                                                                                                          >  2:453688/1‑151 (MQ=255)                                gCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGgc                                                                                                         <  2:520138/151‑1 (MQ=255)                                    ccGCCGGAAGAACCGGAAAAACGCGATAGCCCACGACATGCCTTCCTGATGATAATGTGGGGCATTGTCGGGGGGCTGATCGGGGTTGGTGTCGAATTAACCCCTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCACAGGCGCTc                                                                                                     >  2:107710/2‑151 (MQ=255)                                    aCGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCt                                                                                                      <  2:583305/150‑1 (MQ=255)                                    aCGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCt                                                                                                      <  1:7853/150‑1 (MQ=255)                                    aCGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTc                                                                                                     <  1:46657/151‑1 (MQ=255)                                      gCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTcgc                                                                                                   >  1:509705/1‑151 (MQ=255)                                          gaagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                                                                               >  2:468952/1‑151 (MQ=255)                                           aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCtt                                                                                              <  1:263598/151‑1 (MQ=255)                                             gaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                                                                            <  1:541272/151‑1 (MQ=255)                                                ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCg                                                                                         <  1:743187/151‑1 (MQ=255)                                                  ggTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGa                                                                                        >  2:732129/1‑150 (MQ=255)                                                          gcgATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCg                                                                               >  1:374107/1‑151 (MQ=255)                                                                      gtcgtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTac                                                                   >  1:133158/1‑151 (MQ=255)                                                                        cgtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg                                                                 <  1:27/151‑1 (MQ=255)                                                                         gtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga                                                                <  1:287976/151‑1 (MQ=255)                                                                         gtGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg                                                                 >  2:391898/1‑150 (MQ=255)                                                                           gCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgact                                                              <  2:711679/151‑1 (MQ=255)                                                                                       atTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTAcg                                                  >  1:594505/1‑151 (MQ=255)                                                                                               gggCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCgg                                           >  1:185947/1‑150 (MQ=255)                                                                                                   aTTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGc                                       <  1:453842/150‑1 (MQ=255)                                                                                                     tGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCa                                     <  1:447162/150‑1 (MQ=255)                                                                                                          gggggCTGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCa                                  >  2:639065/1‑148 (MQ=255)                                                                                                                tGATCGGGGCTGGTGTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGgcgc                         >  2:448226/1‑151 (MQ=255)                                                                                                                           ggTCTCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGAAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCAt              <  1:492623/151‑1 (MQ=255)                                                                                                                              gtCGCATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGc            <  1:468901/150‑1 (MQ=255)                                                                                                                                  cATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTgg       <  2:512869/151‑1 (MQ=255)                                                                                                                                  cATTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTgg       >  2:765106/1‑151 (MQ=255)                                                                                                                                   aTTAACCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGc      >  2:65235/1‑151 (MQ=255)                                                                                                                                       aCCCGTCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCaaaa  <  2:15411/151‑1 (MQ=255)                                                                                                                                         ccgtcgTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCaaa   >  2:715838/1‑148 (MQ=255)                                                                                                                                            |                                                                                                                                                  TGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACGCGCCCGGAAGCCATCAAGATGGCGCCGTTGGTGCATGCGTTGGCAAAA  >  NC_000913/3969923‑3970207

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