Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 471,399 G→T 53.8% G255G (GGG→GGT sucD → succinyl‑CoA synthetase, NAD(P)‑binding, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE471,3990GT53.8% ‑2.5 / 57.8 52G255G (GGG→GGTsucDsuccinyl‑CoA synthetase, NAD(P)‑binding, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (18/6);  new base T (8/20);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.91e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGTT  >  minE/471332‑471456
                                                                   |                                                         
aCATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGa                                                          <  1:2957823/69‑1 (MQ=255)
aCATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGa                                                          <  1:747768/69‑1 (MQ=255)
     gCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAggg                                                     <  1:1935498/69‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAACCGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:3352873/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:1116057/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:2860504/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:3212442/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:3295216/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:2562016/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:2516445/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:2304636/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:2170954/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:1986798/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:1800720/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:406810/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:639342/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:1170156/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:1632753/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:1410580/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:1504709/61‑1 (MQ=255)
            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGACGGTGGTAAAgg                                                      <  1:300281/61‑1 (MQ=255)
            tgACTGCCCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:1624947/61‑1 (MQ=255)
              aCTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCgg                                              >  1:1720646/1‑67 (MQ=255)
              aCTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGAATGCgg                                              >  1:2989167/1‑67 (MQ=255)
                tGCGCCGAAAGGCAAACGTATGGGCCCCGCGGGTGCCATCATTGCCGGTGGGAAAgg                                                      >  1:1253324/1‑57 (MQ=255)
                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                                      >  1:3408655/1‑57 (MQ=255)
                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                                      >  1:928445/1‑57 (MQ=255)
                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                                      >  1:2282614/1‑57 (MQ=255)
                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                                      >  1:1526112/1‑57 (MQ=255)
                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                                      >  1:956448/1‑57 (MQ=255)
                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                                      >  1:3305094/1‑57 (MQ=255)
                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                                      >  1:3180582/1‑57 (MQ=255)
                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                                      >  1:3005435/1‑57 (MQ=255)
                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCAACATTGCCGGTGGGAAAgg                                                      >  1:1185660/1‑57 (MQ=255)
                       aaaGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                                      <  1:1880251/50‑1 (MQ=255)
                           gCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAGGGACTGCGGATGAGAAATTCgctg                               >  1:372042/1‑69 (MQ=255)
                           gCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAGGGACTGCGGATGAGAAATTCgctg                               >  1:968694/1‑69 (MQ=255)
                           gCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAGGGACTGCGGATGAGAAATTCgctg                               >  1:2989555/1‑69 (MQ=255)
                           gCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAGGGACTGCGGATGAGAAATTCgctg                               >  1:2754468/1‑69 (MQ=255)
                           gCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAGGGACTGCGGATGAGAAATTCgctg                               >  1:2305955/1‑69 (MQ=255)
                           gCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAGGGACTGCGGATGAGAAATTCgctg                               >  1:1327608/1‑69 (MQ=255)
                           gCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAGGGACTGCGGATGAGAAATTCgctg                               >  1:1161576/1‑69 (MQ=255)
                           gCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAGGGACTGCGGATGAGAAATTCgct                                >  1:2999487/1‑68 (MQ=255)
                                         cACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCa                 >  1:599261/1‑69 (MQ=255)
                                               ggTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTg           <  1:2093062/69‑1 (MQ=255)
                                               ggTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTg           <  1:2089071/69‑1 (MQ=255)
                                               ggTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTg           <  1:1577045/69‑1 (MQ=255)
                                                        atTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGtt  >  1:2735299/1‑69 (MQ=255)
                                                        atTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGtt  >  1:2572122/1‑69 (MQ=255)
                                                        atTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGtt  >  1:1764958/1‑69 (MQ=255)
                                                        atTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGtt  >  1:1733561/1‑69 (MQ=255)
                                                        atTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGtt  >  1:1677908/1‑69 (MQ=255)
                                                                   |                                                         
ACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTGAAAACCGTT  >  minE/471332‑471456

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: