Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,882,489 G→T 53.5% A215A (GCC→GCA relA ← (p)ppGpp synthetase I/GTP pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,882,4890GT53.5% ‑4.0 / 100.4 84A215A (GCC→GCArelA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base G (24/15);  new base T (17/28);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 4.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

GCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTC  >  minE/1882423‑1882552
                                                                  |                                                               
gcgcAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGc                                                               <  1:462879/69‑1 (MQ=255)
 cgcAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGc                                                               <  1:931889/68‑1 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:1499234/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:1960610/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:1958920/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:1953317/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:2527941/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:1748761/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:1744934/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:2685614/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:2918376/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:2953200/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:3033302/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:105451/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:3137745/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:609138/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                              >  1:40370/1‑69 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCg                                               >  1:494796/1‑68 (MQ=255)
                 aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGACGCGTTTGTATTCgg                                              >  1:1684632/1‑69 (MQ=255)
                  ctctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGt                                             >  1:1020897/1‑69 (MQ=255)
                  ctctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGt                                             >  1:1026804/1‑69 (MQ=255)
                  ctctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGt                                             >  1:2556375/1‑69 (MQ=255)
                  ctctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGt                                             >  1:44130/1‑69 (MQ=255)
                  ctctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGt                                             >  1:3075952/1‑69 (MQ=255)
                    ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                           <  1:257554/69‑1 (MQ=255)
                    ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                           <  1:2285505/69‑1 (MQ=255)
                    ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                           <  1:203926/69‑1 (MQ=255)
                    ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                           <  1:572079/69‑1 (MQ=255)
                    ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                           <  1:2942334/69‑1 (MQ=255)
                    ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                           <  1:1227256/69‑1 (MQ=255)
                    ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                           <  1:1587166/69‑1 (MQ=255)
                          tatAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCGGtc                                            <  1:3365970/60‑2 (MQ=255)
                          tatAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCGGtc                                            <  1:635111/60‑2 (MQ=255)
                          tatAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCGGtc                                            <  1:695325/60‑2 (MQ=255)
                          tatAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCGGtc                                            <  1:1258621/60‑2 (MQ=255)
                          tatAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCGGtc                                            <  1:1831722/60‑2 (MQ=255)
                          tatAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATACGTTTGTATTCGGtc                                            <  1:1335950/60‑2 (MQ=255)
                            tAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAg                                    >  1:95214/1‑68 (MQ=255)
                                 ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                               >  1:587110/1‑68 (MQ=255)
                                 ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                               >  1:1015467/1‑68 (MQ=255)
                                 ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                               >  1:1232638/1‑68 (MQ=255)
                                 ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                               >  1:3058751/1‑68 (MQ=255)
                                 ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                               >  1:303488/1‑68 (MQ=255)
                                 ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                               >  1:2492542/1‑68 (MQ=255)
                                 ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                               >  1:2099099/1‑68 (MQ=255)
                                 ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                               >  1:490955/1‑68 (MQ=255)
                                         cGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGt                      >  1:1149590/1‑69 (MQ=255)
                                         cGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGt                      >  1:2202785/1‑69 (MQ=255)
                                         cGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGt                      >  1:718776/1‑69 (MQ=255)
                                          gagaCGCCGTTCATGCAGCAGTTTGTCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTa                     <  1:3299596/69‑1 (MQ=255)
                                          gagaCGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTa                     <  1:727243/69‑1 (MQ=255)
                                          gagaCGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTa                     <  1:1889420/69‑1 (MQ=255)
                                            gaCGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTaa                    <  1:111208/68‑1 (MQ=255)
                                             aCGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTaa                    <  1:704680/67‑1 (MQ=255)
                                                   ttCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAg            <  1:130944/69‑1 (MQ=255)
                                                     cATGCAGCAGTTTTGCGATGCGTTTGTATTCGGtc                                            <  1:3133216/35‑2 (MQ=21)
                                                     cATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtgga                                         >  1:2985487/1‑38 (MQ=255)
                                                         cagcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCa      >  1:1214323/1‑69 (MQ=255)
                                                         cagcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCa      >  1:2380264/1‑69 (MQ=255)
                                                         cagcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCa      >  1:1940723/1‑69 (MQ=255)
                                                         cagcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCa      >  1:3377752/1‑69 (MQ=255)
                                                         cagcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCa      >  1:2261358/1‑69 (MQ=255)
                                                         cagcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCa      >  1:1388132/1‑69 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:969835/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:608927/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:676439/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:932736/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:82555/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:814666/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:41924/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:3425352/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:3378177/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:3357433/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:3325428/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:3250587/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:2906138/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:2720051/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:2676234/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:2043337/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:1816819/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:1752842/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:1737417/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:1582224/69‑1 (MQ=255)
                                                             agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc  <  1:1158111/69‑1 (MQ=255)
                                                                  |                                                               
GCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTC  >  minE/1882423‑1882552

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