Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	418,298	C→G	7.0%	L26L (CTC→CTG)	PP_0345 →	heptose kinase WapP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	418,298	0	C	G	7.0%	172.8 / 4.9	57	L26L (CTC→CTG)	PP_0345	heptose kinase WapP
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base G (2/2);  total (27/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 = 4.24e-01

CTCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCAAA  >  NC_002947/418153‑418446                                                                                                                                                   |                                                                                                                                                     ctCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTAGATTG‑GAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                  <  1:68967/150‑1 (MQ=255) ctCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                  <  2:1000071/150‑1 (MQ=255)    gCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑TGTGACTGTTGGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCg                                                                                                                                               <  1:466326/150‑1 (MQ=255)      ctggctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGaa                                                                                                                                             >  1:949951/1‑150 (MQ=255)         gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACgg                                                                                                                                          >  1:960574/1‑150 (MQ=255)         gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGCGTTCGCCCTGGAGGGCGAAcac                                                                                                                                          <  2:784002/150‑3 (MQ=255)          ctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGc                                                                                                                                         <  2:812278/150‑1 (MQ=255)           tggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCt                                                                                                                                        >  2:353443/1‑150 (MQ=255)             gAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGTGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGNTga                                                                                                                                      <  2:698968/150‑1 (MQ=255)                  cAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgc                                                                                                                                 >  1:1028696/1‑150 (MQ=255)                        ccAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacccg                                                                                                                           >  2:1100381/1‑150 (MQ=255)                          aGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcgacccgcc                                                                                                                         >  2:899407/1‑149 (MQ=255)                            aaaCTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCa                                                                                                                       <  1:376643/150‑1 (MQ=255)                                tctACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                  >  2:1014464/1‑119 (MQ=255)                                tctACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                  <  1:1014464/119‑1 (MQ=255)                                 ctACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGaa                                                                                                                  <  1:929281/150‑1 (MQ=255)                                               aaCGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGa                                                                                                    <  2:949951/150‑1 (MQ=255)                                                aCGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                  >  2:1017773/1‑103 (MQ=255)                                                aCGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                  <  1:1017773/103‑1 (MQ=255)                                                  gCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCg                                                                                                 >  2:12615/1‑150 (MQ=255)                                                     aTGGAGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGa                                                                                              <  1:596755/150‑1 (MQ=255)                                                       ggggATGCACTCT‑GATGGTGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACg                                                                                            <  2:241866/150‑1 (MQ=255)                                                          gATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACNGGa                                                                                         <  2:694886/150‑1 (MQ=255)                                                             cggccgt‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCa                                                                                      >  2:949739/7‑150 (MQ=255)                                                               aCTCT‑GATGGCGAGTTGGACGCTTGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTTGAAGGCTAGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                    <  1:619411/150‑1 (MQ=255)                                                               aCTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTGGAGGGCGAAc                                                                                                                                            >  1:8682/1‑94 (MQ=255)                                                               aCTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                    >  2:168182/1‑150 (MQ=255)                                                               aCTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                    >  2:897101/1‑150 (MQ=255)                                                               aCTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                    <  1:903089/150‑1 (MQ=255)                                                               aCTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACCTCGGCAGCCTCGACGCGGTGTTCGCCCTGGAGGGCGAAc                                                                                                                                            <  2:8682/94‑1 (MQ=255)                                                                              ttGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAgcg                                                                      >  1:1010761/1‑150 (MQ=255)                                                                                 gACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTa                                                                   >  1:433949/1‑150 (MQ=255)                                                                                 gACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGAGACCAACTACTACGTCAAGCGCTa                                                                   >  2:562213/1‑150 (MQ=255)                                                                                    acTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTacac                                                                >  1:161393/1‑150 (MQ=255)                                                                                          gcCGGGGTACGAACACCTGGCGGGGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGc                                                          >  2:551354/1‑150 (MQ=255)                                                                                             ggggTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAACCGGATCGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCgg                                                       <  2:650960/150‑1 (MQ=255)                                                                                                gTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCaa                                                    <  1:12615/150‑1 (MQ=255)                                                                                                           tGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGcc                                         <  2:161393/150‑1 (MQ=255)                                                                                                               ggcggACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTa                                     <  1:13194/150‑1 (MQ=255)                                                                                                                   gACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTg                                 <  2:960574/150‑1 (MQ=255)                                                                                                                    tgtgcGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTgg                                >  2:1060650/5‑150 (MQ=255)                                                                                                                    aCTTCGGCAGCCTCGACGCGGTGTTCGTCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTgg                                <  2:178446/150‑1 (MQ=255)                                                                                                                          gCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACAt                                              >  1:283642/1‑130 (MQ=255)                                                                                                                          gCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcggc                          <  2:222655/150‑1 (MQ=255)                                                                                                                          gCACCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACAt                                              <  2:283642/130‑1 (MQ=255)                                                                                                                                       ggTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGc             <  2:937381/150‑1 (MQ=255)                                                                                                                                         tgtTCGCGCTCCAGGGCGAACGGCTGACCCTCGACCCGCTCAGTTAAATGGTACGCATCGAGCGTGACGGGATCAAGTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGTTGGGACCTGGGGCGGCCACGCATCAAGGCCg           <  1:949739/150‑1 (MQ=255)                                                                                                                                         cgtTCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCg           >  1:47183/2‑150 (MQ=255)                                                                                                                                           ttCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGaa         >  1:315328/1‑150 (MQ=255)                                                                                                                                             cGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGTGGGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGTCCACGCATCAAGGCCGAATg       <  1:562213/150‑1 (MQ=255)                                                                                                                                              gCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACAt                                              >  2:785554/1‑110 (MQ=255)                                                                                                                                              gCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACAt                                              <  1:785554/110‑1 (MQ=255)                                                                                                                                              gCCCTCCAGGGCGAACGGCTGACCCGCAACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCACCTCCTACGTCACGCGCTACCCCGGGGCCGGCCAGCACATACGCCGGTACCCTGGGCGGCCACGCATCAAGGCCGAATgg      >  1:586046/1‑150 (MQ=255)                                                                                                                                                 cTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAg               >  2:760471/1‑138 (MQ=255)                                                                                                                                                 cTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAg               <  1:760471/138‑1 (MQ=255)                                                                                                                                                 cTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCaa   >  2:1100596/1‑150 (MQ=255)                                                                                                                                                  tCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCACCTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCCTCAAGGCCGAATGCCaaa  >  2:925189/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CTCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCT‑GATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCAAA  >  NC_002947/418153‑418446

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