Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 471,378 G→A 54.9% A248A (GCG→GCA 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,3780GA54.9% 4.7 / 39.8 51A248A (GCG→GCAsucDsuccinyl‑CoA synthetase, NAD(P)‑binding, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (15/8);  new base A (8/20);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCAGTTGTGGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCA  >  minE/471316‑471441
                                                              |                                                               
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)
                                                         cACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCa  >  1:599261/1‑69 (MQ=255)
                                                              |                                                               
AAGCCAGTTGTGGGTTACATCGCTGGTGTGACTGCGCCGAAAGGCAAACGTATGGGCCACGCGGGTGCCATCATTGCCGGTGGGAAAGGGACTGCGGATGAGAAATTCGCTGCTCTGGAAGCCGCA  >  minE/471316‑471441

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: