Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,851,242 G→T 100% E96* (GAG→TAG)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,851,2420GT93.3% 32.6 / ‑3.4 15E96* (GAG→TAG) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (8/6);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTCGCA  >  NC_002947/1851095‑1851375
                                                                                                                                                   |                                                                                                                                     
gTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAg                                                                                                                                     >  2:251091/1‑150 (MQ=255)
                 cccGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGc                                                                                                                    <  1:251091/150‑1 (MQ=255)
                                cTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACGCCTCGGTCCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCg                                                                                                     >  2:248901/1‑150 (MQ=255)
                                                                     gAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATc                                                                         <  1:444580/141‑1 (MQ=255)
                                                                     gAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCGGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATc                                                                         >  2:444580/1‑141 (MQ=255)
                                                                                             cTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCggtg                                                                                        >  1:413782/1‑102 (MQ=255)
                                                                                             cTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCggtg                                                                                        <  2:413782/102‑1 (MQ=255)
                                                                                                  gCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTcagcc                                                                                                         >  1:538996/1‑80 (MQ=255)
                                                                                                  gCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTcagcc                                                                                                         <  2:538996/80‑1 (MQ=255)
                                                                                                        ccAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGccc                             >  2:272790/1‑150 (MQ=255)
                                                                                                        ccAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGccc                             >  1:267315/1‑150 (MQ=255)
                                                                                                           ggACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATc                                                                         >  1:120494/1‑103 (MQ=255)
                                                                                                           ggACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATc                                                                         <  2:120494/103‑1 (MQ=255)
                                                                                                                                  tCTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTcgc   >  2:44376/1‑150 (MQ=255)
                                                                                                                                   cTACCGCTCGGTGCGCTAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTcgca  <  1:248901/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                     
GTAGAGAAAAAGGGCAACCCGGCCAAGCATTCCTGGGCGGATGGTACGTCCAGCTTCCAGGCAGGCCTGGAAATCGCCGAAATTCTGGCTGACCTCAAGCTCGACCAGGACTCCCTGGTGGCTGCGGTCATCTACCGCTCGGTGCGCGAGGGCAAGGTCACCCTCGCCGAGGTCAGCCAGCGGTTTGGCCCGGTGGTGTCCAAGCTGATCGACGGTGTGCTGCGCATGGCCGCCATCAGTGCCAGCCTCAGCCCACGACAGTCGCTGGTGCTGGGCTCGCA  >  NC_002947/1851095‑1851375

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: