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

CAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTCACACGACTTCT  >  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)
                                                                gaccgacTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACa                                                                             >  1:395670/1‑149 (MQ=255)
                                                                  ccgacTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGa                                                                           <  1:127735/149‑1 (MQ=255)
                                                                     acTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCACTGCCCATCTAGAAGTTCGAGAGGATCTCCACCACCCCAACACCACGCCACAAT‑AACAGACTg                                                                        >  1:321612/1‑149 (MQ=255)
                                                                       ttACCGACCCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGc                                                                      <  2:651731/149‑1 (MQ=255)
                                                                        tACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGGCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCa                                                                     <  1:505270/149‑1 (MQ=255)
                                                                              cGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGa                                                               <  1:614348/149‑1 (MQ=255)
                                                                              cGCTCTCACGCAGCACATCTCCGCTGACTACGTTGGCATCCGGCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGa                                                               <  1:682539/149‑1 (MQ=255)
                                                                               gCTCTGACGAACCACCTCTTCGCTGACTCCGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCACAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGTATCCGGAt                                                              <  2:93043/149‑1 (MQ=255)
                                                                                     cccccgcaCATCTTCGCGGACTACGTTGCCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCt                                                        <  2:341387/144‑1 (MQ=255)
                                                                                         agcaCATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGt                                                    <  1:393443/149‑1 (MQ=255)
                                                                                            acaTCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTAAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTg                                                 <  2:261478/149‑1 (MQ=255)
                                                                                            acaTCGTCGCTGACTACGTGGGCATCCGCCAGAGGGAAAAGCCGCGCCGCCCGCAGCGTGGCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTAAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTg                                                 <  2:261377/149‑1 (MQ=255)
                                                                                                cTTCGGTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGc                                             <  2:272568/149‑1 (MQ=255)
                                                                                                cTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGc                                             >  1:168383/1‑149 (MQ=255)
                                                                                                cTTCGCTGACTACGTTGGCATCCGCAAGAGGGAAAAGCAGCGCCGCCCACAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGc                                             <  2:640654/149‑1 (MQ=255)
                                                                                                 ttCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCa                                            <  2:395670/149‑1 (MQ=255)
                                                                                                        acgaCGTTGGCATCCGCCAGAGGGAAAAGCCGCGCCGCCCGGAGCGTGTCATTGTCCATACTCAATGGCGAGAGAATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCTGACAATa                                     <  2:310636/146‑1 (MQ=255)
                                                                                                             gTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCGTGCGTCTGTTGCAGACAATACGccc                                <  1:366596/149‑1 (MQ=255)
                                                                                                              gtGGCATCCGCCAGAGGGAACAGCAGCGCCGCCCACAGCGTGTCAATGTACATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCa                               <  2:481152/148‑1 (MQ=255)
                                                                                                                  cATCCGCCAGAGGGAAAAGCAGCGCCGCCCACAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGGtt                           >  2:538171/1‑149 (MQ=255)
                                                                                                                    tCCGCCAGAGGGAAAAGCAGCGCCGCCCGAAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCCTCCGGATGCCCCTGCGCCTGTTGCAGACACAACTCCCAgcgttc                         >  2:333997/1‑144 (MQ=255)
                                                                                                                          agagGGAAAAGCAGCGCCGCCCACAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAGCACCACGCCACAAT‑AACCGACTGGCATCCGGATGCCCCTGCGTCTGTTGCTGACAATACGCCCAGGTTTCGGCTaa                   <  2:344090/149‑1 (MQ=255)
                                                                                                                               gAAAAGCAGCGCCGCCCGCCGCGTGTCAATGTCCATACTTAATGTCGTGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACt              <  2:624862/149‑1 (MQ=255)
                                                                                                                               gAAAAGCAGCGCCACCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGATTTCGGCTAAGCAcc              >  2:384451/1‑148 (MQ=255)
                                                                                                                                  aaGCAGCGCCGCCCACAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTcac           <  1:60499/149‑1 (MQ=255)
                                                                                                                                  aaGCAGCCCCGCCCGCAGCGTGTCAATGTCCATACTTCATGTCGCGAGGCTCTCCACCATCTCAACACCACGCCACAATCCACCG‑GTGCCATCCGGATGCTCCTGCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTcac           <  1:426899/149‑1 (MQ=255)
                                                                                                                                      agcGCCCCCCGCAGCCTGTCAGTGTCCATATTTAATTTCGTGAGGATATCTGCCATCTCAACACCCCGCGACAAT‑AGGAGTCTGGCAGCCTGATGCCCTTGCTTCTGTTGCAGACAATACGCACAGGTTTCGGCTAAGCACTCACACGa       <  2:38903/149‑1 (MQ=255)
                                                                                                                                          ccgcccacAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCCACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTCCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTCACACGACTTc   >  1:702917/1‑149 (MQ=255)
                                                                                                                                           ggaccgcAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTAAACACCACGCCACAAT‑CACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTCACACGACttat  <  2:210594/146‑3 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACGACCGACTTACCGACGCTCTCACGCAGCACATCTTCGCTGACTACGTTGGCATCCGCCAGAGGGAAAAGCAGCGCCGCCCGCAGCGTGTCAATGTCCATACTTAATGTCGAGAGGATCTCCACCATCTCAACACCACGCCACAAT‑AACAGACTGGCATCCGGATGCCCCTGCGTCTGTTGCAGACAATACGCCCAGGTTTCGGCTAAGCACTCACACGACTTCT  >  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: ‑