Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 471,384 C→G 57.1% A250A (GCC→GCG 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,3840CG57.1% 15.1 / 49.9 49A250A (GCC→GCGsucDsuccinyl‑CoA synthetase, NAD(P)‑binding, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (15/6);  new base G (8/20);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 4.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

CCAGTTGTGGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTG  >  minE/471319‑471447
                                                                 |                                                               
ccAGTTGTGGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCcatc                                                              <  1:2645137/69‑1 (MQ=255)
      gtgGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTTCCATCATTGcc                                                        >  1:804620/1‑69 (MQ=255)
      gtgGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGc                                                         >  1:528124/1‑68 (MQ=255)
             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:406810/61‑1 (MQ=255)
                         tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                             <  1:2860504/61‑1 (MQ=255)
                         tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                             <  1:2562016/61‑1 (MQ=255)
                         tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                             <  1:2516445/61‑1 (MQ=255)
                         tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                             <  1:639342/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:1632753/61‑1 (MQ=255)
                         tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                             <  1:3212442/61‑1 (MQ=255)
                         tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                             <  1:3295216/61‑1 (MQ=255)
                         tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                             <  1:1504709/61‑1 (MQ=255)
                         tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                             <  1:1410580/61‑1 (MQ=255)
                         tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                                             <  1:1170156/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:3180582/1‑57 (MQ=255)
                             tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                             >  1:3005435/1‑57 (MQ=255)
                             tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                             >  1:3305094/1‑57 (MQ=255)
                             tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                             >  1:3408655/1‑57 (MQ=255)
                             tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                             >  1:2282614/1‑57 (MQ=255)
                             tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                             >  1:928445/1‑57 (MQ=255)
                             tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                             >  1:1526112/1‑57 (MQ=255)
                             tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                                             >  1:956448/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)
                                                                 |                                                               
CCAGTTGTGGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCAGGCGTG  >  minE/471319‑471447

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: