Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,913,196 T→C K152K (AAA→AAG relA ← (p)ppGpp synthetase I/GTP pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,913,1960TC95.1% 137.6 / ‑0.8 42K152K (AAA→AAGrelA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (20/19);  total (21/21)
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.59e-01

AGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGTTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAAC  >  NC_000913/2913048‑2913342
                                                                                                                                                    |                                                                                                                                                  
aGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGctt                                                                                                                                                  >  2:552667/1‑151 (MQ=255)
       ttCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACt                                                                                                                                           >  1:261130/1‑151 (MQ=255)
                 cATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGaa                                                                                                                                 <  1:228733/151‑1 (MQ=255)
                    ttCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGaaaa                                                                                                                               <  1:218519/150‑1 (MQ=255)
                      cAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAatca                                                                                                                            <  2:70336/151‑1 (MQ=255)
                          tGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCg                                                                                                                         >  2:670516/1‑150 (MQ=255)
                                aTTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCg                                                                                                                  >  2:573908/1‑151 (MQ=255)
                                aTTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCg                                                                                                                  <  2:277843/151‑1 (MQ=255)
                                   ccGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCa                                                                                                               <  1:191579/151‑1 (MQ=255)
                                          ggTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACAt                                                                                                        >  1:501762/1‑151 (MQ=255)
                                            gtAGCCGTCGGTGCGTAGGTGTTGGTACTCTCTTTGGCCTCCAGTACACGTTCATCTTCCGGCTCATCTTGTCTGTCGCGCAGATGAGCAATACTCTCCGCCAGGTTCATGACTACGCAGCGAAAATCAGCGACCATCGCCAATAACATc                                                                                                       <  1:671659/149‑1 (MQ=255)
                                                 caccagTGCGTAAATGTTGGTACGCTCTTTTGCCGCCAGTACACGTTCATCCTCCGGCGCACCTTCTCCTTCGCCCAGATTCGCAACGCGCCCCACAAGCTTGATGCCTACGCAGCGAAAATCGTCGCCCATCGCCAATAACATCCGGCg                                                                                                  <  1:395864/145‑1 (MQ=255)
                                                   gCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGTTTGATGAc                                                                                                                                            >  2:608393/1‑106 (MQ=255)
                                                   gCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGTTTGATGAc                                                                                                                                            <  1:608364/106‑1 (MQ=255)
                                                        gCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTAt                                                                                          <  1:413674/151‑1 (MQ=255)
                                                            aGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGAc                                                                                      >  2:453520/1‑151 (MQ=255)
                                                               tGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTg                                                                                   >  2:502307/1‑151 (MQ=255)
                                                                 ttGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTg                                                                                   >  2:176111/1‑149 (MQ=255)
                                                                   ggTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTc                                                                                <  2:71031/150‑1 (MQ=255)
                                                                      acacTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCAGCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCggag                                                                            <  1:212853/151‑1 (MQ=255)
                                                                             tttGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGaaa                                                                        <  2:163927/148‑1 (MQ=255)
                                                                              ttGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGaaa                                                                        <  1:190052/147‑1 (MQ=255)
                                                                                gccgccAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAAt                                                                  <  2:238296/151‑1 (MQ=255)
                                                                                          cacGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGt                                                        >  2:179370/1‑151 (MQ=255)
                                                                                            cGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCg                                                      >  2:368452/1‑151 (MQ=255)
                                                                                            cGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCg                                                      >  1:356458/1‑151 (MQ=255)
                                                                                                    ttCCGGCGCACCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGAACTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGc                                              >  1:215759/1‑151 (MQ=255)
                                                                                                     tCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCa                                                <  1:670481/148‑1 (MQ=255)
                                                                                                        ggCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGc                                          >  2:658848/1‑151 (MQ=255)
                                                                                                            cATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGAt                                      <  2:319610/151‑1 (MQ=255)
                                                                                                             aTCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGAt                                      <  2:112451/150‑1 (MQ=255)
                                                                                                               cTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATcgc                                   >  2:469418/1‑151 (MQ=255)
                                                                                                               cTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATcgc                                   >  2:50132/1‑151 (MQ=255)
                                                                                                                ttttACTTCGCGCAGACGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATcgcc                                  >  2:456915/1‑151 (MQ=255)
                                                                                                                 tttACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATcgcc                                  <  1:170601/150‑1 (MQ=255)
                                                                                                                 tttACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATc                                     <  1:362714/147‑1 (MQ=255)
                                                                                                                          cgcAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATcac                        >  2:14744/1‑151 (MQ=255)
                                                                                                                                tGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCacgcacgc                  >  2:387308/1‑151 (MQ=255)
                                                                                                                                   gCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCg                <  2:332296/150‑1 (MQ=255)
                                                                                                                                       tACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAAt           >  2:431687/1‑151 (MQ=255)
                                                                                                                                         cGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGCCTTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATa          <  1:544020/150‑1 (MQ=255)
                                                                                                                                               gCCAGCTTGATGACTACGCAGCGACAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGGGTGCGTCGCTTTCAGCTGGCGGATCGGCGCCATAGCACGCAGGCCGTGAATAAGGtt     >  1:176677/1‑149 (MQ=255)
                                                                                                                                                ccAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAAc  >  2:433957/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGTTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAAC  >  NC_000913/2913048‑2913342

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: