Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,518,676 G→A 61.0% I133I (ATC→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,6760GA61.0% 47.4 / 68.5 100I133I (ATC→ATTspoTbifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (15/24);  new base A (42/19);  total (57/43)
Fisher's exact test for biased strand distribution p-value = 3.73e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCT  >  minE/2518617‑2518731
                                                           |                                                       
tCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGAcgcg                                                 <  1:2373070/68‑1 (MQ=255)
tCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGAcgcg                                                 <  1:1101275/68‑1 (MQ=255)
tCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGAcgcg                                                 <  1:2764175/68‑1 (MQ=255)
tCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGAcgcg                                                 <  1:3055391/68‑1 (MQ=255)
tCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGAcgcg                                                 <  1:1187393/68‑1 (MQ=255)
tCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGAcgcg                                                 <  1:1367708/68‑1 (MQ=255)
tCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGAcgcg                                                 <  1:3297651/68‑1 (MQ=255)
tCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGAcgcg                                                 <  1:3196839/68‑1 (MQ=255)
tCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGAcgcg                                                 <  1:695060/68‑1 (MQ=255)
tCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGAcgag                                                 <  1:2692290/68‑3 (MQ=255)
    ggCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGatat                                            <  1:203580/69‑1 (MQ=255)
    ggCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGatat                                            <  1:168804/69‑1 (MQ=255)
    ggCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGatat                                            <  1:1740086/69‑1 (MQ=255)
    ggCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGatat                                            <  1:2462189/69‑1 (MQ=255)
    ggCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGatat                                            <  1:1906874/69‑1 (MQ=255)
      cGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGAcgcg                                                 <  1:1668234/62‑1 (MQ=255)
          gTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGc                                       >  1:2247645/1‑68 (MQ=255)
          gTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGc                                       >  1:2283252/1‑68 (MQ=255)
                 cAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGAcgcg                                                 <  1:863327/51‑1 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGTACCATCGc                               >  1:439670/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:2232682/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:3033331/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:2458462/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:3019822/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:2417721/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:2421270/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:244544/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:1053258/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:3004185/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:2504659/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:2507697/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:2514555/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:254077/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:2732939/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:2846550/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:1418332/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:1161824/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:866043/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:1210978/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:772727/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:707905/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:646798/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:1466618/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:425108/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:394322/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:1685362/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:1700373/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:3169098/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGc                               >  1:1102679/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGc                               >  1:2632868/1‑68 (MQ=255)
                  aGCGTGCGCATGTTGTGGGTACGGTCGGAAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGc                               >  1:38319/1‑68 (MQ=255)
                        cgcATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGa                                                     >  1:1250923/1‑40 (MQ=255)
                        cgcATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatc                        >  1:1389888/1‑69 (MQ=255)
                        cgcATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatc                        >  1:598563/1‑69 (MQ=255)
                        cgcATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatc                        >  1:2476495/1‑69 (MQ=255)
                        cgcATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatc                        >  1:207257/1‑69 (MQ=255)
                        cgcATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatc                        >  1:2001837/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:1064008/1‑69 (MQ=255)
                          cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatcat                      >  1:1270565/1‑69 (MQ=255)
                          cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAatcat                      >  1:1275770/1‑69 (MQ=255)
                          cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGCCGCGGATATCCTg                                        >  1:126152/1‑51 (MQ=255)
                          cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                      >  1:1506263/1‑69 (MQ=255)
                          cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                      >  1:959450/1‑69 (MQ=255)
                          cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                      >  1:870608/1‑69 (MQ=255)
                          cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                      >  1:110304/1‑69 (MQ=255)
                          cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                      >  1:2572002/1‑69 (MQ=255)
                          cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                      >  1:1855609/1‑69 (MQ=255)
                          cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                      >  1:2904489/1‑69 (MQ=255)
                          cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                      >  1:2302667/1‑69 (MQ=255)
                          cATGTTGTGGGTACGGTCGGCAAGTTTGATGAGAATGACGCGGATATCCTGCACCATCGCCATAatcat                      >  1:1966294/1‑69 (MQ=255)
                              ttgtGGGTACGGTCGGCAAGTTTGATGAGGATGAcgcg                                                 <  1:1521253/38‑1 (MQ=255)
                                         gTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAat                         <  1:3410214/51‑1 (MQ=255)
                                              cAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGcc   <  1:1887775/68‑1 (MQ=255)
                                              cAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGcc   <  1:728649/68‑1 (MQ=255)
                                              cAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGcc   <  1:1474542/68‑1 (MQ=255)
                                              cAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGcc   <  1:2037040/68‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCTAAAGTTTTCGGCCt  <  1:1327588/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCTAAAGTTTTCGGCCt  <  1:118749/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:1881032/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:3104263/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:937937/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:2983283/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:246102/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:856235/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:3163006/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:379258/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:2183165/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:575793/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:558700/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:3175448/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:410394/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:409305/69‑1 (MQ=255)
                                              cAAGTTTGATGAGAATGACGCGGATATCCTGCACAATCGCCATAATCATCTTGCGAAAGTTTTCGGCCt  <  1:1755182/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)
                                                           |                                                       
TCCGGGCGAAGTGAGCCCAGCGTGCGCATGTTGTGGGTACGGTCGGCAAGTTTGATGAGGATGACGCGGATATCCTGCACCATCGCCATAATCATCTTGCGAAAGTTTTCGGCCT  >  minE/2518617‑2518731

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: