Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,951,024 G→A 7.1% F663F (TTC→TTT ppk ← polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,951,0240GA 7.1% 156.6 / 3.0 57F663F (TTC→TTTppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (26/26);  new base A (2/2);  total (28/29)
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.39e-01

CTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCA  >  NC_002947/5950877‑5951157
                                                                                                                                                   |                                                                                                                                     
cTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGaa                                                                                                                                     >  2:1270774/1‑150 (MQ=255)
 tGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAg                                                                                                                                    <  2:335640/150‑1 (MQ=255)
  gAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGc                                                                                                                                   >  1:255358/1‑150 (MQ=255)
   aaCGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCa                                                                                                                                  >  1:1078671/1‑150 (MQ=255)
   aaCGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCa                                                                                                                                  <  2:1250030/150‑1 (MQ=255)
              ggTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGacac                                                                                                                       >  2:770863/1‑150 (MQ=255)
              ggTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGacac                                                                                                                       <  2:936399/150‑1 (MQ=255)
               gTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACg                                                                                                                      <  1:152906/150‑1 (MQ=255)
                    ggTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGt                                                                                                                 <  1:770863/150‑1 (MQ=255)
                          cGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCCGAGCTACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGt                                                                                                           <  2:816220/150‑1 (MQ=255)
                           ggTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGtt                                                                                                          <  1:613458/150‑1 (MQ=255)
                           ggTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGtt                                                                                                          >  1:564925/1‑150 (MQ=255)
                               ggggTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCTTTACCACGGTAAGCAGGTCTCGACACGCTTGTCGAGGTTgcgc                                                                                                      <  1:1210471/150‑1 (MQ=255)
                                     ctgcGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtccatc                                                                                                <  1:791142/150‑1 (MQ=255)
                                             gTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCGTCCACCGGAAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGgc                                                                                        >  2:1168732/1‑150 (MQ=255)
                                             gTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCGTCCACCGGAAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGgc                                                                                        >  2:1168785/1‑150 (MQ=255)
                                                gcgcCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCGTCCACCGGAAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCt                                                                                     <  1:1168785/150‑1 (MQ=255)
                                                gcgcCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCGTCCACCGGAAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCt                                                                                     <  1:1168732/150‑1 (MQ=255)
                                                   cccGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGtt                                                                                                          <  1:355921/126‑1 (MQ=255)
                                                   cccGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGtt                                                                                                          >  2:355921/1‑126 (MQ=255)
                                                       tCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGCTCCACAGATCGGAAGGGGAAGCATGTCTCGACACGCATGTCGCCGGACGGCTCTTTCCAGTCGGCGCTGGACAgg                                                                              >  1:665385/1‑150 (MQ=255)
                                                        cTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGt                                                                             <  1:1177429/150‑1 (MQ=255)
                                                          ggCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAg                                                                           <  1:1270774/150‑1 (MQ=255)
                                                                 aGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGACGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGct                                                                    <  2:922599/150‑1 (MQ=255)
                                                                 aGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGct                                                                    <  2:624407/150‑1 (MQ=255)
                                                                   gTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGc                                                                                                                     >  2:303552/1‑99 (MQ=255)
                                                                   gTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGc                                                                                                                     <  1:303552/99‑1 (MQ=255)
                                                                                tATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGa                                                     <  2:564925/150‑1 (MQ=255)
                                                                                tATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATATGCTCTTCGCCGCCATTGa                                                     >  2:733923/1‑150 (MQ=255)
                                                                                 aTTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAg                                                    <  2:1078671/150‑1 (MQ=255)
                                                                                        gTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAg                                             >  1:659655/1‑150 (MQ=255)
                                                                                         tCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTACTTGCATTCCACCAGGAAGCAGGTCTGGGCAGGCTTATCGACAGAAGTCGCCAACCGGTCGCCGCGGTACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGa                                            >  2:786121/1‑150 (MQ=255)
                                                                                            ggTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGt                                               <  2:1180625/144‑1 (MQ=255)
                                                                                            ggTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGt                                               >  1:1180625/1‑144 (MQ=255)
                                                                                                      cAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCCCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGtctg                               >  1:67412/1‑150 (MQ=255)
                                                                                                       aGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtcc                                                                                                   >  2:1152669/1‑81 (MQ=255)
                                                                                                       aGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtcc                                                                                                   <  1:1152669/81‑1 (MQ=255)
                                                                                                        gCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtcc                                                                                                   <  2:1052902/80‑1 (MQ=255)
                                                                                                        gCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtcc                                                                                                   >  1:1052902/1‑80 (MQ=255)
                                                                                                        gCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCt                             >  1:1239969/1‑150 (MQ=255)
                                                                                                           ccTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGa                          <  1:1189357/150‑1 (MQ=255)
                                                                                                                ttcACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGATTGCGCTCGATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGTCGCCATTGAGGACGTAGAACACCCGCGTGTGCTCGAGGAAg                     >  1:998929/1‑150 (MQ=255)
                                                                                                                  cACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCg                   >  2:606225/1‑150 (MQ=255)
                                                                                                                     ccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGAc                                                                                 >  2:995122/1‑85 (MQ=255)
                                                                                                                     ccGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGAc                                                                                 <  1:995122/85‑1 (MQ=255)
                                                                                                                       gcagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAaca                                         >  2:72513/1‑123 (MQ=255)
                                                                                                                       gcagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAaca                                         <  1:72513/123‑1 (MQ=255)
                                                                                                                       gcagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCaa              <  2:659655/150‑1 (MQ=255)
                                                                                                                        cagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGcc                <  2:277298/147‑1 (MQ=255)
                                                                                                                        cagcaACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGcc                >  1:277298/1‑147 (MQ=255)
                                                                                                                            aaCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTTGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATc         >  1:380272/1‑150 (MQ=255)
                                                                                                                            aaCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATc         >  2:1000241/1‑150 (MQ=255)
                                                                                                                            aaCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGATGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATc         <  2:255358/150‑1 (MQ=255)
                                                                                                                             aCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCg        >  1:810533/1‑150 (MQ=255)
                                                                                                                                   tCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAACTGCCTTCTACGCCGCCATTGAGGAAGTAGAACACCAGCGTGTGCTCGATGAAGCGGCCAATGATCGAGCGCa  >  1:988189/1‑150 (MQ=255)
                                                                                                                                       gCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACAcc                                       >  2:1270486/1‑109 (MQ=255)
                                                                                                                                       gCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACAcc                                       <  1:1270486/109‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                     
CTGAACGCTGCGCGGGTTCTGGTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCAACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCA  >  NC_002947/5950877‑5951157

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