Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,851,651 G→T 25.0% R232L (CGG→CTG)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,851,6510GT25.0% 45.2 / 10.4 24R232L (CGG→CTG) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (10/8);  new base T (3/3);  total (13/11)
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.73e-01

TTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGG  >  NC_002947/1851524‑1851794
                                                                                                                               |                                                                                                                                               
ttCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCTGCGGCTGGACCGCGAGCGCTtc                                                                                                                           >  1:426268/1‑150 (MQ=255)
 tCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTTCATGAGCGGCGGCTGGa                                                                                                                                        >  2:295338/1‑136 (MQ=255)
 tCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGa                                                                                                                                        <  1:295338/136‑1 (MQ=255)
    aCATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTtcatca                                                                                                                       <  2:301814/150‑1 (MQ=255)
         tACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCTGCGGCTGGACCGCGAGCGCTTCATCAGCGAc                                                                                                                  <  2:392834/150‑1 (MQ=255)
                                 ggTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAg                                                                                          >  2:290438/1‑150 (MQ=255)
                                     tCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAg                                                                                          <  2:334797/146‑1 (MQ=255)
                                     tCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAg                                                                                          >  1:334797/1‑146 (MQ=255)
                                         tCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCTGCGGCTTGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGc                                                                                  <  1:103490/150‑1 (MQ=255)
                                              tCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGtgat                                                                                                             <  1:17142/118‑1 (MQ=255)
                                              tCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGtgat                                                                                                             >  2:17142/1‑118 (MQ=255)
                                              tCAAGTGGGAGCTGGAAGACCTGTCCTTACGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGa                                                                                                                                        >  2:409592/1‑91 (MQ=255)
                                              tCAAGTGGGAGCTGGAAGACCTGTCCTTACGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGa                                                                                                                                        <  1:409592/91‑1 (MQ=255)
                                                   tGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTg                                                                        >  2:461464/1‑150 (MQ=255)
                                                          tGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCg                                                                 <  1:148171/150‑1 (MQ=255)
                                                                    gTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCTGCGGCTGGACCGCGAGCGCTTCATCAGc                                                                                                                     <  1:390685/88‑1 (MQ=255)
                                                                    gTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCTGCGGCTGGACCGCGAGCGCTTCATCAGc                                                                                                                     >  2:390685/1‑88 (MQ=255)
                                                                        ttCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGGCGACATCAGCTGCCgg                                                   >  2:201793/1‑150 (MQ=255)
                                                                                  tGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCTGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACata                                         >  1:190601/1‑150 (MQ=255)
                                                                                    gAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATc                                       <  1:290438/150‑1 (MQ=255)
                                                                                                            gCCAAGCTGACGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCAACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAg               >  2:20692/1‑150 (MQ=255)
                                                                                                                         aTGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATCTAt                                    <  2:157526/116‑1 (MQ=255)
                                                                                                                         aTGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATCTAt                                    >  1:157526/1‑116 (MQ=255)
                                                                                                                         aTGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCTGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTgg  >  1:202177/1‑150 (MQ=255)
                                                                                                                               |                                                                                                                                               
TTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGG  >  NC_002947/1851524‑1851794

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: