Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,049,348 G→A 10.3% V52M (GTG→ATG) ‡ spoT → bifunctional (p)ppGpp synthase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,049,3480GA10.3% 151.9 / 8.1 58V52M (GTG→ATG) ‡spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base G (23/29);  new base A (3/3);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

CTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGT  >  NC_002947/6049212‑6049493
                                                                                                                                        |                                                                                                                                                 
cTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgggtg                                                                                                                                      <  1:258538/150‑6 (MQ=255)
 ttGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTc                                                                                                                                     <  2:549851/150‑1 (MQ=255)
 ttGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTc                                                                                                                                     >  2:1117951/1‑150 (MQ=255)
       gAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGAc                                                                                                                               >  1:1242011/1‑150 (MQ=255)
        aaCGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACa                                                                                                                              <  2:1242011/150‑1 (MQ=255)
         aCGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACAt                                                                                                                             >  1:321856/1‑150 (MQ=255)
           ggCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGc                                                                                                                           >  2:1406604/1‑150 (MQ=255)
              tttCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGcaca                                                                                                                        >  1:616826/1‑150 (MQ=255)
              tttCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGCCGTGCCTATTTCTACCCCGAACAGGCCCACAACGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGcaca                                                                                                                        >  1:320384/1‑150 (MQ=255)
                   aCCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCCCCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGAc                                                                                                                   >  1:683392/1‑150 (MQ=255)
                                        gTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGcc                                                                                              >  2:807088/1‑150 (MQ=255)
                                           aaCCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCCGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCAtg                                                                                           >  2:711340/1‑150 (MQ=255)
                                           aaCCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgggtg                                                                                                                                      >  2:101855/1‑102 (MQ=255)
                                           aaCCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgggtg                                                                                                                                      <  1:101855/107‑6 (MQ=255)
                                                  ttCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACg                                                                                    >  2:1019531/1‑150 (MQ=255)
                                                          gCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc                                                                            <  1:1019531/150‑1 (MQ=255)
                                                           ccTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCAGATGGCGGCCATGCTGCACGATGTGATCg                                                                           >  1:767061/1‑150 (MQ=255)
                                                               tttCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc                                                                            >  1:737759/1‑145 (MQ=255)
                                                               tttCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc                                                                            <  2:737759/145‑1 (MQ=255)
                                                                              ggCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCaa                                                        <  2:683392/150‑1 (MQ=255)
                                                                               gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTg                                                                                                                                                 <  1:850791/60‑1 (MQ=255)
                                                                               gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTg                                                                                                                                                 >  2:850791/1‑60 (MQ=255)
                                                                                 ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtg                                                                                 <  2:1161017/122‑1 (MQ=255)
                                                                                 ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtg                                                                                 >  1:1161017/1‑122 (MQ=255)
                                                                                 ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCg                                                                  >  2:1166054/1‑137 (MQ=255)
                                                                                 ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCg                                                                  <  1:1166054/137‑1 (MQ=255)
                                                                                 ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGa                                                     >  2:183641/1‑150 (MQ=255)
                                                                                 ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGAGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGa                                                     >  2:920427/1‑150 (MQ=255)
                                                                                   aCGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAg                                                   <  2:1174290/150‑1 (MQ=255)
                                                                                    cGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGc                                                  <  1:477735/150‑1 (MQ=255)
                                                                                     gATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGcc                                                 <  1:371568/150‑1 (MQ=255)
                                                                                         ggCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTca                                             <  1:881737/150‑1 (MQ=255)
                                                                                         ggCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTca                                             <  2:616826/150‑1 (MQ=255)
                                                                                          gCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTcag                                            >  2:1060105/1‑150 (MQ=255)
                                                                                             aCGCCGCCGGAGTGGCGAGCCCTACGGGACGCAACCGCTGGCCGTGGGCAGCATCGTGGCCGACATGCACATGGACAGTCAGAGGCGGATGGCGGCCATGCTGCACGATGTGATCGAGGTCACCGGCATCGCGAAGGAAGCCCTcagcca                                         <  2:320384/150‑1 (MQ=255)
                                                                                             aCGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTcagcca                                         <  2:899288/150‑1 (MQ=255)
                                                                                               gccgccgcAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTcagccagc                                       >  1:250688/1‑150 (MQ=255)
                                                                                                ccgccgcAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCcagca                                      <  1:1060105/150‑1 (MQ=255)
                                                                                                   ccgcAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATggcg                                                                                                 <  1:869740/88‑1 (MQ=255)
                                                                                                   ccgcAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATggcg                                                                                                 >  2:869740/1‑88 (MQ=255)
                                                                                                               gCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCa                                                                                                              <  2:434278/63‑1 (MQ=255)
                                                                                                               gCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCa                                                                                                              >  1:434278/1‑63 (MQ=255)
                                                                                                               gCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGt                       <  2:946738/150‑1 (MQ=255)
                                                                                                                cccTACGTGACCCACCCGCTGGCCGTGGCCAGCATCACCGCACCAATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCgcg                      <  2:944719/150‑3 (MQ=255)
                                                                                                                           cccccccCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGTTATTCTGCACGATGTGATCGAAGACACCGGCATCGCCAATGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGt           <  1:734604/147‑1 (MQ=255)
                                                                                                                           ccacccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGt           <  2:250688/150‑1 (MQ=255)
                                                                                                                              gatgCTGGCCACGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGa        >  1:306777/4‑150 (MQ=255)
                                                                                                                              cccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAATCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGa        <  1:524049/150‑1 (MQ=255)
                                                                                                                              cccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGa        <  1:1406604/150‑1 (MQ=255)
                                                                                                                               ccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGAt       <  2:321856/150‑1 (MQ=255)
                                                                                                                                tgCTGGCCACGGCCAGCATCCTCGCCGACATGCACATGGAc                                                                                                                   <  2:174442/40‑1 (MQ=255)
                                                                                                                                tgCTGGCCACGGCCAGCATCCTCGCCGACATGCACATGGAc                                                                                                                   >  1:174442/2‑41 (MQ=255)
                                                                                                                                 gCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtgt                                                                                <  1:1303107/75‑1 (MQ=255)
                                                                                                                                 gCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATgg     >  1:546842/1‑150 (MQ=255)
                                                                                                                                 gCTGGACGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCTGCCATGCTGCACGAtgt                                                                                >  2:1303107/1‑75 (MQ=255)
                                                                                                                                  cTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATTCTGCACTATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATggg    <  1:807088/150‑1 (MQ=255)
                                                                                                                                    ggCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGt  <  1:1117951/150‑1 (MQ=255)
                                                                                                                                    ggCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGt  <  1:616749/150‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                 
CTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGT  >  NC_002947/6049212‑6049493

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: