Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,518,697 C→A 26.3% M126I (ATG→ATT spoT ← bifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,518,6970CA26.3% 126.6 / 39.9 87M126I (ATG→ATTspoTbifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (58/6);  new base A (2/21);  total (60/27)
Fisher's exact test for biased strand distribution p-value = 8.23e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

AGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTT  >  minE/2518635‑2518762
                                                              |                                                                 
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGTACCATCGc                                                              >  1:439670/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:2458462/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:3004185/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:3019822/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:2232682/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:3033331/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:2846550/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:2417721/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:2421270/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:3169098/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:1700373/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:1685362/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:244544/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:394322/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:1466618/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:1418332/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:425108/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:2732939/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:2504659/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:2507697/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:646798/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:707905/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:1210978/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:772727/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:1161824/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:866043/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:1102679/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:2514555/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:1053258/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                                                              >  1:254077/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGc                                                              >  1:2632868/1‑68 (MQ=255)
aGCGTGCGCATGTTGTGGGTACGGTCGGAAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGc                                                              >  1:38319/1‑68 (MQ=255)
      cgcATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatc                                                       >  1:207257/1‑69 (MQ=255)
      cgcATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatc                                                       >  1:2001837/1‑69 (MQ=255)
      cgcATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatc                                                       >  1:1389888/1‑69 (MQ=255)
      cgcATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatc                                                       >  1:598563/1‑69 (MQ=255)
      cgcATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatc                                                       >  1:2476495/1‑69 (MQ=255)
      cgcATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAat                                                        >  1:3130066/1‑68 (MQ=255)
      cgaaTGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatc                                                       >  1:3077055/4‑69 (MQ=255)
      cgaaTGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatc                                                       >  1:3080988/4‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:2896553/1‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:1275770/1‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:1270565/1‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:1064008/1‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:2904489/1‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:870608/1‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:959450/1‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:2572002/1‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:1506263/1‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:1855609/1‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:1966294/1‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:2302667/1‑69 (MQ=255)
        cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                                                     >  1:110304/1‑69 (MQ=255)
                       gTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAat                                                        <  1:3410214/51‑1 (MQ=255)
                            cAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGcc                                  <  1:728649/68‑1 (MQ=255)
                            cAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGcc                                  <  1:2037040/68‑1 (MQ=255)
                            cAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGcc                                  <  1:1887775/68‑1 (MQ=255)
                            cAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGcc                                  <  1:1474542/68‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCTAAAGTTTTCGGCCt                                 <  1:118749/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCTAAAGTTTTCGGCCt                                 <  1:1327588/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:856235/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:575793/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:558700/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:937937/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:410394/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:409305/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:246102/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:3175448/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:1755182/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:3163006/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:2983283/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:2183165/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:379258/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:3104263/69‑1 (MQ=255)
                            cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:1881032/69‑1 (MQ=255)
                             aaGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGcc                                  <  1:1245433/67‑1 (MQ=255)
                             aaGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:129385/68‑1 (MQ=255)
                              aGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:2362925/67‑1 (MQ=255)
                                                    tatCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:2808026/45‑1 (MQ=255)
                                                      tCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACtt        >  1:3153164/1‑68 (MQ=255)
                                                      tCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTg       >  1:1409663/1‑69 (MQ=255)
                                                      tCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTg       >  1:2363969/1‑69 (MQ=255)
                                                      tCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTg       >  1:1027086/1‑69 (MQ=255)
                                                      tCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTg       >  1:1010710/1‑69 (MQ=255)
                                                      tCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt                                 <  1:742398/43‑1 (MQ=255)
                                                           cACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGttt  >  1:473469/1‑69 (MQ=255)
                                                           cACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGttt  >  1:1242627/1‑69 (MQ=255)
                                                              |                                                                 
AGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCTGCGCCTCTTTCTTATCGCGGAACTTGAGTTT  >  minE/2518635‑2518762

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: