Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,480,246 C→T 40.5% G16G (GGC→GGT glpA → sn‑glycerol‑3‑phosphate dehydrogenase (anaerobic), large subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,480,2460CT40.5% 29.4 / 39.8 64G16G (GGC→GGTglpAsn‑glycerol‑3‑phosphate dehydrogenase (anaerobic), large subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (28/10);  new base T (16/10);  total (44/20)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAACTTCAGAGGGATAACAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGCCCGCGACTGTGCCCTGCGCGGGCTGCGCGTGATTTTGGTTGAGC  >  minE/1480179‑1480313
                                                                   |                                                                   
cGAACTTCAGAGGGATAACAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTggcggct                                                                    >  1:1587845/1‑68 (MQ=255)
cGAACTTCAGAGGGATAACAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTggcggcg                                                                    >  1:2838512/1‑69 (MQ=255)
cGAACTTCAGAGGGATAACAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTggcggcg                                                                    >  1:2509464/1‑69 (MQ=255)
cGAACTTCAGAGGGATAACAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTggcggcg                                                                    >  1:1770474/1‑69 (MQ=255)
cGAACTTCAGAGGGATAACAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTggcggcg                                                                    >  1:2273823/1‑69 (MQ=255)
cGAACTTCAGAGGGATAACAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTggcggcg                                                                    >  1:702018/1‑69 (MQ=255)
            ggATAACAATGAAAACTCGCGACTCTCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAg                                                        >  1:3355764/1‑69 (MQ=255)
            ggATAACAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAg                                                        >  1:895027/1‑69 (MQ=255)
            ggATAACAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAg                                                        >  1:260176/1‑69 (MQ=255)
                 aCAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCggg                                                   >  1:1431792/1‑69 (MQ=255)
                 aCAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCgg                                                    >  1:1773981/1‑68 (MQ=255)
                 aCAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCgg                                                    >  1:752899/1‑68 (MQ=255)
                 aCAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCgg                                                    >  1:2269627/1‑68 (MQ=255)
                 aCAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCg                                                     >  1:2671398/1‑67 (MQ=255)
                 aCAATCAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCgg                                                    <  1:2544358/68‑1 (MQ=255)
                 aCAATCAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCgg                                                    <  1:1153049/68‑1 (MQ=255)
                 aCAATCAAAACCCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCgg                                                    <  1:3263568/68‑1 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCggg                                                   >  1:2876083/1‑68 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:277710/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:1819857/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:747877/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:2090430/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:3153620/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:1635610/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:3222703/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:1415459/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:1392281/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:2571564/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:756595/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:3335924/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:2390224/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:2484132/1‑69 (MQ=255)
                  cAATGAAAACTCGCGACTCGCAATCAAGAGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGa                                                  >  1:484015/1‑69 (MQ=255)
                   aaTGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGTGAGCCGGGa                                                  >  1:671671/1‑68 (MQ=255)
                   aaTGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGAt                                                 >  1:2881839/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             >  1:3328957/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             >  1:1806525/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             >  1:2803119/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             >  1:1810921/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             >  1:3018417/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             >  1:1558514/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             >  1:2644857/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             >  1:3061738/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             >  1:2296977/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             >  1:1046944/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             >  1:2348412/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             >  1:3150354/1‑69 (MQ=255)
                       aaaaCTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGc                                              >  1:2943795/1‑68 (MQ=255)
                              gcgACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGCCCGCGAc                                       <  1:1770821/68‑1 (MQ=255)
                                gACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCgg                                                    <  1:847725/53‑1 (MQ=255)
                                gACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGcc                                             <  1:244818/60‑1 (MQ=255)
                                gACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGACgg                                                    <  1:2731633/53‑1 (MQ=255)
                                 aCTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCgg                                                    <  1:2177387/52‑1 (MQ=255)
                                         tCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGATTGCCCGCGACtg                                     <  1:1559379/59‑1 (MQ=255)
                                         tCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGATTGCCCGCGACtg                                     <  1:2492333/59‑1 (MQ=255)
                                         tCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGATTGCCCGCGACtg                                     <  1:1662472/59‑1 (MQ=255)
                                         tCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGATTGCCCGCGACtg                                     <  1:2755077/59‑1 (MQ=255)
                                         tCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGATTGCCCGCGACtg                                     <  1:2674552/59‑1 (MQ=255)
                                         tCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGATTGCCCGCGACtg                                     <  1:2152903/59‑1 (MQ=255)
                                         tCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGATTGCCCGCGACtg                                     <  1:2347488/59‑1 (MQ=255)
                                         tCAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGATTGCCCGCGACtg                                     <  1:59892/59‑1 (MQ=255)
                                          cAAGTGACGTGATTATCATTGGCGGTGGCGCAACGGGAGCCGGGATTGCCCGCGACtg                                     <  1:1093270/58‑1 (MQ=255)
                                                    gATTATCATTGGCGGCGGCGCAACGTGAGCCGGGATAGCCCGCGACtgtg                                   >  1:2117212/1‑50 (MQ=39)
                                                    gATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGCCCGCGACTGTGCCCTGCGCGGGCTGCGCGt                >  1:557883/1‑69 (MQ=255)
                                                    gATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGCCCGCGACTGTGCCCTGCGCGGGCTGCGCGt                >  1:733874/1‑69 (MQ=255)
                                                    gATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGCCCGCGACTGTGCCCTGCGCGGGCTGCGCGt                >  1:1195372/1‑69 (MQ=255)
                                                     aTTATCATTGGCGGTGGCGCAACGGGAGCCGGGATTGCCCGCGACtg                                     <  1:1261886/47‑1 (MQ=255)
                                                     aTTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGCCCGCGACTGTGCCCTGCGCGGGCTGCGCGTg               >  1:863736/1‑69 (MQ=255)
                                                           aTTGGCGGCGGCGCAACGGGAGCCGGGATTGCCCGCGACTGTGCCCTGCGCGGGCTGCGCGTGATTTTg         <  1:2041198/69‑1 (MQ=255)
                                                            ttGGCGGCGGCGCAACGGGAGCCGGGATTGCCCGCGACTGTGCCCTGCGCGGGCTGCGCGTGATTTTg         <  1:2520318/68‑1 (MQ=255)
                                                                  gcggcgCAACGGGAGCCGGGATTGCCCGCGACTGTGCCCTGCGTGGGCTGCGCGTGATTTTGGTTGAgc  <  1:1445026/69‑1 (MQ=255)
                                                                  gcggcgCAACGGGAGCCGGGATTGCCCGCGACTGTGCCCTGCGCGGGCTGCGCGTGATTTTGGTTGAgc  <  1:417972/69‑1 (MQ=255)
                                                                  gcggcgCAACGGGAGCCGGGATTGCCCGCGACTGTGCCCTGCGCGGGCTGCGCGTGATTTTGGTTGAgc  <  1:2543925/69‑1 (MQ=255)
                                                                  gcggcgCAACGGGAGCCGGGATTGCCCGCGACTGTGCCCTGCGCGGGCTGCGCGTGATTTTGGTTGAgc  <  1:117274/69‑1 (MQ=255)
                                                                   |                                                                   
CGAACTTCAGAGGGATAACAATGAAAACTCGCGACTCGCAATCAAGTGACGTGATTATCATTGGCGGCGGCGCAACGGGAGCCGGGATTGCCCGCGACTGTGCCCTGCGCGGGCTGCGCGTGATTTTGGTTGAGC  >  minE/1480179‑1480313

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: