Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	417,375	C→T	8.3%	P361L (CCG→CTG)	waaG →	lipopolysaccharide glucosyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	417,375	0	C	T	8.3%	173.4 / 5.2	60	P361L (CCG→CTG)	waaG	lipopolysaccharide glucosyltransferase I
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base T (2/3);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCCGACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCG  >  NC_002947/417227‑417521                                                                                                                                                     |                                                                                                                                                    gCCGACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCg                                                                                                                                                    >  1:774148/1‑150 (MQ=255)   cGACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCgca                                                                                                                                                  >  1:693200/1‑150 (MQ=255)           gCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTgcggccgtctg                                                                                                                                          <  2:599961/150‑7 (MQ=255)           gCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCAtgctg                                                                                                                                          <  1:758588/150‑1 (MQ=255)               ggTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGa                                                                                                                                      <  2:875616/150‑1 (MQ=255)                  gCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGt                                                                                                                                   >  1:395393/1‑150 (MQ=255)                    tGGAAGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGGGCAGGCGCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCGTTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGa                                                                                                                                 >  2:676526/1‑150 (MQ=255)                      gACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATc                                                                                                                               >  2:1140264/1‑150 (MQ=255)                      gACGAGCAGTCAGAACAAGCGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTGACACAATGCCCTGGCATTTGCTGACAGCGCCGATCTGTACAGCATCCCGCACCATGCGGCCGACGTGATc                                                                                                                               <  2:513073/150‑1 (MQ=255)                           gCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgtacagat                                                                                                                                   >  2:525825/1‑134 (MQ=255)                           gCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgtacagat                                                                                                                                   <  1:525825/141‑8 (MQ=255)                           gCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTggg                                                                                                                          >  1:288434/1‑150 (MQ=255)                           gCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTggg                                                                                                                          <  1:715569/150‑1 (MQ=255)                                     cAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGa                                                                                                                                      <  2:498193/128‑1 (MQ=255)                                     cAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGa                                                                                                                                      >  1:498193/1‑128 (MQ=255)                                       agagCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGc                                                                                                             <  1:1011403/150‑1 (MQ=255)                                               tCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCt                                                                                                     <  1:956072/150‑1 (MQ=255)                                                 aaCGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGGTCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGa                                                                                                   <  1:553832/150‑1 (MQ=255)                                                    ggCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCa‑aaaaa                                                                                                                 >  2:912861/1‑133 (MQ=255)                                                    ggCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCa‑aaaaa                                                                                                                 <  1:912861/133‑1 (MQ=255)                                                      cTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTg                                                                                              <  2:774148/150‑1 (MQ=255)                                                        aTCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGc                                                                                            <  1:1140264/150‑1 (MQ=255)                                                        aTCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTAAAAAAAAAAAA‑AAAAACCGCATGAAGCTGATATTGGc                                                                                            <  2:1027282/150‑1 (MQ=255)                                                           tGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCa                                                                                                                       <  2:776165/121‑1 (MQ=255)                                                           tGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCa                                                                                                                       >  1:776165/1‑121 (MQ=255)                                                             cAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAAc                                                                                       <  1:352750/150‑1 (MQ=255)                                                                               ctgCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTggg                                                                     >  2:1141966/3‑150 (MQ=255)                                                                               ccGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTggg                                                                     >  2:907642/1‑150 (MQ=255)                                                                                       gcgTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGa                                                                                                                                 <  2:490933/83‑1 (MQ=255)                                                                                       gcgTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGa                                                                                                                                 >  1:490933/1‑83 (MQ=255)                                                                                       gcgTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTgg                                                                                                                           >  1:398621/1‑89 (MQ=255)                                                                                       gcgTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTgg                                                                                                                           <  2:398621/89‑1 (MQ=255)                                                                                        cgTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTg                                                            <  2:395393/150‑1 (MQ=255)                                                                                           gCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATg                                                         >  2:290969/1‑150 (MQ=255)                                                                                             cAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGcc                                                                   >  1:579279/1‑138 (MQ=255)                                                                                             cAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGcc                                                                   <  2:579279/138‑1 (MQ=255)                                                                                                      gcgcAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGc                                                                                            <  2:965992/104‑1 (MQ=255)                                                                                                      gcgcAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGc                                                                                            >  1:965992/1‑104 (MQ=255)                                                                                                      gcgcAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAAcc                                                                                      >  2:904595/1‑110 (MQ=255)                                                                                                      gcgcAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAAcc                                                                                      <  1:904595/110‑1 (MQ=255)                                                                                                      gcgcAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGcc                                              >  2:426074/1‑150 (MQ=255)                                                                                                                   gCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCACAAAGACTGC‑GCAAGCTTATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGc                                 >  2:347210/1‑150 (MQ=255)                                                                                                                   gCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGcc                                                                                           <  2:471363/92‑1 (MQ=255)                                                                                                                   gCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGcc                                                                                           >  1:471363/1‑92 (MQ=255)                                                                                                                   gCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACTCCGTGGAGGCGCTGc                                 >  1:467985/1‑150 (MQ=255)                                                                                                                     aTTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCtt                                                                                   >  1:1007981/1‑98 (MQ=255)                                                                                                                     aTTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCtt                                                                                   <  2:1007981/98‑1 (MQ=255)                                                                                                                          cTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGa                          <  1:536309/150‑1 (MQ=255)                                                                                                                          cTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGa                          <  2:693200/150‑1 (MQ=255)                                                                                                                           tGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg                         <  1:1141966/150‑1 (MQ=255)                                                                                                                             aTAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGt                       <  1:426074/150‑1 (MQ=255)                                                                                                                                gcgcCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTa                    >  1:1145203/1‑150 (MQ=255)                                                                                                                                  gcCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATc                  >  1:506438/1‑150 (MQ=255)                                                                                                                                    cGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATcgc                >  2:323993/1‑150 (MQ=255)                                                                                                                                       cgtcggcAGCATGCTGCGGCATGCTGTCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg                         <  2:1014336/132‑1 (MQ=255)                                                                                                                                       cgtcggcAGCATGCTGCGGCATGCTGTCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg                         >  1:1014336/7‑138 (MQ=255)                                                                                                                                           tACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGcc                                                       <  2:1193616/104‑1 (MQ=255)                                                                                                                                           tACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGcc                                                       >  1:1193616/1‑104 (MQ=255)                                                                                                                                                 aTGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGgc   >  2:1086540/1‑150 (MQ=255)                                                                                                                                                  tGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGgcg  >  1:946066/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCCGACAGCGGCCTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCGCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCA‑AAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAAGGGCG  >  NC_002947/417227‑417521

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