Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,109,347 C→T 20.3% L88F (CTC→TTC)  ydiS → predicted oxidoreductase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,109,3470CT20.3% 129.5 / 26.6 59L88F (CTC→TTC) ydiSpredicted oxidoreductase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (25/22);  new base T (8/4);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

TCAGCGCCGGTAGAACGCAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCT  >  minE/1109278‑1109414
                                                                     |                                                                   
tCAGCGCCGGTAGAACGCAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCt                                                                    >  1:2348715/1‑71 (MQ=255)
tCAGCGCCGGTAGAACGCAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCt                                                                    <  1:3375275/71‑1 (MQ=255)
      ccGGTAGAACGCAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGAttt                                                              <  1:1434793/71‑1 (MQ=255)
        ggTAGAACGCAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGAtttt                                                             >  1:2061463/1‑70 (MQ=255)
        ggTAGAACGCAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGAtttt                                                             >  1:3706638/1‑70 (MQ=255)
           aGAACGCAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCAcc                                                         >  1:2537312/1‑71 (MQ=255)
              aCGCAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCt                                                                    >  1:2625885/1‑57 (MQ=255)
              aCGCAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACcgcg                                                      <  1:1192394/71‑1 (MQ=255)
               cGCAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACcgc                                                       >  1:3245355/1‑69 (MQ=255)
                gCAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGATACCCTGGATTTTCACCGCGag                                                    >  1:1715726/1‑71 (MQ=255)
                 cAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGag                                                    >  1:2645358/1‑70 (MQ=255)
                   aGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGagcag                                                 <  1:2480813/71‑1 (MQ=255)
                    ggTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGagcagc                                                >  1:4098260/1‑71 (MQ=255)
                     gTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGcc                                               <  1:2760034/71‑1 (MQ=255)
                      tCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGAtttt                                                             >  1:2390123/1‑56 (MQ=255)
                      tCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGagcagc                                                >  1:1605134/1‑69 (MQ=255)
                      tCACACCCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCa                                              <  1:1313832/71‑1 (MQ=255)
                       cacaCGCGAGAAATTCTCCTTCTTACCCGAAGATAGCGCCGTTACCTTCGATTTTCACCGCGagcagc                                                >  1:2180143/1‑68 (MQ=255)
                       cacaCGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGagcagc                                                >  1:2956266/1‑68 (MQ=255)
                               agaAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCTTCGATTTTCACCGCGAGCAGCCAGAt                                           <  1:2629311/65‑1 (MQ=255)
                               agaAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCTTCGATTTTCACCGCGAGCAGCCAGAt                                           <  1:1216003/65‑1 (MQ=255)
                               agaAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCg                                     >  1:681932/1‑71 (MQ=255)
                                gaAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCTTCGATTTTCACCGCGAGCAGCCAGATGTTCCg                                     >  1:2708561/1‑70 (MQ=255)
                                gaAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGc                                    <  1:2110869/71‑1 (MQ=255)
                                 aaaaTCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGAt                                           <  1:2274911/63‑1 (MQ=255)
                                 aaaaTCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGAt                                           <  1:3231397/63‑1 (MQ=255)
                                 aaaaTCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTc                                       <  1:1013882/67‑1 (MQ=255)
                                  aaaTCTCCTTCTTAACCGAAGAAAGCGCCGTTAGCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCaa                                  <  1:460339/71‑1 (MQ=255)
                                  aaaTCTCCTTCTTAACCGAAGAAAGCGCCGTTACCTTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCaa                                  <  1:3274973/71‑1 (MQ=255)
                                   aaTCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGa                                            >  1:3988180/1‑60 (MQ=255)
                                     tctcCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAacac                               >  1:1443019/1‑71 (MQ=255)
                                       tcCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAacac                               >  1:2215875/1‑69 (MQ=255)
                                         cttcttAACCGAAGAAAGCGCCGTTACCCTCGATTTTc                                                            >  1:1302458/1‑38 (MQ=255)
                                          ttcttAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACcgc                                                       <  1:384522/42‑1 (MQ=255)
                                           tcttAACCGAAGAAAGCGCCGTTACCTTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCaa                                  >  1:604256/1‑62 (MQ=255)
                                           tcttAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCaa                                  <  1:1993678/62‑1 (MQ=255)
                                            cttAACCGAAGAAAGCGCCGTTACCTTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATc                          <  1:2897168/69‑1 (MQ=255)
                                            cttAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCt                         <  1:1012366/70‑1 (MQ=255)
                                             ttAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCAt                           <  1:2521455/67‑1 (MQ=255)
                                              tAACCGAAGAAAGCGCCGTTACCTTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTta                       >  1:552584/1‑70 (MQ=255)
                                              tAACCGAAGAAAGCGCCGTTACCTTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTta                       >  1:699465/1‑70 (MQ=255)
                                               aaCCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCt                         >  1:567321/1‑67 (MQ=255)
                                               aaCCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCt                         >  1:2351605/1‑67 (MQ=255)
                                               aaCCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTta                       <  1:1739391/69‑1 (MQ=255)
                                                aCCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAg                                             <  1:2213869/46‑1 (MQ=255)
                                                   gaagaaAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAg                                             <  1:285656/43‑1 (MQ=255)
                                                   gaagaaAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGc                                    <  1:2528898/52‑1 (MQ=255)
                                                    aagaaAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGa                                            >  1:1220033/1‑43 (MQ=255)
                                                      gaaAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTtat                      <  1:1969776/63‑1 (MQ=255)
                                                        aaGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCaa                                  <  1:3881534/49‑1 (MQ=255)
                                                         aGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCg           >  1:2872269/1‑71 (MQ=255)
                                                         aGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCg           >  1:3886572/1‑71 (MQ=255)
                                                         aGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGAAGTTCCGCAACACGCATCTTATACCGTATTGCg           >  1:3318022/1‑71 (MQ=255)
                                                             ccGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAac                                 <  1:4054945/45‑1 (MQ=255)
                                                               gTTACCTTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCg     >  1:1324841/1‑71 (MQ=255)
                                                               gTTACCTTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATAACGTATTGCGTAATCg     >  1:2891636/1‑71 (MQ=255)
                                                                ttACCTTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACg                              >  1:3383972/1‑45 (MQ=255)
                                                                 tACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTc   >  1:717436/1‑71 (MQ=255)
                                                                  aCCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCt  >  1:831354/1‑71 (MQ=255)
                                                                    ccTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCAt                           <  1:4107891/44‑1 (MQ=255)
                                                                    ccTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCAt                           <  1:1488796/44‑1 (MQ=255)
                                                                    ccTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCAt                           <  1:2074388/44‑1 (MQ=255)
                                                                    ccTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCAt                           <  1:954482/44‑1 (MQ=255)
                                                                     |                                                                   
TCAGCGCCGGTAGAACGCAAGGTCACACGCGAGAAAATCTCCTTCTTAACCGAAGAAAGCGCCGTTACCCTCGATTTTCACCGCGAGCAGCCAGATGTTCCGCAACACGCATCTTATACCGTATTGCGTAATCGTCT  >  minE/1109278‑1109414

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: