Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,882,402 G→A 30.7% G244G (GGC→GGT 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,4020GA30.7% 92.5 / 41.9 94G244G (GGC→GGTrelA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base G (20/45);  new base A (23/6);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 1.52e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAG  >  minE/1882341‑1882467
                                                             |                                                                 
tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                             <  1:2200608/68‑1 (MQ=255)
tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                             <  1:1078158/68‑1 (MQ=255)
tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                             <  1:2900846/68‑1 (MQ=255)
tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                             <  1:196666/68‑1 (MQ=255)
tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                             <  1:809253/68‑1 (MQ=255)
tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                             <  1:1710023/68‑1 (MQ=255)
tCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCa                                                             <  1:643976/68‑1 (MQ=255)
tCTGAATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGAACATACACTTCCGCTTTAACGCCTTCa                                                             <  1:2664282/68‑1 (MQ=255)
 cTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                                            <  1:2521838/68‑1 (MQ=255)
 cTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                                            <  1:3031318/68‑1 (MQ=255)
 cTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                                            <  1:688625/68‑1 (MQ=255)
      ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTc                                                       <  1:1082936/68‑1 (MQ=255)
      ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTc                                                       <  1:2172306/68‑1 (MQ=255)
      ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTc                                                       <  1:199599/68‑1 (MQ=255)
      ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTc                                                       <  1:3267771/68‑1 (MQ=255)
      ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                                                      >  1:1578815/1‑69 (MQ=255)
      ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                                                      >  1:2714847/1‑69 (MQ=255)
      ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                                                      >  1:2781318/1‑69 (MQ=255)
      ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                                                      >  1:247757/1‑69 (MQ=255)
      ttttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCa                                                      >  1:1583189/1‑69 (MQ=255)
        ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCAt                                                     >  1:1484043/1‑68 (MQ=255)
        ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCAt                                                     >  1:2563348/1‑68 (MQ=255)
        ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCAt                                                     >  1:2415964/1‑68 (MQ=255)
        ttACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCAt                                                     >  1:2252410/1‑68 (MQ=255)
             ccAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAgc                                               >  1:3351888/1‑69 (MQ=255)
             ccAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAgc                                               >  1:3315393/1‑69 (MQ=255)
             ccAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAgc                                               >  1:1276040/1‑69 (MQ=255)
                 aTGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAg                                                            <  1:1875347/52‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:2499759/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:2594890/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:1429997/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:863065/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:531738/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:2312177/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:262711/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:875258/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:3024413/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:1278626/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:1882531/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:1878815/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:1285844/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:1798290/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:175995/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:1744129/69‑1 (MQ=255)
                  tGCTGTAGATGTGTTTCGGACGACCATACACTTACGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:2426769/69‑1 (MQ=255)
                   gCTGTAGATGTGTTTCGGTCGACCATACACTTCAGCTTTAACACCTACAGCTTTCATCTCAGCGCGCa                                          <  1:1157505/68‑1 (MQ=255)
                   gCTGTAGATGTGTTTCGGACGACCATAGACTTCCGCTTTAACGCCTTCa                                                             <  1:2582384/49‑1 (MQ=255)
                            gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:622342/59‑1 (MQ=255)
                            gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:965747/59‑1 (MQ=255)
                            gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:3369300/59‑1 (MQ=255)
                            gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:969798/59‑1 (MQ=255)
                            gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:2515459/59‑1 (MQ=255)
                            gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:2339363/59‑1 (MQ=255)
                            gtgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCa                                          <  1:1625744/59‑1 (MQ=255)
                              gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                              >  1:1486866/1‑67 (MQ=255)
                              gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                              >  1:3334062/1‑67 (MQ=255)
                              gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                              >  1:163754/1‑67 (MQ=255)
                              gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                              >  1:2025154/1‑67 (MQ=255)
                              gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                              >  1:3396844/1‑67 (MQ=255)
                              gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                              >  1:3423460/1‑67 (MQ=255)
                              gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                              >  1:1224723/1‑67 (MQ=255)
                              gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                              >  1:2155422/1‑67 (MQ=255)
                              gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAcaa                              >  1:942004/1‑67 (MQ=255)
                              gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAca                               >  1:1220239/1‑67 (MQ=255)
                              gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGAca                               >  1:3308541/1‑67 (MQ=255)
                              gtTTCGGTCGACCATACACTTCAGCTTTAACACCTTCAGCTTTCAACTCAGCGCGCAGATGGCCGAcaa                              >  1:1750404/1‑67 (MQ=255)
                                     acgacCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCa                                  <  1:1371645/58‑1 (MQ=255)
                                              cacTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAgtg              >  1:1543971/1‑69 (MQ=255)
                                              cacTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAgtg              >  1:923619/1‑69 (MQ=255)
                                                  tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTcg          >  1:1205101/1‑69 (MQ=255)
                                                  tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTcg          >  1:3196501/1‑69 (MQ=255)
                                                  tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTcg          >  1:2560801/1‑69 (MQ=255)
                                                  tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTcg          >  1:2114627/1‑69 (MQ=255)
                                                  tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTcg          >  1:3053387/1‑69 (MQ=255)
                                                       tttAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCgg       >  1:1245528/1‑67 (MQ=255)
                                                         tAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGACAAACTCTTCGATATAGTGTTCGCGGTCga   >  1:3274440/1‑69 (MQ=255)
                                                         tAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGACAAACTCTTCGATATAGTGTTCGCGGTCga   >  1:2473459/1‑69 (MQ=255)
                                                         tAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGACAAACTCTTCGATATAGTGTTCGCGGTCga   >  1:1307668/1‑69 (MQ=255)
                                                         tAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGACAAACTCTTCGATATAGTGTTCGCGGTCga   >  1:798432/1‑69 (MQ=255)
                                                         tAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGACAAACTCTTCGATATAGTGTTCGCGGTCga   >  1:2774766/1‑69 (MQ=255)
                                                         tAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGACAAACTCTTCGATATAGTGTTCGCGGTCga   >  1:2800208/1‑69 (MQ=255)
                                                         tAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGACAAACTCTTCGATATAGTGTTCGCGGTCga   >  1:473241/1‑69 (MQ=255)
                                                         tAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGACAAACTCTTCGATATAGTGTTCGCGGTCga   >  1:300113/1‑69 (MQ=255)
                                                         tAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGACAAACTCTTCGATATAGTGTTCGCGGTCga   >  1:1810400/1‑69 (MQ=255)
                                                         tAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGACAAACTCTTCGATATAGTGTTCGCGGTCga   >  1:3186839/1‑69 (MQ=255)
                                                         tAACACCTTCAGCTTTCATCTCAGCGCGCAGATGGCCGACAAACTCTTCGATATAGTGTTCGCGGTCg    >  1:3069598/1‑68 (MQ=255)
                                                           aCGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCgat  <  1:2170942/68‑2 (MQ=255)
                                                           aCGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCgag  <  1:1600824/68‑1 (MQ=255)
                                                           aCGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCgag  <  1:2061310/68‑1 (MQ=255)
                                                           aCGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCgag  <  1:2047790/68‑1 (MQ=255)
                                                           aCGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCgag  <  1:817783/68‑1 (MQ=255)
                                                           aCGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCgag  <  1:2116478/68‑1 (MQ=255)
                                                           aCGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCgag  <  1:945470/68‑1 (MQ=255)
                                                           aCGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCgag  <  1:1668544/68‑1 (MQ=255)
                                                             |                                                                 
TCTGCATTTTACGCCAGATGCTGTAGATGTGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAG  >  minE/1882341‑1882467

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: