Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,851,912 A→G 100% D319G (GAC→GGC)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,851,9120AG97.2% 88.1 / ‑5.9 29D319G (GAC→GGC) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (16/12);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGC  >  NC_002947/1851771‑1852057
                                                                                                                                             |                                                                                                                                                 
gcAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACa                                                                                                                                           >  2:519592/1‑150 (MQ=255)
         aGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCCCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTACAGATCGTCTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCCAcc                                                                                                                                  <  1:573801/150‑1 (MQ=255)
          gcgcAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAAccc                                                                                                                                 >  2:416848/1‑150 (MQ=255)
            gcAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTACGGTGCGCGTGCTGGTGCCGGAAATCCGGGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAAccccc                                                                                                                               >  2:356539/1‑149 (MQ=255)
              aaaGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAg                                                                                                                             >  2:57069/1‑150 (MQ=255)
                     tGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACg                                                                                                                      >  2:500317/1‑150 (MQ=255)
                     tGGACTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACg                                                                                                                      >  2:500322/1‑150 (MQ=255)
                      ggAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCggcgac                                                                                                                                               <  1:207141/124‑1 (MQ=255)
                      ggAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCggcgac                                                                                                                                               >  2:207141/1‑124 (MQ=255)
                         gTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTa                                                                                                                  <  1:331711/150‑1 (MQ=255)
                            cAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTACCg                                                                                                               <  1:519592/150‑1 (MQ=255)
                             aGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTACCGc                                                                                                              >  2:115039/1‑150 (MQ=255)
                               ccAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTc                                                                                                            <  1:122340/150‑1 (MQ=255)
                                  gATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCCCCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGtt                                                                                                         <  1:500322/150‑1 (MQ=255)
                                  gATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCAAAGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGtt                                                                                                         <  1:500317/150‑1 (MQ=255)
                                   aTCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCTCCAACCCCAAg                                                                                                                             <  2:183111/129‑1 (MQ=255)
                                   aTCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCCACCCCAAGGAGAACGGCTACCGCTCGTTg                                                                                                        <  1:420837/150‑1 (MQ=255)
                                   aTCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAg                                                                                                                             >  1:183111/1‑129 (MQ=255)
                                      tACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAg                                                                                                                             >  2:198948/1‑126 (MQ=255)
                                      tACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAg                                                                                                                             <  1:198948/126‑1 (MQ=255)
                                             tgcgtgcgGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCgg                                                                                              >  1:254313/1‑150 (MQ=255)
                                                           gTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAg                                                                                <  1:115039/150‑1 (MQ=255)
                                                                              gcgACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCt                                                                                                                   <  1:184938/96‑1 (MQ=255)
                                                                              gcgACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCt                                                                                                                   >  2:184938/1‑96 (MQ=255)
                                                                                      tACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCgg                                                                                       <  2:91498/116‑1 (MQ=255)
                                                                                      tACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCgg                                                                                       >  1:91498/1‑116 (MQ=255)
                                                                                                                cTTGTGGCGGCATATTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCACGGTGCACACTCCGGCAACCCGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGa                           >  1:44096/1‑150 (MQ=255)
                                                                                                                            tatTCCCAAGGAGTTCGGCGACTACATCGCCAACCCCGAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATg               >  2:215475/1‑150 (MQ=255)
                                                                                                                              ttCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATgcg             >  1:52038/1‑150 (MQ=255)
                                                                                                                                 ccAAGGAGTTCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGGGTATgggcgc          >  2:141931/1‑148 (MQ=255)
                                                                                                                                         ttCGGCGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGAGCAAGGTGGTGGGGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGgcgc  >  1:382603/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
GCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTCCCAAGGAGTTCGACGACTACATCGCCAACCCCAAGGAGAACGGCTACCGCTCGTTGCACACTGCGGTAATCGGCCCCGAGGGCAAGGTGCTGGAGGTGCAGATCCGTACCCACGGCATGCACGAAGAGGCCGAACTTGGCGTATGCGCCCACTGGCGC  >  NC_002947/1851771‑1852057

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: