Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,913,421 G→A 17.3% R77W (CGG→TGG)  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,4210GA17.3% 83.1 / 8.0 46R77W (CGG→TGG) relA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base G (12/25);  new base A (3/5);  total (15/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTCACACGACTTCT  >  NC_000913/2913277‑2913564
                                                                                                                                                |                                                                                                                                                
cAGTGTGCGTCGCTTTCAGCTGGTGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGc                                                                                                                                              <  2:376395/149‑1 (MQ=255)
cAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCGGAGGGAAAAGCAGCGCCGCCCGCAGc                                                                                                                                              >  1:641806/1‑149 (MQ=255)
  gtgtGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCgt                                                                                                                                            <  1:693433/149‑1 (MQ=255)
  gtgtGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCgt                                                                                                                                            <  1:74508/149‑1 (MQ=255)
    gtgCGTCGCTTTCAGCTGGCGGATCGCCGCCATAACACGAACGCCCTGAATAAGGGGAACGACCGACTCACCGACGCTCTCACGCAGCACATCTCCGCTGACTACGTTGGCATCCGCAAGAGGGAAAACCAGCGCCCCCCACAGCgctt                                                                                                                                          >  2:28292/1‑146 (MQ=255)
        gtgcgtttCAGCTGGCTGATCCCCGCCATATACGGCACGCCGTGAATAAGGTTAACTACCGAGTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTGGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAAt                                                                                                                                      <  2:20590/144‑1 (MQ=255)
        gTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCACAGCGTGTCAAt                                                                                                                                      >  1:82095/1‑149 (MQ=255)
             tttCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCa                                                                                                                                 >  2:366596/1‑149 (MQ=255)
             tttCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACCACGCTGGCATCCGCCATAGGGAAACGCAGCGCCGCCCGCAGCGTGTCAATGTCCa                                                                                                                                 <  1:548632/149‑1 (MQ=255)
                cAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATAc                                                                                                                              <  1:578563/149‑1 (MQ=255)
                 aGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGCCCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAAGGTCCATACt                                                                                                                             >  1:93721/1‑149 (MQ=255)
                        cgATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTc                                                                                                                      <  1:124854/148‑1 (MQ=255)
                              ccgccATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAgg                                                                                                                >  1:575525/1‑149 (MQ=255)
                                  cATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGAtct                                                                                                            <  1:302827/149‑1 (MQ=255)
                                              gcCGTTAATGAGGTTAACTCCCGACTTCCCGACGCCCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCTCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCaa                                                                                                <  2:30064/149‑1 (MQ=255)
                                              gcCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCCCCACCATCTCaa                                                                                                >  1:135387/1‑149 (MQ=255)
                                                      tAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGCGGATCTCCACCATCTCAACACCACGc                                                                                        >  1:2312/1‑149 (MQ=255)
                                                         ggTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCCCAACACCCCGCcac                                                                                     >  2:243466/1‑149 (MQ=255)
                                                         ggTGAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCcac                                                                                     <  1:236384/149‑1 (MQ=255)
                                                                gaccgacTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACa                                                                             >  1:395670/1‑149 (MQ=255)
                                                                  ccgacTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGa                                                                           <  1:127735/149‑1 (MQ=255)
                                                                     acTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCACTGCCCATCTAGAAGTTCGAGAGGATCTCCACCACCCCAACACCACGCCACAATAACAGACTg                                                                        >  1:321612/1‑149 (MQ=255)
                                                                       ttACCGACCCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGc                                                                      <  2:651731/149‑1 (MQ=255)
                                                                        tACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGGCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCa                                                                     <  1:505270/149‑1 (MQ=255)
                                                                              cGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGa                                                               <  1:614348/149‑1 (MQ=255)
                                                                              cGCTCTCACGCAGCACATCTCCGCTGACTACGTTGGCATCCGGCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGa                                                               <  1:682539/149‑1 (MQ=255)
                                                                               gCTCTGACGAACCACCTCTTCGCTGACTCCGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCACAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGTATCCGGAt                                                              <  2:93043/149‑1 (MQ=255)
                                                                                     cccccgcaCATCTTCGCGGACTACGTTGCCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCt                                                        <  2:341387/144‑1 (MQ=255)
                                                                                         agcaCATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGt                                                    <  1:393443/149‑1 (MQ=255)
                                                                                            acaTCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTAAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTg                                                 <  2:261478/149‑1 (MQ=255)
                                                                                            acaTCGTCGCTGACTACGTGGGCATCCGCCAGAGGGAAAAGCCGCGCCGCCCGCAGCGTGGCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTAAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTg                                                 <  2:261377/149‑1 (MQ=255)
                                                                                                cTTCGGTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGc                                             <  2:272568/149‑1 (MQ=255)
                                                                                                cTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGc                                             >  1:168383/1‑149 (MQ=255)
                                                                                                cTTCGCTGACTACGTTGGCATCCGCAAGAGGGAAAAGCAGCGCCGCCCACAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGc                                             <  2:640654/149‑1 (MQ=255)
                                                                                                 ttCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCa                                            <  2:395670/149‑1 (MQ=255)
                                                                                                        acgaCGTTGGCATCCGCCAGAGGGAAAAGCCGCGCCGCCCGGAGCGTGTCATTGTCCATACTCAATGGCGAGAGAATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCTGACAATa                                     <  2:310636/146‑1 (MQ=255)
                                                                                                             gTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCGTGCGTCTGTTGCAGACAATACGccc                                <  1:366596/149‑1 (MQ=255)
                                                                                                              gtGGCATCCGCCAGAGGGAACAGCAGCGCCGCCCACAGCGTGTCAATGTACATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCa                               <  2:481152/148‑1 (MQ=255)
                                                                                                                  cATCCGCCAGAGGGAAAAGCAGCGCCGCCCACAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGGtt                           >  2:538171/1‑149 (MQ=255)
                                                                                                                    tCCGCCAGAGGGAAAAGCAGCGCCGCCCGAAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCCTCCGGATGCCCCTGCGCCTGTTGCAGACACAACTCCCAgcgttc                         >  2:333997/1‑144 (MQ=255)
                                                                                                                          agagGGAAAAGCAGCGCCGCCCACAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAGCACCACGCCACAATAACCGACTGGCATCCGGATGCCCCTGCGTCTGTTGCTGACAATACGCCCAGGTTTCGGCTaa                   <  2:344090/149‑1 (MQ=255)
                                                                                                                               gAAAAGCAGCGCCGCCCGCCGCGTGTCAATGTCCATACTTAATGTCGTGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACt              <  2:624862/149‑1 (MQ=255)
                                                                                                                               gAAAAGCAGCGCCACCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGATTTCGGCTAAGCAcc              >  2:384451/1‑148 (MQ=255)
                                                                                                                                  aaGCAGCGCCGCCCACAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTcac           <  1:60499/149‑1 (MQ=255)
                                                                                                                                  aaGCAGCCCCGCCCGCAGCGTGTCAATGTCCATACTTCATGTCGCGAGGCTCTCCACCATCTCAACACCACGCCACAATCCACCGGTGCCATCCGGATGCTCCTGCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTcac           <  1:426899/149‑1 (MQ=255)
                                                                                                                                      agcGCCCCCCGCAGCCTGTCAGTGTCCATATTTAATTTCGTGAGGATATCTGCCATCTCAACACCCCGCGACAATAGGAGTCTGGCAGCCTGATGCCCTTGCTTCTGTTGCAGACAATACGCACAGGTTTCGGCTAAGCACTCACACGa       <  2:38903/149‑1 (MQ=255)
                                                                                                                                          ccgcccacAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCCACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTCCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTCACACGACTTc   >  1:702917/1‑149 (MQ=255)
                                                                                                                                           ggaccgcAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTAAACACCACGCCACAATCACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTCACACGACttat  <  2:210594/146‑3 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAATAACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTCACACGACTTCT  >  NC_000913/2913277‑2913564

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