Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,856,513 A→G 100% D203G (GAT→GGT)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,856,5130AG96.0% 83.0 / ‑4.1 25D203G (GAT→GGT) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (11/13);  total (11/14)
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.87e-01
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.

GTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTT  >  NC_002947/1856380‑1856621
                                                                                                                                     |                                                                                                            
gTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTc                                                                                                            <  1:256998/136‑1 (MQ=255)
gTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTc                                                                                                            <  1:462416/136‑1 (MQ=255)
 gggCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTca                                                                                                           <  1:41110/135‑1 (MQ=255)
          tCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTg                                                                                                  <  2:18830/136‑1 (MQ=255)
                     gAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAg                                                                                       >  1:336992/1‑136 (MQ=255)
                      aaCGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGa                                                                                      >  2:436335/1‑136 (MQ=255)
                               gcgcAATCCGGGCGGTCAAGTCCGCCGATGACTAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCtt                                                                             <  1:397433/136‑1 (MQ=255)
                                              tCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGcc                                                              <  1:541274/136‑1 (MQ=255)
                                               cAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGccc                                                             <  1:55823/136‑1 (MQ=255)
                                                  gTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGAc                                                          >  2:129783/1‑136 (MQ=255)
                                                          aTGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACaa                                                  <  2:402680/136‑1 (MQ=255)
                                                                            gtgtCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGtt                                <  1:489577/136‑1 (MQ=255)
                                                                                   gcgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAg                         <  2:480877/136‑1 (MQ=255)
                                                                                   gcgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAg                         >  2:317100/1‑136 (MQ=255)
                                                                                   gcgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAg                         >  2:314088/1‑136 (MQ=255)
                                                                                   gcgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAg                         <  2:266572/136‑1 (MQ=255)
                                                                                       gAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAgcggc                     <  2:477295/136‑1 (MQ=255)
                                                                                                gACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGc                                    >  1:460467/1‑112 (MQ=255)
                                                                                                gACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGc                                    <  2:460467/112‑1 (MQ=255)
                                                                                                gACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGAcc            >  2:329413/1‑136 (MQ=255)
                                                                                                    tCTACGCGCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGa        <  2:336992/136‑1 (MQ=255)
                                                                                                        cgcgCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAgcgc    >  2:32895/1‑136 (MQ=255)
                                                                                                        cgcgCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAgcgc    >  2:335892/1‑136 (MQ=255)
                                                                                                        cgcgCCGCTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAgcgc    >  2:385275/1‑136 (MQ=255)
                                                                                                          cgccgcTCGCGCACCGCCTGGGTATCGGTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCtt  >  1:236174/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                            
GTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTT  >  NC_002947/1856380‑1856621

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