Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,218,861 C→T 37.3% V32V (GTC→GTT yhdH → predicted oxidoreductase, Zn‑dependent and NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,218,8610CT37.3% 44.2 / 45.8 75V32V (GTC→GTTyhdHpredicted oxidoreductase, Zn‑dependent and NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (33/14);  new base T (15/13);  total (48/27)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATC  >  minE/2218800‑2218924
                                                             |                                                               
aaaCTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTc                                                          <  1:2745605/69‑1 (MQ=255)
   ctctCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCAc                                                        >  1:3277612/1‑68 (MQ=255)
   ctctCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACt                                                       >  1:88030/1‑69 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:3195312/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:2707262/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:3333393/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:3427410/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:2553197/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:2936667/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:2374458/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:2099487/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:1715226/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:1712728/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:3257297/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:1159794/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:1279881/1‑68 (MQ=255)
      tcGCATCAGTACAGACTCTGGACGAAAGTCGCATGCCGGAAGGCGATGTCACGGTTGATGTTCACTgg                                                     >  1:24138/1‑68 (MQ=255)
      tcGCATAAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAAGGCGATCTCACGGTTGATGTTCACTgg                                                     >  1:435752/1‑68 (MQ=255)
          aTCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGa                                                 <  1:1183312/68‑1 (MQ=255)
          aTCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGa                                                 <  1:110158/68‑1 (MQ=255)
          aTCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGa                                                 <  1:3018111/68‑1 (MQ=255)
          aTCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGa                                                 <  1:2048228/68‑1 (MQ=255)
          aTCAGTACAGACTCAGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGa                                                 <  1:1350478/68‑1 (MQ=255)
             aGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGcc                                              >  1:1080582/1‑68 (MQ=255)
             aGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGcc                                              >  1:327354/1‑68 (MQ=255)
               tACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTg                                            >  1:507033/1‑68 (MQ=255)
               tACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTg                                            <  1:307489/68‑1 (MQ=255)
               tACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTg                                            >  1:3270196/1‑68 (MQ=255)
               tACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGa                                           >  1:3259234/1‑69 (MQ=255)
                     ctctGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACtata                                     <  1:1311818/69‑1 (MQ=255)
                     ctctGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACtata                                     <  1:1380096/69‑1 (MQ=255)
                       ctGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATaaa                                   <  1:3388759/69‑1 (MQ=255)
                       ctGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATaaa                                   <  1:149885/69‑1 (MQ=255)
                       ctGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATaaa                                   <  1:1705254/69‑1 (MQ=255)
                        tGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTg                                            >  1:250394/1‑59 (MQ=255)
                        tGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTg                                            >  1:1256559/1‑59 (MQ=255)
                        tGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTg                                            >  1:2156729/1‑59 (MQ=255)
                        tGGACGAAAGTCGCCTGACGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTg                                            >  1:642602/1‑59 (MQ=255)
                                aGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCt                           >  1:928564/1‑68 (MQ=255)
                                aGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTg                          >  1:1033388/1‑69 (MQ=255)
                                aGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTg                          >  1:422805/1‑69 (MQ=255)
                                aGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTg                          >  1:193254/1‑69 (MQ=255)
                                aGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTg                          >  1:728583/1‑69 (MQ=255)
                                                     gTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGaaaa     >  1:2675324/1‑69 (MQ=255)
                                                     gTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGaaaa     >  1:3276489/1‑69 (MQ=255)
                                                     gTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGaaaa     >  1:2355515/1‑69 (MQ=255)
                                                     gTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGaaaa     >  1:446427/1‑69 (MQ=255)
                                                     gTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGaaaa     >  1:3220809/1‑69 (MQ=255)
                                                        aCGGTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:384042/69‑1 (MQ=255)
                                                        aCGGTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:557046/69‑1 (MQ=255)
                                                        aCGGTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:910860/69‑1 (MQ=255)
                                                        aCGGTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:1044751/69‑1 (MQ=255)
                                                        aCGGTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:3033123/69‑1 (MQ=255)
                                                        aCGGTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:1142363/69‑1 (MQ=255)
                                                        aCGGTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:1225785/69‑1 (MQ=255)
                                                        aCGGTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:1227746/69‑1 (MQ=255)
                                                        aCGGTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:145096/69‑1 (MQ=255)
                                                        aCGGTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:1740844/69‑1 (MQ=255)
                                                        aCGGTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:2571467/69‑1 (MQ=255)
                                                        aCGGTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:2837472/69‑1 (MQ=255)
                                                        aCGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCt                           >  1:2562144/1‑44 (MQ=255)
                                                        aCGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGAt                     >  1:2095894/1‑50 (MQ=255)
                                                        aCGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTa             >  1:653735/1‑58 (MQ=255)
                                                        aCGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  >  1:2603792/1‑69 (MQ=255)
                                                        aCGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  >  1:2514702/1‑69 (MQ=255)
                                                        aCGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  >  1:2247634/1‑69 (MQ=255)
                                                        aCGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  >  1:473321/1‑69 (MQ=255)
                                                        aCGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  >  1:2074157/1‑69 (MQ=255)
                                                        aCGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  >  1:69369/1‑69 (MQ=255)
                                                        aCGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  >  1:2940021/1‑69 (MQ=255)
                                                        aCGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  >  1:3073832/1‑69 (MQ=255)
                                                        aCGGTCGATGTTCACTGGTCGAGCCTGAACTAAAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  >  1:3013106/1‑69 (MQ=255)
                                                          ggTTGATGTTCACTGGTCGAGCCTGAACTACAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatc  <  1:2564038/67‑1 (MQ=255)
                                                           gTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCt                           <  1:177163/41‑1 (MQ=255)
                                                           gTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCt                           <  1:2785446/41‑1 (MQ=255)
                                                             |                                                               
AAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATC  >  minE/2218800‑2218924

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: