Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,950,987 T→G 7.9% K657T (AAG→ACG)  ppk ← polyphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,950,9870TG 7.9% 143.0 / 5.3 51K657T (AAG→ACG) ppkpolyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base G (2/2);  total (25/26)
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.26e-01

GTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAATGCCCG  >  NC_002947/5950842‑5951134
                                                                                                                                                  |                                                                                                                                                   
gTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑GCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGt                                                                                                                                                  >  8:76727/1‑149 (MQ=255)
gTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCCGGAAGCAGGTCTCGACACGCTTg                                                                                                                                                   >  2:244536/1‑148 (MQ=255)
    ccGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGa                                                                                                                                               <  5:287521/148‑1 (MQ=255)
    ccGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTACAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGGTCCTTCTTCCCCCGCAGCA‑ACAGTTTCTTGCCTTCCCCCGGGAAGCAGGTCTCGACACGCTTGTCGAg                                                                                                                                              >  6:242506/1‑149 (MQ=255)
     cGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAgg                                                                                                                                             <  5:139199/149‑1 (MQ=255)
     cGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGAga                                                                                                                                             <  1:77006/149‑2 (MQ=255)
              tgctgcGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑GCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtcc                                                                                                                                    <  7:76727/149‑1 (MQ=255)
                  gcgcACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtccatc                                                                                                                                 >  8:86626/1‑148 (MQ=255)
                   cgcACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCtccatc                                                                                                                                 <  4:71736/147‑1 (MQ=255)
                          aGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGATAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGgcg                                                                                                                        >  4:695/1‑149 (MQ=255)
                           gcgcCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGgcgc                                                                                                                       <  4:117315/149‑1 (MQ=255)
                            cgcCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGTGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGgcgc                                                                                                                       <  6:156892/148‑1 (MQ=255)
                            cgcCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCt                                                                                                                      >  6:71722/1‑149 (MQ=255)
                                     ggCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTGTTCACCCGCAGGA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGt                                                                                                              <  6:30509/148‑1 (MQ=255)
                                      gCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGt                                                                                                              <  5:71722/147‑1 (MQ=255)
                                      gCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGt                                                                                                              <  1:244536/147‑1 (MQ=255)
                                        tGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGa                                                                                                           >  7:91943/1‑148 (MQ=255)
                                        tGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGa                                                                                                           >  1:2770/1‑148 (MQ=255)
                                        tGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGAgacct                                                                                                                                          >  3:327055/1‑113 (MQ=255)
                                        tGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACAAGCGTGTCGAgacct                                                                                                                                          <  4:327055/117‑5 (MQ=255)
                                              gTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTAGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGctc                                                                                                    <  3:73557/149‑1 (MQ=255)
                                              gTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTATTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGc                                                                                                      >  5:279042/1‑147 (MQ=255)
                                              gTCCAGGAGTGGGTATTGTCTGTCAGTTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTTCCTTCCACCGGGAAGCAGGTCTCGACACGATTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGc                                                                                                      <  6:279042/147‑1 (MQ=255)
                                                ccAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCtt                                                                                                  >  1:247621/1‑149 (MQ=255)
                                                     cGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTcgccg                                                                                             <  7:86626/149‑1 (MQ=255)
                                                      gTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTcgccgc                                                                                            >  5:230855/1‑149 (MQ=255)
                                                           tATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTg                                                                                       <  2:2770/149‑1 (MQ=255)
                                                               gTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAgg                                                                                    <  2:247621/148‑1 (MQ=255)
                                                                          aGCCTTCCAGCTCCTTCTTCACCCGCAGCATACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGaa                                                                            >  1:31333/1‑146 (MQ=255)
                                                                          aGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGaa                                                                            <  2:31333/145‑1 (MQ=255)
                                                                              ttCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCgg                                                                                                                          <  7:129882/95‑1 (MQ=255)
                                                                              ttCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCgg                                                                                                                          >  8:129882/1‑95 (MQ=255)
                                                                                     tCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTc                                                             >  1:227830/1‑149 (MQ=255)
                                                                                      ccTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTc                                                             >  8:150590/1‑148 (MQ=255)
                                                                                          cttcACCAGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAg                                                          <  7:80496/147‑1 (MQ=255)
                                                                                                  gcagca‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGcc                                                 >  8:297251/1‑148 (MQ=255)
                                                                                                       a‑aCAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGAt                                           >  5:191690/1‑149 (MQ=255)
                                                                                                          cAGTTGCTTGCCTTCCACCGCTAAGCATGTCTCGACACGCTTGTCGAGGTTGCGCTACATCCCGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCg                                         <  8:134797/149‑1 (MQ=255)
                                                                                                           aGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATc                                          >  1:56090/1‑147 (MQ=255)
                                                                                                           aGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGa                                        <  6:230855/149‑1 (MQ=255)
                                                                                                            gTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAg                                       <  8:91943/149‑1 (MQ=255)
                                                                                                                      ttCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGa                             >  8:224145/1‑149 (MQ=255)
                                                                                                                        ccACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAt                            <  3:695/148‑1 (MQ=255)
                                                                                                                         cACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtgt                          >  8:64101/1‑149 (MQ=255)
                                                                                                                            cGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtgttgt                       <  6:191690/149‑1 (MQ=255)
                                                                                                                            cGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGAtgttgt                       <  8:34098/149‑1 (MQ=255)
                                                                                                                                     aacctcGACAAGCGTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAAccc              >  6:296127/4‑149 (MQ=255)
                                                                                                                                         tcGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAAcccc             <  2:107556/146‑1 (MQ=255)
                                                                                                                                         tcGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAAcccc             >  1:107556/1‑146 (MQ=255)
                                                                                                                                                gCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGACCACCCGCGTGTGCTCGAGGAAGCGGCCAATGAACGCGCGCCCCTGGATGATGTGCGAAACCCCCGGAATGccc   >  4:62461/1‑149 (MQ=255)
                                                                                                                                                 cTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTTGCCGCCATTGAGGAAGTAGAACCCCCGCGTGTGCTCGAGGAAGCGGCCAATGAGCGAGCGCACCTGGATGTTGTGCGACACCCCCGGAATGCCCg  >  3:65464/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTTGCCGGTCGGGGTGCTGCGCACGTAGCGCCCGTCTGGCTGCAAGGTCCAGGCGTGGGTATTGTCGGTCAGGTAGCCTTCCAGCTCCTTCTTCACCCGCAGCA‑ACAGTTTCTTGCCTTCCACCGGGAAGCAGGTCTCGACACGCTTGTCGAGGTTGCGCTCCATCCAGTCGGCGCTGGACAGGTAGATCTGCTCTTCGCCGCCATTGAGGAAGTAGAACACCCGCGTGTGCTCGAGGAAGCGGCCAATGATCGAGCGCACCTGGATGTTGTGCGAAACCCCCGGAATGCCCG  >  NC_002947/5950842‑5951134

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