Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,852,005 G→A R350H (CGT→CAT)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,852,0050GA92.9% 163.6 / 0.6 56R350H (CGT→CAT) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (25/27);  total (27/29)
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.06e-01

CACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGC  >  NC_002947/1851859‑1852150
                                                                                                                                                  |                                                                                                                                                 
cacCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATAc                                                                                                                                                <  1:67931/150‑1 (MQ=255)
 acCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGGGGTGCAGATCCGTAcc                                                                                                                                               >  2:1611543/1‑150 (MQ=255)
  ccGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGAACACTGCGGTCATCCACCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATAccc                                                                                                                                              >  2:139606/1‑150 (MQ=255)
     cgcTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACg                                                                                                                                           <  2:1086542/150‑1 (MQ=255)
     cgcTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGATTTAATCCTCCCCGAGGGCAAGGTGCTGGGGGTGCAGATCCGTACCCACg                                                                                                                                           <  1:1611543/150‑1 (MQ=255)
             aTCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCAc                                                                                                                                   >  2:1096563/1‑150 (MQ=255)
             aTCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACACCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCAc                                                                                                                                   >  2:1400287/1‑150 (MQ=255)
              tCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACg                                                                                                                                  >  1:1009391/1‑150 (MQ=255)
                    acacCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAgg                                                                                                                            >  2:123753/1‑150 (MQ=255)
                             ggcggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTg                                                                                                                   <  2:1414502/150‑1 (MQ=255)
                              gcggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCACGGAACCCGGCCACCCTGCGTTGCACACTGCGTTAATCGTCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTgg                                                                                                                  <  2:635894/150‑1 (MQ=255)
                              gcggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTgg                                                                                                                  >  2:1032028/1‑150 (MQ=255)
                              gcggcATATTCCCAAGGAGTTCGACGACTACACCGCCAGCCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTgg                                                                                                                  <  1:1400287/150‑1 (MQ=255)
                                ggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTg                                                                                                                   >  1:594181/1‑147 (MQ=255)
                                ggcATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTg                                                                                                                   <  2:594181/147‑1 (MQ=255)
                                    tatTCCCAAGGAGTTCGACGACTACATGGGCAACCCCAAGGAGAACGGCTACCGCTCGTTGCGCACTGCGGTAATCGGCCCGGAGGGCAAGGGGCTGGAGGGGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCatata                                                                                                            >  2:33521/1‑149 (MQ=255)
                                    tatTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGc                                                                                                                 <  2:1357537/145‑1 (MQ=255)
                                    tatTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGc                                                                                                                 >  1:1357537/1‑145 (MQ=255)
                                      ttCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATgcg                                                                                                          >  2:1339698/1‑150 (MQ=255)
                                           aaGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCAc                                                                                                     >  2:829676/1‑150 (MQ=255)
                                           aaGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCAc                                                                                                     <  1:1096563/150‑1 (MQ=255)
                                               aGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGgc                                                                                                 <  2:313968/150‑1 (MQ=255)
                                                     acgacTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACa                                                                                           >  2:689898/1‑150 (MQ=255)
                                                      cgacTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACaa                                                                                          <  2:1663220/150‑1 (MQ=255)
                                                       gacTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAg                                                                                         >  1:1108632/1‑150 (MQ=255)
                                                       gacTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAg                                                                                         <  1:689898/150‑1 (MQ=255)
                                                       gacTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAg                                                                                         <  1:137930/150‑1 (MQ=255)
                                                       gacTACATCCCTACCCGCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCGCTGGTTGCAAGGTGCTGCAGGTGCAGATCCATCCCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAg                                                                                         <  1:123753/150‑1 (MQ=255)
                                                                gCCAACCCCAAGGAGAACGGCCACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGAc                                                                                >  2:56841/1‑150 (MQ=255)
                                                                       ccAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGc                                                                         >  2:650827/1‑150 (MQ=255)
                                                                        cAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGcc                                                                        <  1:1623240/150‑1 (MQ=255)
                                                                        cAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGcc                                                                        <  1:1339698/150‑1 (MQ=255)
                                                                        cAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGcc                                                                        >  2:1623240/1‑150 (MQ=255)
                                                                         aaGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGccc                                                                       <  1:539200/150‑1 (MQ=255)
                                                                           ggAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAg                                                                     <  2:805191/150‑1 (MQ=255)
                                                                                 cGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCaa                                                               <  2:1008671/150‑1 (MQ=255)
                                                                                  ggCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAAc                                                              <  1:99717/150‑1 (MQ=255)
                                                                                       ccGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGc                                                                         >  2:970623/1‑134 (MQ=255)
                                                                                       ccGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGc                                                                         <  1:970623/134‑1 (MQ=255)
                                                                                          cTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACga                                                      >  1:418803/1‑150 (MQ=255)
                                                                                                gCACACTGCGGTAATCGGCCCCGAGGTCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCCACGACGTCAAGCCCAGCTCCAACCACTACGAAGaaaa                                                <  1:139606/150‑1 (MQ=255)
                                                                                                      tGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTc                                          >  2:1032841/1‑150 (MQ=255)
                                                                                                           tAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGt                                     >  1:21745/1‑150 (MQ=255)
                                                                                                              tCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGc                                  <  1:829676/150‑1 (MQ=255)
                                                                                                               cGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCg                                 <  1:56841/150‑1 (MQ=255)
                                                                                                                 gCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTCTGTGCTCACTGGCGCTACAAGTTCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTc                               <  1:650827/150‑1 (MQ=255)
                                                                                                                   cccGAGGGCAAGGTGCTGGAGGTGCCGATCCATCCCGATGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAg                             >  1:139589/1‑150 (MQ=255)
                                                                                                                   cccGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAg                             <  1:1032841/150‑1 (MQ=255)
                                                                                                                    ccGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAgg                            >  2:1302627/1‑150 (MQ=255)
                                                                                                                    ccGAGGGCAAGGGGCTGGAGGTGCCGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAgg                            >  1:617400/1‑150 (MQ=255)
                                                                                                                         ggCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTg                       <  1:1592482/150‑1 (MQ=255)
                                                                                                                                tGCTGGGGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGc                <  1:1554041/150‑1 (MQ=255)
                                                                                                                                   tGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAg                                               <  2:187640/116‑1 (MQ=255)
                                                                                                                                   tGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAg                                               >  1:187640/1‑116 (MQ=255)
                                                                                                                                    ggAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTCGCACGa            >  2:170565/1‑150 (MQ=255)
                                                                                                                                              aTCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGc  >  1:563753/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGGC  >  NC_002947/1851859‑1852150

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