Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,514,504 G→A 24.5% R329R (CGC→CGT purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,514,5040GA24.5% 107.3 / 27.8 85R329R (CGC→CGTpurFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (43/21);  new base A (13/8);  total (56/29)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

CAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATT  >  minE/1514438‑1514555
                                                                  |                                                   
cAGTTTGCGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCg                                                   >  1:439275/1‑69 (MQ=255)
cAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAGGGTGCGGCCAACATAGCg                                                   >  1:795112/1‑69 (MQ=255)
cAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCg                                                   >  1:1005773/1‑69 (MQ=255)
cAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCg                                                   >  1:3050994/1‑69 (MQ=255)
cAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCg                                                   >  1:2981604/1‑69 (MQ=255)
cAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCg                                                   >  1:3405561/1‑69 (MQ=255)
cAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCg                                                   >  1:37423/1‑69 (MQ=255)
cAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCg                                                   >  1:2058423/1‑69 (MQ=255)
cAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCg                                                   >  1:473939/1‑69 (MQ=255)
cAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCg                                                   >  1:513813/1‑69 (MQ=255)
cAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCg                                                   >  1:723794/1‑69 (MQ=255)
cAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCg                                                   >  1:849744/1‑69 (MQ=255)
  gTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCgg                                                  <  1:3289844/68‑1 (MQ=255)
     tACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGttt                                               >  1:1085613/1‑68 (MQ=255)
     tACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGttt                                               >  1:2910099/1‑68 (MQ=255)
     tACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGttt                                               >  1:747487/1‑68 (MQ=255)
     tACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGttt                                               >  1:2696448/1‑68 (MQ=255)
     tACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGttt                                               >  1:2368440/1‑68 (MQ=255)
     tACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGttt                                               >  1:1153366/1‑68 (MQ=255)
              cGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGc                                     <  1:3084279/69‑1 (MQ=255)
              cGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGc                                     <  1:215394/69‑1 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAg                                      >  1:2363968/1‑68 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGc                                     >  1:2299822/1‑69 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGc                                     >  1:2152108/1‑69 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGc                                     >  1:2383278/1‑69 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGc                                     >  1:1122243/1‑69 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGc                                     >  1:571533/1‑69 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGc                                     >  1:2554709/1‑69 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGc                                     >  1:1873110/1‑69 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGc                                     >  1:703878/1‑69 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGc                                     >  1:1783858/1‑69 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGc                                     >  1:258230/1‑69 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGc                                     >  1:2610367/1‑69 (MQ=255)
              cGGACTTACGACGCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAAAATAACGGTTTTTAACGAAGc                                     >  1:1194744/1‑69 (MQ=255)
               ggACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGc                                     <  1:2445857/68‑1 (MQ=255)
                       gacgCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGCCCTGGCgg                             <  1:39067/68‑1 (MQ=255)
                       gacgCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGCCCTGGCgg                             <  1:2018453/68‑1 (MQ=255)
                       gacgCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGCCCTGGCgg                             <  1:1158219/68‑1 (MQ=255)
                       gacgCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGCCCTGGCgg                             <  1:3113635/68‑1 (MQ=255)
                       gacgCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGCCCTGGCgg                             <  1:1568739/68‑1 (MQ=255)
                       gacgCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGCCCTGGCgg                             <  1:758319/68‑1 (MQ=255)
                       gacgCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGCCCTGGCgg                             <  1:1629795/68‑1 (MQ=255)
                       gacgCAGCTGCTGACCCGGCATGATGAAGGTGCGGCCAACATAACGGTTTTTAACGAAGCCCTGGCgg                             <  1:2879268/68‑1 (MQ=255)
                                      ccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAg               <  1:1122308/67‑1 (MQ=255)
                                      ccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAg               <  1:1138310/67‑1 (MQ=255)
                                      ccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAg               <  1:862958/67‑1 (MQ=255)
                                      ccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAg               <  1:1709569/67‑1 (MQ=255)
                                      ccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGaa             <  1:1166213/69‑1 (MQ=255)
                                      ccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGaa             <  1:1439484/69‑1 (MQ=255)
                                      ccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGaa             <  1:3302288/69‑1 (MQ=255)
                                      ccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGaa             <  1:1736031/69‑1 (MQ=255)
                                      ccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGaa             <  1:2236930/69‑1 (MQ=255)
                                      ccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGaa             <  1:607528/69‑1 (MQ=255)
                                      ccGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGaa             <  1:2206109/69‑1 (MQ=255)
                                        tgCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATa           <  1:1076379/68‑1 (MQ=255)
                                        ggCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATa           <  1:1949757/69‑1 (MQ=255)
                                        ggCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATa           <  1:501968/69‑1 (MQ=255)
                                        ggCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATa           <  1:1958111/69‑1 (MQ=255)
                                        ggCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATa           <  1:2423907/69‑1 (MQ=255)
                                           aTGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATAc          >  1:3243563/1‑67 (MQ=255)
                                              aTAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCt                                  >  1:650472/1‑40 (MQ=255)
                                              aTAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCa     >  1:3400126/1‑69 (MQ=255)
                                              aTAAAGGTCCCGCCAACATAGCGGTTTTTAACGAAGCCCTcgcg                              >  1:3423979/1‑44 (MQ=255)
                                                   ggTGCGGCCAACATAGGGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:1920710/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:1789511/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:975374/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:1666037/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:697490/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:1261429/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:1241209/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:1097652/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:894585/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:1945089/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:2369864/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:3330652/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:3305222/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:3299/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:3135424/1‑67 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAt   >  1:241865/1‑66 (MQ=255)
                                                   ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAt   >  1:1615783/1‑66 (MQ=255)
                                                   ggTGCGGCCAACATAGAGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAAtt  >  1:1770352/1‑67 (MQ=255)
                                                                tAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTg                    >  1:2150437/1‑36 (MQ=255)
                                                                tAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTg                    >  1:2964689/1‑36 (MQ=255)
                                                                tAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTg                    >  1:2990184/1‑36 (MQ=255)
                                                                 aGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATa           <  1:2774382/44‑1 (MQ=255)
                                                                  |                                                   
CAGTTTACGGCGCACGGACTTACGACGCAGCTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATT  >  minE/1514438‑1514555

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: