Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,852,053 G→T W366L (TGG→TTG)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,852,0530GT98.0% 147.6 / ‑5.3 51W366L (TGG→TTG) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (26/24);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

TCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTG  >  NC_002947/1851909‑1852196
                                                                                                                                                |                                                                                                                                                
tCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCTCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCt                                                                                                                                             >  2:615940/1‑150 (MQ=255)
 cgacgaCTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTa                                                                                                                                            <  1:750666/150‑1 (MQ=255)
     gacTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACGGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAg                                                                                                                                        <  1:28766/150‑1 (MQ=255)
       cgaCATCGCCAACCCCAAGGAGCACGGCTACCGCTCGTTGCACACCGCGTTAATCGGCCCCGAGGGCAAGGGGCTGGAGGAGCAGATCAGTGCCCACGGAATGCAAGAAGAGGCCGAACTTGGCGTATGCGCGCACTGGCGCTACAAggg                                                                                                                                      >  2:265231/3‑150 (MQ=255)
        tACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCAGACTGCGGTAATCAGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAAGTGCCGAACATAGCGAATGCGCCCCATTGCGCTACAAGGGc                                                                                                                                     >  2:116588/1‑150 (MQ=255)
        tACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGc                                                                                                                                     >  2:693280/1‑150 (MQ=255)
           aTCGCCAACCCCAAGGAGAACGGCTACGGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCAcc                                                                                                                                  >  1:578410/1‑150 (MQ=255)
              gCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGAc                                                                                                                               >  2:411239/1‑150 (MQ=255)
                     ccAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGc                                                                                                                        >  2:649728/1‑150 (MQ=255)
                        aGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCa                                                                                                                     <  1:411239/150‑1 (MQ=255)
                         ggAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAg                                                                                                                    >  1:941290/1‑150 (MQ=255)
                            gaACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTc                                                                                                                 >  2:160413/1‑150 (MQ=255)
                                                      cGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAGAAGATTTCCt                                                                                      <  1:160413/150‑1 (MQ=255)
                                                        gTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTgg                                                                                    <  1:649728/150‑1 (MQ=255)
                                                          aaTCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGaaaa                                                                                              <  1:73306/138‑1 (MQ=255)
                                                          aaTCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGaaaa                                                                                              >  2:73306/1‑138 (MQ=255)
                                                          aaTCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGtt                                                                                  <  1:693280/150‑1 (MQ=255)
                                                           aTCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTg                                                                                 >  1:389220/1‑150 (MQ=255)
                                                                  ccGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCACACGTTCGCTACAAGGGCACCGACGTCAAGACAACATCTAACTACCACTAAGAAAAGATTTCCTGGTTGCGTCAgg                                                                          >  1:707612/1‑150 (MQ=255)
                                                                      gggCAAGGTGCTGGAGGCGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCt                                                                      >  1:641776/1‑150 (MQ=255)
                                                                        gCAAGGTGCTGGATCTGCAGAGCCGTACCCACGGCATGCACGAAGAGGCCGAACTTTGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTgg                                                                    <  2:389220/150‑1 (MQ=255)
                                                                         cAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGa                                                                   <  1:318925/150‑1 (MQ=255)
                                                                          aaGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCAATTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAg                                                                  >  2:359905/1‑150 (MQ=255)
                                                                            ggTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTg                                                                <  2:336143/150‑1 (MQ=255)
                                                                                 tGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACg                                                           >  1:732152/1‑150 (MQ=255)
                                                                                        gCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCt                                                    <  1:615940/150‑1 (MQ=255)
                                                                                        gCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCt                                                    >  1:166490/1‑150 (MQ=255)
                                                                                           gATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACgaa                                                                                                    >  2:774100/1‑100 (MQ=255)
                                                                                           gATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACgaa                                                                                                    <  1:774100/100‑1 (MQ=255)
                                                                                            aTCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGc                                                <  2:472107/150‑1 (MQ=255)
                                                                                                    ccACGGCATACACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCgg                                        >  1:273791/1‑150 (MQ=255)
                                                                                                     cACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGTACGAAGAGCTGGGCGACATCGGt                                       <  1:670814/150‑1 (MQ=255)
                                                                                                        ggCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTg                                                                                     <  1:544936/101‑1 (MQ=255)
                                                                                                        ggCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTg                                                                                     >  2:544936/1‑101 (MQ=255)
                                                                                                                 gAAGAGGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAg                           >  1:663961/1‑150 (MQ=255)
                                                                                                                   agagGCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCa                         <  2:273791/150‑1 (MQ=255)
                                                                                                                       gCCGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGGCGAAAGAGATTCTTGTTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTg                     >  2:526937/1‑150 (MQ=255)
                                                                                                                         cGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTTGGCGACATCGGTGGCCTGGCCGAGCAGTTGCg                   >  2:176344/1‑150 (MQ=255)
                                                                                                                         cGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGt                                       >  2:377935/1‑130 (MQ=255)
                                                                                                                         cGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGt                                       <  1:377935/130‑1 (MQ=255)
                                                                                                                         cGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCg                   >  1:970217/1‑150 (MQ=255)
                                                                                                                         cGAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCg                   >  2:234661/1‑150 (MQ=255)
                                                                                                                          gAACTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCgg                  <  2:970217/150‑1 (MQ=255)
                                                                                                                            aCTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGt                <  2:732152/150‑1 (MQ=255)
                                                                                                                             cTTGGCGTATGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTc               <  2:166490/150‑1 (MQ=255)
                                                                                                                                     aTGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGa       <  1:639481/150‑1 (MQ=255)
                                                                                                                                     aTGCGCCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGa       <  1:682344/150‑1 (MQ=255)
                                                                                                                                       gcgcCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGc     <  2:641776/150‑1 (MQ=255)
                                                                                                                                        cgcCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGc     <  1:972192/149‑1 (MQ=255)
                                                                                                                                        cgcCCACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCAATACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGc     >  2:972192/1‑149 (MQ=255)
                                                                                                                                          cccACTTGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTg  >  2:164572/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGCGACATCGGTGGCCTGGCCGAGCAGTTGCGGGTCGACATCGAGCCTG  >  NC_002947/1851909‑1852196

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: