Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,851,762 C→T 43.7% S269L (TCG→TTG)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,851,7620CT43.7% 0.3 / 15.2 16S269L (TCG→TTG) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base T (2/5);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

AAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAG  >  NC_002947/1851623‑1851904
                                                                                                                                           |                                                                                                                                              
aaGCAGATCGCTAAGCTGTTGCATGAGCGGCGGCTAGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGgcgc                                                                                                                                      >  2:339262/1‑150 (MQ=255)
          ccAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAgc                                                                                                                            <  1:20179/150‑1 (MQ=255)
                       tGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGa                                                                                                               <  1:420945/150‑1 (MQ=255)
                       tGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGa                                                                                                               <  1:42421/150‑1 (MQ=255)
                                                           cGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGt                                                                           <  2:20962/150‑1 (MQ=255)
                                                                            aGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTTGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATcc                                                          <  1:160324/150‑1 (MQ=255)
                                                                                      gCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTTGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTa                                                <  1:156280/150‑1 (MQ=255)
                                                                                        tCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTaca                                              >  1:92871/1‑150 (MQ=255)
                                                                                                   gTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTTGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGc                                   >  2:432868/1‑150 (MQ=255)
                                                                                                   gTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGc                                   >  2:412321/1‑150 (MQ=255)
                                                                                                      aaGGCCGACATCAGCGGTCGGGCGAAACATATCTATTTGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATc                                <  2:27298/150‑1 (MQ=255)
                                                                                                                   gCGGCCGGGCGAAACATATCTATTTGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTtgt                   <  2:256220/150‑1 (MQ=255)
                                                                                                                             gAAACATATCTATTTGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCatat         >  2:201901/1‑150 (MQ=255)
                                                                                                                               aaCATATCTATTTGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTc       <  2:227407/150‑1 (MQ=255)
                                                                                                                                aCATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGGGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTcc      >  2:474928/1‑150 (MQ=255)
                                                                                                                                    atCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAg  <  1:339262/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
AAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAG  >  NC_002947/1851623‑1851904

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