Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,882,495 T→G 52.2% R213R (CGA→CGC 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,4950TG52.2% 15.5 / 92.7 86R213R (CGA→CGCrelA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base T (28/13);  new base G (17/27);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 8.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-01

ACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGT  >  minE/1882440‑1882558
                                                       |                                                               
aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                                    >  1:609138/1‑69 (MQ=255)
aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                                    >  1:1744934/1‑69 (MQ=255)
aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                                    >  1:1748761/1‑69 (MQ=255)
aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                                    >  1:1953317/1‑69 (MQ=255)
aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                                    >  1:1958920/1‑69 (MQ=255)
aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                                    >  1:1499234/1‑69 (MQ=255)
aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                                    >  1:1960610/1‑69 (MQ=255)
aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                                    >  1:2527941/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:3137745/1‑69 (MQ=255)
aCTCTTCGATATAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCgg                                                    >  1:105451/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:2556375/1‑69 (MQ=255)
 ctctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGt                                                   >  1:3075952/1‑69 (MQ=255)
 ctctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGt                                                   >  1:44130/1‑69 (MQ=255)
 ctctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGt                                                   >  1:1026804/1‑69 (MQ=255)
   ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                                 <  1:572079/69‑1 (MQ=255)
   ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                                 <  1:1227256/69‑1 (MQ=255)
   ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                                 <  1:2942334/69‑1 (MQ=255)
   ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                                 <  1:1587166/69‑1 (MQ=255)
   ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                                 <  1:257554/69‑1 (MQ=255)
   ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                                 <  1:203926/69‑1 (MQ=255)
   ctTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTtg                                                 <  1:2285505/69‑1 (MQ=255)
         tatAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCGGtc                                                  <  1:3365970/60‑2 (MQ=255)
         tatAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCGGtc                                                  <  1:1258621/60‑2 (MQ=255)
         tatAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCGGtc                                                  <  1:635111/60‑2 (MQ=255)
         tatAGTGTTCGCGGTCGAGACGGCGTTCATGCAGCAGTTTTGCGATGCGTTTGTATTCGGtc                                                  <  1:695325/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:2099099/1‑68 (MQ=255)
                ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                                     >  1:587110/1‑68 (MQ=255)
                ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                                     >  1:1015467/1‑68 (MQ=255)
                ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                                     >  1:2492542/1‑68 (MQ=255)
                ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                                     >  1:3058751/1‑68 (MQ=255)
                ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                                     >  1:303488/1‑68 (MQ=255)
                ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                                     >  1:490955/1‑68 (MQ=255)
                ttCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAAc                                     >  1:1232638/1‑68 (MQ=255)
                        cGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGt                            >  1:2202785/1‑69 (MQ=255)
                        cGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGt                            >  1:1149590/1‑69 (MQ=255)
                        cGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGt                            >  1:718776/1‑69 (MQ=255)
                         gagaCGCCGTTCATGCAGCAGTTTGTCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTa                           <  1:3299596/69‑1 (MQ=255)
                         gagaCGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTa                           <  1:1889420/69‑1 (MQ=255)
                         gagaCGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTa                           <  1:727243/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:3377752/1‑69 (MQ=255)
                                        cagcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCa            >  1:2261358/1‑69 (MQ=255)
                                        cagcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCa            >  1:2380264/1‑69 (MQ=255)
                                        cagcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCa            >  1:1214323/1‑69 (MQ=255)
                                        cagcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCa            >  1:1388132/1‑69 (MQ=255)
                                        cagcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCa            >  1:1940723/1‑69 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:82555/69‑1 (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:814666/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:2676234/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:1158111/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:1582224/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:1737417/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:1752842/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:1816819/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:2043337/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:41924/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:2720051/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:2906138/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:3250587/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:3325428/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:3357433/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:3378177/69‑1 (MQ=255)
                                            agTTTTGCGATGCGTTTGTATTCGGTCGGGTGCAGATAACGGAAGCAGTAATCTTCCAGCTCCCATTTc        <  1:3425352/69‑1 (MQ=255)
                                                  gCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGt  >  1:2763577/1‑69 (MQ=255)
                                                  gCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGt  >  1:2573965/1‑69 (MQ=255)
                                                  gCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGt  >  1:1259529/1‑69 (MQ=255)
                                                  gCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTACAGTTCCCATTTCAGTTg   >  1:2082731/1‑68 (MQ=255)
                                                       |                                                               
ACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGT  >  minE/1882440‑1882558

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: