Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 225,827 T→C 8.2% intergenic (+14/‑78) gpt → / → frsA guanine‑hypoxanthine phosphoribosyltransferase/hydrolase, binds to enzyme IIA(Glc)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE225,8270TC 8.2% 133.8 / 4.4 49intergenic (+14/‑78)gpt/frsAguanine‑hypoxanthine phosphoribosyltransferase/hydrolase, binds to enzyme IIA(Glc)
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base C (2/2);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

GAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTG  >  minE/225760‑225892
                                                                   |                                                                 
gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTgg                                                                 <  1:1765466/70‑1 (MQ=255)
gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTgg                                                                 <  1:5028230/70‑1 (MQ=255)
 aaCAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGc                                                                <  1:1201328/70‑1 (MQ=255)
 aaCAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGc                                                                <  1:4798797/70‑1 (MQ=255)
 aaCAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCgg                                                                 >  1:348639/1‑69 (MQ=255)
 aaCAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCgg                                                                 >  1:1643718/1‑69 (MQ=255)
      ccGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGcc                                                          <  1:6727811/71‑1 (MQ=255)
      ccGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGcc                                                          >  1:3704418/1‑71 (MQ=255)
      ccGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGcc                                                          <  1:5539184/71‑1 (MQ=255)
        gTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGccag                                                        <  1:5762608/71‑3 (MQ=255)
        gTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGccag                                                        <  1:1477079/71‑3 (MQ=255)
         tGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCggg                                                       <  1:2042658/71‑1 (MQ=255)
                  ggCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTc                                               <  1:4111681/70‑1 (MQ=255)
                   gCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                             <  1:3218264/71‑1 (MQ=255)
                        gtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGc                                                      >  1:1390820/1‑57 (MQ=255)
                        gtATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTaa                                        >  1:6914322/1‑71 (MQ=255)
                         tATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgttgg                                                 >  1:5593940/1‑60 (MQ=255)
                         tATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAAc                                       >  1:4656433/1‑71 (MQ=255)
                           ttCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttttt                                          <  1:1585936/66‑1 (MQ=255)
                              gTCCCGCCAATCTCCGGTCGCTTATCTTTTCAACGCCTg                                                                  >  1:4774387/1‑39 (MQ=255)
                              gTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTg                                                                  >  1:6398691/1‑39 (MQ=255)
                              gTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCa                                                               >  1:3459226/1‑42 (MQ=255)
                                 ccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCtt                                             >  1:3731137/1‑57 (MQ=255)
                                 ccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCg                               >  1:1543348/1‑71 (MQ=255)
                                       aTCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTAc                         <  1:5859689/71‑1 (MQ=255)
                                              gTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGtt                  <  1:869613/71‑1 (MQ=255)
                                               tCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTAAGGCGGGTTACAATAGttt                 >  1:2164411/1‑71 (MQ=255)
                                                  cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTa                                         <  1:907809/44‑1 (MQ=255)
                                                    aaTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACtt                                     >  1:3108631/1‑46 (MQ=255)
                                                    aaTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCag             >  1:6336373/1‑70 (MQ=255)
                                                    aaTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCa              >  1:674856/1‑69 (MQ=255)
                                                     aTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCag             >  1:6830335/1‑69 (MQ=255)
                                                      tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttttt                                          <  1:4025559/39‑1 (MQ=255)
                                                      tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttttt                                          <  1:2338840/39‑1 (MQ=255)
                                                      tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttttt                                          <  1:6338431/39‑1 (MQ=255)
                                                      tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttttt                                          <  1:559651/39‑1 (MQ=255)
                                                       cTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag         >  1:6640391/1‑71 (MQ=255)
                                                       cTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag         >  1:4128186/1‑71 (MQ=255)
                                                       cTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag         >  1:5334213/1‑71 (MQ=255)
                                                       cTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag         >  1:4906610/1‑71 (MQ=255)
                                                          ttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaagta       <  1:2383408/70‑1 (MQ=255)
                                                          ttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaagt        <  1:2932939/69‑1 (MQ=255)
                                                          ttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaa          <  1:537430/67‑1 (MQ=255)
                                                          ttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaa          <  1:1232537/67‑1 (MQ=255)
                                                          ttCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAt      >  1:2085448/1‑71 (MQ=255)
                                                           tCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTAt      <  1:6206191/70‑1 (MQ=255)
                                                             aaCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTc    <  1:5023337/70‑1 (MQ=255)
                                                             aaCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTc    <  1:3268487/70‑1 (MQ=255)
                                                              aCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg  <  1:3958153/71‑1 (MQ=255)
                                                                   |                                                                 
GAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTG  >  minE/225760‑225892

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: