Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,004,772 G→T 36.4% G326G (GGC→GGA hybC ← hydrogenase 2, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,004,7720GT36.4% 88.0 / 69.4 96G326G (GGC→GGAhybChydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base T (4/31);  total (37/59)
Fisher's exact test for biased strand distribution p-value = 2.74e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

GTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATT  >  minE/2004704‑2004834
                                                                    |                                                              
gtggtgCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATg                                                                >  1:3244030/1‑69 (MQ=255)
gtggtgCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATg                                                                >  1:2925050/1‑69 (MQ=255)
gtggtgCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATg                                                                >  1:2880988/1‑69 (MQ=255)
      ccTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTg                                                          <  1:2902299/69‑1 (MQ=255)
                 tGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTc                                               >  1:3000248/1‑69 (MQ=255)
                  gCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCa                                              <  1:1048737/69‑1 (MQ=255)
                  gCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCa                                              <  1:2794460/69‑1 (MQ=255)
                  gCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCa                                              <  1:2342933/69‑1 (MQ=255)
                  gCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCa                                              >  1:223352/1‑69 (MQ=255)
                  gCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCa                                              <  1:2202425/69‑1 (MQ=255)
                      gcggcgCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCgg                                          <  1:3417239/69‑1 (MQ=255)
                      gcggcgCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCgg                                          <  1:882088/69‑1 (MQ=255)
                      gcggcgCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCgg                                          <  1:1720432/69‑1 (MQ=255)
                         gcgcTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGa                                         <  1:3279744/67‑1 (MQ=255)
                         gcgcTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGa                                         <  1:2257632/67‑1 (MQ=255)
                         gcgcTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGa                                         <  1:801478/67‑1 (MQ=255)
                              tCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGaa                                  >  1:2351130/1‑69 (MQ=255)
                              tCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGaa                                  >  1:3358068/1‑69 (MQ=255)
                              tCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGaa                                  >  1:3045030/1‑69 (MQ=255)
                              tCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGaa                                  >  1:2543018/1‑69 (MQ=255)
                              tCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATgaga                                   >  1:3029794/1‑68 (MQ=255)
                              tCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGaa                                  >  1:1023177/1‑69 (MQ=255)
                              tCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGaa                                  >  1:2258050/1‑69 (MQ=255)
                              tCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGaa                                  >  1:1558994/1‑69 (MQ=255)
                               cGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCa                                              <  1:846055/56‑1 (MQ=255)
                               cGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCgg                                          <  1:2349832/60‑1 (MQ=255)
                                  cTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTg                               <  1:1487692/68‑1 (MQ=255)
                                   tttATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGtt                               <  1:187034/67‑2 (MQ=255)
                                     tataCCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCg                           >  1:2532770/1‑69 (MQ=255)
                                     tataCCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCg                           >  1:86041/1‑69 (MQ=255)
                                     tataCCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCg                           >  1:833279/1‑69 (MQ=255)
                                               gTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACg                  >  1:2652404/1‑68 (MQ=255)
                                               gTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACg                  >  1:2747587/1‑68 (MQ=255)
                                               gTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACg                  >  1:2527284/1‑68 (MQ=255)
                                               gTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACg                  >  1:2410598/1‑68 (MQ=255)
                                               gTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACg                  >  1:1528982/1‑68 (MQ=255)
                                               gTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACg                  >  1:1084793/1‑68 (MQ=255)
                                               gTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACg                  >  1:2282996/1‑68 (MQ=255)
                                               gTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACg                  >  1:819812/1‑68 (MQ=255)
                                               gTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACg                  >  1:1081557/1‑68 (MQ=255)
                                               gTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACg                  >  1:2089178/1‑68 (MQ=255)
                                               gTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACg                  >  1:498446/1‑68 (MQ=255)
                                               gTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACg                  >  1:64544/1‑68 (MQ=255)
                                                tGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAg                >  1:2586371/1‑69 (MQ=255)
                                                tGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAg                >  1:300676/1‑69 (MQ=255)
                                                tGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAg                >  1:710422/1‑69 (MQ=255)
                                                tGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAg                >  1:1144213/1‑69 (MQ=255)
                                                tGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAg                >  1:2870242/1‑69 (MQ=255)
                                                tGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAg                >  1:653507/1‑69 (MQ=255)
                                                tGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAg                >  1:456992/1‑69 (MQ=255)
                                                tGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAg                >  1:2101247/1‑69 (MQ=255)
                                                tGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAg                >  1:21957/1‑69 (MQ=255)
                                                tGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAg                >  1:3100258/1‑69 (MQ=255)
                                                tGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAg                >  1:472566/1‑69 (MQ=255)
                                                tGCTTCGCGCTTTCCTGAATGCCATTGATCAGGTATTCATCgg                                          >  1:1089587/1‑43 (MQ=255)
                                                    tCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGtt                               <  1:183067/50‑2 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCATGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:2552685/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:3364075/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:989716/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:896180/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:3434336/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:385665/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:465961/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:540274/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:560855/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:641423/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:694102/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:709661/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:3234267/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:1280553/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:1307439/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:1405017/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:1438349/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:1566097/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:1635602/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:1666154/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:2211518/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:2273894/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:2344673/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:1117142/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:326932/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:2751489/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:2901259/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:3108792/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:3202973/60‑1 (MQ=255)
                                                          tttCCTGAATTCCTTTGATCAGGCATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:856310/60‑1 (MQ=255)
                                                          tttCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCg      <  1:2341721/69‑1 (MQ=255)
                                                          tttCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCg      <  1:527095/69‑1 (MQ=255)
                                                           ttCCTGAATTCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAgg               <  1:2338377/59‑1 (MQ=255)
                                                               tGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCAtt  <  1:1141613/68‑1 (MQ=255)
                                                               tGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCAtt  <  1:117100/68‑1 (MQ=255)
                                                               tGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCAtt  <  1:853561/68‑1 (MQ=255)
                                                               tGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCAtt  <  1:109694/68‑1 (MQ=255)
                                                               tGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCAtt  <  1:149255/68‑1 (MQ=255)
                                                               tGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCAtt  <  1:709918/68‑1 (MQ=255)
                                                               tGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCAtt  <  1:2028491/68‑1 (MQ=255)
                                                               tGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCAtt  <  1:2377985/68‑1 (MQ=255)
                                                               tGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCAtt  <  1:2660959/68‑1 (MQ=255)
                                                               tGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCAtt  <  1:456933/68‑1 (MQ=255)
                                                                    |                                                              
GTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTATTCATCGGAATGAGAAGTGATCGGACGATACGAGGACAGATCCGCATT  >  minE/2004704‑2004834

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: