Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	417,378	A→G	9.5%	Q362R (CAG→CGG)	waaG →	lipopolysaccharide glucosyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	417,378	0	A	G	9.5%	109.7 / 3.2	42	Q362R (CAG→CGG)	waaG	lipopolysaccharide glucosyltransferase I
Reads supporting (aligned to +/- strand):  ref base A (20/18);  new base G (2/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

CTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAA  >  NC_002947/417239‑417516                                                                                                                                             |                                                                                                                                           cTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGcc                                                                                                                                  <  1:25163/150‑1 (MQ=255)  tGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCTGCGGCAtgctgtac                                                                                                                                 <  2:380477/150‑4 (MQ=255)                    tCGAACAAGAGCAGCTCAACGGCTCGCGGCAGAGCATGCTCGAAGACCTGCAGGCGCGTGCCAGGTGGTCGCG‑AAATGGCCTGGCATTAGCTGGTAGCGCCGATCTGTACAGCACGCGGCAGCATGCTGCCGACGTGATCCTGGGACaaa                                                                                                               <  2:62033/150‑1 (MQ=255)                         cAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAAcc                                                                                                          <  2:323305/150‑1 (MQ=255)                                   tCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCt                                                                                                >  1:140010/1‑150 (MQ=255)                                          cTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTg                                                                                         >  2:741110/1‑150 (MQ=255)                                                gCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGaa                                                                                   >  1:354473/1‑150 (MQ=255)                                                 cAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGCT‑GCGTCAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAAc                                                                                  <  1:7250/150‑1 (MQ=255)                                                          cTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAAc                                                                         >  2:195799/1‑150 (MQ=255)                                                           tCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACg                                                                        >  1:221031/1‑150 (MQ=255)                                                                   ccGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTggg                                                                >  1:236822/1‑150 (MQ=255)                                                                   ccGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTggg                                                                >  1:402505/1‑150 (MQ=255)                                                                   ccGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTggg                                                                >  2:657054/1‑150 (MQ=255)                                                                   ccGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTggg                                                                >  1:560849/1‑150 (MQ=255)                                                                    cGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGc                                                               >  1:162442/1‑150 (MQ=255)                                                                        ggCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAg                                                                                                  <  1:297750/111‑1 (MQ=255)                                                                        ggCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAg                                                                                                  >  2:297750/1‑111 (MQ=255)                                                                         gcgcGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGc                                                          <  1:741110/150‑1 (MQ=255)                                                                              tGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGAt                                                     <  2:345697/150‑1 (MQ=255)                                                                                       gTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACg                                                                                                                              >  2:506472/1‑68 (MQ=255)                                                                                       gTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACg                                                                                                                              <  1:506472/68‑1 (MQ=255)                                                                                            g‑cAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGt                                       >  2:97260/1‑150 (MQ=255)                                                                                                      tGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCt                              <  1:28657/150‑1 (MQ=255)                                                                                                       ggCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTg                             >  1:517204/1‑150 (MQ=255)                                                                                                           tttGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAg                         <  1:510347/150‑1 (MQ=255)                                                                                                                tGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg                    <  2:402505/150‑1 (MQ=255)                                                                                                                tGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCGGTGGNCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGac                    <  1:195799/150‑2 (MQ=255)                                                                                                                      cgcCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAgg                   <  2:441071/145‑1 (MQ=255)                                                                                                                      cgcCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAgg                   >  1:441071/1‑145 (MQ=255)                                                                                                                         cGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTc                                               <  2:482868/114‑1 (MQ=255)                                                                                                                         cGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTc                                               >  1:482868/1‑114 (MQ=255)                                                                                                                           aTCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg                    >  2:73857/1‑139 (MQ=255)                                                                                                                           aTCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg                    <  1:73857/139‑1 (MQ=255)                                                                                                                              tcggcAGCATGCTGCGGCATGCTGTCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg                    >  2:712187/5‑136 (MQ=255)                                                                                                                              tcggcAGCATGCTGCGGCATGCTGTCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAg                    <  1:712187/132‑1 (MQ=255)                                                                                                                              tGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGGGGCGCTGCAAGGCGAGGTCTATCGCGAGCt      <  2:236822/150‑1 (MQ=255)                                                                                                                              tGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTa               >  2:107625/1‑141 (MQ=255)                                                                                                                              tGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTa               <  1:107625/141‑1 (MQ=255)                                                                                                                               gTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTggg                                                                                                                    >  2:111382/1‑39 (MQ=255)                                                                                                                               gTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTggg                                                                                                                    <  1:111382/39‑1 (MQ=255)                                                                                                                                tACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTgg    >  2:22161/1‑150 (MQ=255)                                                                                                                                  cAGCATGCTGCGGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGaa  >  1:507104/1‑150 (MQ=255)                                                                                                                                             |                                                                                                                                           CTGGTGCTGGACGAGCCGTTCGAACAAGAGCAGCTCAACGGCTATCTGCAGCGCATGCTCGAAGACCCGCAGGCGCGTGCCAGCTGGTCGCG‑CAATGGCCTGGCATTTGCTGATAGCGCCGATCTGTACAGCATGCCGCAGCATGCTGCCGACGTGATCCTGGGGCAAAAAACCGCATGAAGCTGATATTGGCCGAACCCTTCAAACGCCTGTGGGCCGGGCGTGATGCCTTCGACGCCGTGGAGGCGCTGCAAGGCGAGGTCTATCGCGAGCTGGAA  >  NC_002947/417239‑417516

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