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,0050GA96.4% 175.5 / ‑2.3 56R350H (CGT→CAT) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (29/24);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

CATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGG  >  NC_002947/1851871‑1852149
                                                                                                                                      |                                                                                                                                                
cATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCa                                                                                                                                   <  2:726152/150‑1 (MQ=255)
     tGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAag                                                                                                                              >  2:253123/1‑150 (MQ=255)
        acacCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAgg                                                                                                                           >  2:1039329/1‑150 (MQ=255)
        acacCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAgg                                                                                                                           <  1:321213/150‑1 (MQ=255)
        acacCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAgg                                                                                                                           >  2:203057/1‑150 (MQ=255)
                       atatTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTAt                                                                                                            <  1:592311/150‑1 (MQ=255)
                             ccAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCa                                                                                                                                       <  1:268701/117‑1 (MQ=255)
                             ccAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCa                                                                                                                                       >  2:268701/1‑117 (MQ=255)
                                           gacTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAg                                                                                        >  1:1180430/1‑150 (MQ=255)
                                           gacTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAg                                                                                        <  1:203057/150‑1 (MQ=255)
                                            acTACATCGACAACCCAACGGAGCAGGGCTAGCGAGGGGGCCATACTGCTGGTAGCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAgg                                                                                       <  2:759111/150‑1 (MQ=255)
                                             cTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAggg                                                                                      >  1:715456/1‑150 (MQ=255)
                                               aCATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCa                                                                                    >  2:267294/1‑150 (MQ=255)
                                                cATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCAc                                                                                   >  2:995326/1‑150 (MQ=255)
                                                 aTCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCAcc                                                                                  >  2:29753/1‑150 (MQ=255)
                                                 aTCGCCAACCCCAAGGAGAACGGCTACCGCTCGGTGCACACTGCGGTAATCAGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGCCATATACGCCCACTGGCGATaaaaccggcca                                                                                  >  2:688752/1‑143 (MQ=255)
                                                     ccAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACg                                                                              >  1:1417846/1‑150 (MQ=255)
                                                         ccccAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATTACTTGCCATGGCAGGCACGAAGAGGCCGAACTTGGCGTATGAGCCCACTGGCGCTACAAGGGCACCGACGTCaa                                                                          >  1:403842/1‑150 (MQ=255)
                                                           ccAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGc                                                                        >  2:176922/1‑150 (MQ=255)
                                                           ccAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGc                                                                        >  1:958312/1‑150 (MQ=255)
                                                                     cGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCaa                                                              <  1:267294/150‑1 (MQ=255)
                                                                     cGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCaa                                                              <  2:715456/150‑1 (MQ=255)
                                                                         tACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCAc                                                          <  1:162371/150‑1 (MQ=255)
                                                                           ccGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGATGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGCATGTGCGTACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTa                                                        <  1:1150474/150‑1 (MQ=255)
                                                                           ccGCTCGTTGCACACTGCAGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTa                                                        >  2:1191191/1‑150 (MQ=255)
                                                                               tCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACgaa                                                    >  2:1350759/1‑150 (MQ=255)
                                                                                 gTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGaa                                                                                                                      <  2:1361786/82‑1 (MQ=255)
                                                                                 gTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGaa                                                                                                                      >  1:1361786/1‑82 (MQ=255)
                                                                                            cGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGc                                                                                                                >  2:1311503/1‑77 (MQ=255)
                                                                                            cGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGc                                                                                                                <  1:1311503/77‑1 (MQ=255)
                                                                                              gTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGcc                                                                       >  1:76564/1‑116 (MQ=255)
                                                                                              gTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGcc                                                                       <  2:76564/116‑1 (MQ=255)
                                                                                                 aTCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTa                                                                                            >  1:207522/1‑92 (MQ=255)
                                                                                                 aTCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTa                                                                                            <  2:207522/92‑1 (MQ=255)
                                                                                                  tCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGc                                 <  2:958312/150‑1 (MQ=255)
                                                                                                     gCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGTCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTc                                         <  2:1312763/139‑1 (MQ=255)
                                                                                                     gCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACgaag                                                   >  1:586985/1‑129 (MQ=255)
                                                                                                     gCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACgaag                                                   <  2:586985/129‑1 (MQ=255)
                                                                                                     gCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTc                                         >  1:1312763/1‑139 (MQ=255)
                                                                                                       cccGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAg                            >  1:1322959/1‑150 (MQ=255)
                                                                                                        ccGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAgg                           >  2:645886/1‑150 (MQ=255)
                                                                                                         cGAGGGCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGt                          >  1:182131/1‑150 (MQ=255)
                                                                                                             ggCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTg                      <  1:253123/150‑1 (MQ=255)
                                                                                                             ggCAAGGTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTg                      <  2:1417846/150‑1 (MQ=255)
                                                                                                               cAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGcc                                                                       >  1:364742/1‑99 (MQ=255)
                                                                                                               cAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGcc                                                                       <  2:364742/99‑1 (MQ=255)
                                                                                                                  ggTGCTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTg                 >  1:1003867/1‑150 (MQ=255)
                                                                                                                      cTGGAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCAc             >  1:560050/1‑150 (MQ=255)
                                                                                                                        ggAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCAc                                                                                                    >  2:136091/1‑61 (MQ=255)
                                                                                                                        ggAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCAc                                                                                                    <  1:136091/61‑1 (MQ=255)
                                                                                                                         gAGGTGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGaa          <  2:182131/150‑1 (MQ=255)
                                                                                                                             tGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCt                                       <  1:1319462/117‑1 (MQ=255)
                                                                                                                             tGCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCt                                       >  2:1319462/1‑117 (MQ=255)
                                                                                                                              gCAGATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCt     <  2:403842/150‑1 (MQ=255)
                                                                                                                                 gATCCATACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTggg  <  2:1180430/150‑1 (MQ=255)
                                                                                                                                 gATCCATACACACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTggg  >  2:985215/1‑150 (MQ=255)
                                                                                                                                      |                                                                                                                                                
CATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGCTACAAGGGCACCGACGTCAAGCCCAGCTCCAACCACTACGAAGAAAAGATTTCCTGGTTGCGTCAGGTGCTGGAGTGGCACGAAGAGCTGGG  >  NC_002947/1851871‑1852149

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