Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,969,685	A→G	I219V (ATC→GTC)	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,969,685	0	A	G	100.0%	173.9 / NA	49	I219V (ATC→GTC)	wzzE	Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/29);  total (20/29)

ACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCG  >  NC_000913/3969539‑3969830                                                                                                                                                   |                                                                                                                                                  aCGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCt                                                                                                                                                >  2:706184/1‑150 (MQ=255)  cGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTAc                                                                                                                                              >  2:93405/1‑151 (MQ=255)   gCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACg                                                                                                                                             <  2:438991/151‑1 (MQ=255)   gCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACg                                                                                                                                             <  2:889885/151‑1 (MQ=255)            cTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCgg                                                                                                                                     >  2:344692/1‑150 (MQ=255)              gTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATg                                                                                                                                  <  1:264874/151‑1 (MQ=255)                     cAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCa                                                                                                                           <  1:317496/151‑1 (MQ=255)                       gTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCAt                                                                                                                          <  1:289696/150‑1 (MQ=255)                             tGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCa                                                                                                                    <  1:749060/150‑1 (MQ=255)                              gCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAgg                                                                                                                  >  2:669329/1‑151 (MQ=255)                               cTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGgc                                                                                                                 >  1:870321/1‑151 (MQ=255)                                 tttGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGgc                                                                                                                 <  1:486078/149‑1 (MQ=255)                                               tGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTg                                                                                                   >  1:855919/1‑149 (MQ=255)                                                 cagccagcCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGagca                                                                                               <  1:409276/151‑1 (MQ=255)                                                  agccagccATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGagcag                                                                                              <  2:329838/151‑1 (MQ=255)                                                     cagccaTCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCat                                                                                           <  1:140704/151‑1 (MQ=255)                                                     cagccaTCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCat                                                                                           >  1:909045/1‑151 (MQ=255)                                                              gAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATAttt                                                                                    <  2:730981/149‑1 (MQ=255)                                                               aaTGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTcgc                                                                                 >  1:401311/1‑151 (MQ=255)                                                                  gatgaGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTcgcg                                                                                <  1:546055/149‑1 (MQ=255)                                                                       gCTGAAAGGCGCATGGGCCGCGCGGACCAGCCAGGTGAAAGCTCAGGTGAAGCTTCAGGAAGTGGTGGCGAAGGCCGTTTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGGGCAGCATATGATTGCGTGCAGTGc                                                                          <  2:520490/150‑1 (MQ=255)                                                                          gAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGAc                                                                       >  2:690330/1‑150 (MQ=255)                                                                          gAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACa                                                                      <  1:591106/151‑1 (MQ=255)                                                                           aaaGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAg                                                                     <  2:556195/151‑1 (MQ=255)                                                                             aGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGAt                                                                   <  1:434450/151‑1 (MQ=255)                                                                             aGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGAt                                                                   <  1:207392/151‑1 (MQ=255)                                                                              ggCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATg                                                                  <  1:603051/151‑1 (MQ=255)                                                                                cgcATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGt                                                                 >  2:87760/1‑150 (MQ=255)                                                                                 gcATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTa                                                                <  1:802210/150‑1 (MQ=255)                                                                                   aTGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCt                                                             <  2:469759/151‑1 (MQ=255)                                                                                     gggCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTg                                                            <  2:70825/150‑1 (MQ=255)                                                                                         ggcgCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCg                                                         >  1:385055/1‑149 (MQ=255)                                                                                            gcgTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGa                                                     <  1:1184/150‑1 (MQ=255)                                                                                                 ccATCCAGGTGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTAcc                                               >  1:342604/1‑151 (MQ=255)                                                                                                    tCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTg                                             >  1:644064/1‑150 (MQ=255)                                                                                                     ccAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGa                                            <  2:123191/150‑1 (MQ=255)                                                                                                     ccAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGAt                                           >  2:668959/1‑151 (MQ=255)                                                                                                            aaaGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGaaa                                    <  2:114840/151‑1 (MQ=255)                                                                                                             aaGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAAt                                   <  1:617914/151‑1 (MQ=255)                                                                                                              aGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATg                                  >  1:507427/1‑151 (MQ=255)                                                                                                                  cAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTc                               <  2:48739/150‑1 (MQ=255)                                                                                                                   aGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCt                             >  1:676197/1‑151 (MQ=255)                                                                                                                         aGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTgg                       <  1:159941/151‑1 (MQ=255)                                                                                                                         aGCGTCAGGAAGAGGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTgg                       <  1:641246/151‑1 (MQ=255)                                                                                                                            gTCAGGAAGAGGTGGCGACATCCGTCTACGACCGCCTTATTAACAGCATTGAGCAGTCGCTTAAAATTTCTGAGCATCATAATATTTCGCGCAGTGCGACAGATTTACCTGCCGAGTAATTACCTGATTCAGAAATTTTCCTGCtttttcg                    >  1:358556/1‑146 (MQ=255)                                                                                                                                   agagGTGGCGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGTAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATg             <  2:891091/151‑1 (MQ=255)                                                                                                                                           cGAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGGAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTAACTGCCGAGGAATTACCTGATTAAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAgg      >  2:762128/1‑150 (MQ=255)                                                                                                                                            gAAAGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCt    >  1:701062/1‑151 (MQ=255)                                                                                                                                              aaGCCGTCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCg  >  2:176114/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                  ACGCTAATAACCTGTTACGTCAGTATGTTGCTTTTGCCAGCCAGCGTGCAGCCAGCCATCTGAATGATGAGCTGAAAGGCGCATGGGCGGCGCGTACCATCCAGATGAAAGCTCAGGTGAAGCGTCAGGAAGAGGTGGCGAAAGCCATCTACGACCGCCGGATGAACAGCATTGAGCAGGCGCTGAAAATTGCTGAGCAGCATAATATTTCGCGCAGTGCGACAGATGTACCTGCCGAGGAATTACCTGATTCAGAAATGTTCCTGCTTGGGCGTCCAATGCTTCAGGCTCG  >  NC_000913/3969539‑3969830

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