Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 471,372 C→T 56.0% G246G (GGC→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,3720CT56.0% 6.3 / 37.8 50G246G (GGC→GGTsucDsuccinyl‑CoA synthetase, NAD(P)‑binding, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (14/8);  new base T (8/20);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 2.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

AAGCCAGTTGTGGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTG  >  minE/471316‑471427
                                                        |                                                       
aaGCCAGTTGTGGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGc                                              <  1:1088038/68‑1 (MQ=255)
aaGCCAGTTGTGGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGc                                              <  1:2418232/68‑1 (MQ=255)
aaGCCAGTTGTGGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGc                                              <  1:2069439/68‑1 (MQ=255)
aaGCCAGTTGTGGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGc                                              <  1:2419921/68‑1 (MQ=255)
   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:3295216/61‑1 (MQ=255)
                            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                         <  1:3212442/61‑1 (MQ=255)
                            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                         <  1:1116057/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:406810/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:1170156/61‑1 (MQ=255)
                            tgACTGCGCCGAAAGGCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAgg                         <  1:1800720/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:2282614/1‑57 (MQ=255)
                                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                         >  1:3408655/1‑57 (MQ=255)
                                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                         >  1:3305094/1‑57 (MQ=255)
                                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                         >  1:1526112/1‑57 (MQ=255)
                                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                         >  1:3180582/1‑57 (MQ=255)
                                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                         >  1:3005435/1‑57 (MQ=255)
                                tGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAgg                         >  1:928445/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:1161576/1‑69 (MQ=255)
                                           gCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAGGGACTGCGGATGAGAAATTCgctg  >  1:968694/1‑69 (MQ=255)
                                           gCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAGGGACTGCGGATGAGAAATTCgctg  >  1:1327608/1‑69 (MQ=255)
                                           gCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAGGGACTGCGGATGAGAAATTCgctg  >  1:372042/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)
                                           gCAAACGTATGGGTCACGCAGGTGCGATCATTGCCGGTGGTAAAGGGACTGCGGATGAGAAATTCgct   >  1:2999487/1‑68 (MQ=255)
                                                        |                                                       
AAGCCAGTTGTGGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTG  >  minE/471316‑471427

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: