Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,202,586 G→A 51.4% A82T (GCC→ACC)  yoaD → predicted phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,202,5860GA51.4% ‑5.1 / 62.0 68A82T (GCC→ACC) yoaDpredicted phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (8/25);  new base A (9/26);  total (17/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

ACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATC  >  minE/1202522‑1202635
                                                                |                                                 
aCTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCa                                               <  1:972433/69‑1 (MQ=255)
aCTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCa                                               <  1:2679783/69‑1 (MQ=255)
aCTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCa                                               <  1:2679460/69‑1 (MQ=255)
aCTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCa                                               <  1:267266/69‑1 (MQ=255)
aCTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCa                                               <  1:2595282/69‑1 (MQ=255)
aCTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCa                                               <  1:2448008/69‑1 (MQ=255)
aCTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCa                                               <  1:2155684/69‑1 (MQ=255)
aCTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGGCTGCCCTGCTCTGTCGCCCa                                               <  1:3339850/69‑1 (MQ=255)
      tCCCCTACAGGCAGTTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGc                                          <  1:312436/68‑1 (MQ=255)
      tCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGc                                          <  1:3244347/68‑1 (MQ=255)
      tCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGc                                          <  1:2225623/68‑1 (MQ=255)
      tCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGc                                          <  1:2007374/68‑1 (MQ=255)
         ccTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCAt                                       <  1:3309051/68‑1 (MQ=255)
         ccTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCAt                                       <  1:3210649/68‑1 (MQ=255)
               ggcaggTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTa                                 <  1:1039395/68‑1 (MQ=255)
               ggcaggTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTa                                 <  1:1541011/68‑1 (MQ=255)
                 caggTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGc                                          <  1:1454050/57‑1 (MQ=255)
                  aggTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTTCTCTGTCGCCCATTTGc                                          <  1:1970313/56‑1 (MQ=255)
                     tAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAggc                          >  1:323994/1‑69 (MQ=255)
                     tAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAggc                          >  1:1813993/1‑69 (MQ=255)
                     tAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAggc                          >  1:607559/1‑69 (MQ=255)
                      aGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAggc                          >  1:2370722/1‑68 (MQ=255)
                      aGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAggc                          >  1:1367640/1‑68 (MQ=255)
                      aGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAggc                          >  1:1503345/1‑68 (MQ=255)
                      aGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAggc                          >  1:2017409/1‑68 (MQ=255)
                          aaGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCaa                     <  1:621840/69‑1 (MQ=255)
                          aaGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCaa                     <  1:1377379/69‑1 (MQ=255)
                          aaGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCaa                     <  1:2634786/69‑1 (MQ=255)
                          aaGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCaa                     <  1:2697380/69‑1 (MQ=255)
                              cTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAAc                  >  1:2787484/1‑68 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGTCTGCAAAACTc                <  1:1901326/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAAccc                <  1:1372970/69‑3 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:3097307/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:3227414/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:3158035/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:323362/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:1438002/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:2885389/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:1220755/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:3432097/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:417772/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:645866/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:835322/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:953058/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:980682/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:3026839/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:2960414/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:2958477/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:1691619/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:1596793/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:2258055/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:2356869/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:1501624/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:2652269/69‑1 (MQ=255)
                               ttGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:2589281/69‑1 (MQ=255)
                                tGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTc                <  1:158955/68‑1 (MQ=255)
                                 gCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTcc               >  1:847771/1‑68 (MQ=255)
                                 gCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTcc               >  1:2439774/1‑68 (MQ=255)
                                 gCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTCCa              >  1:181076/1‑69 (MQ=255)
                                 gCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTCCa              >  1:3233924/1‑69 (MQ=255)
                                 gCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTCCa              >  1:2778633/1‑69 (MQ=255)
                                 gCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTCCa              >  1:2848548/1‑69 (MQ=255)
                                 gCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTCCa              >  1:3151184/1‑69 (MQ=255)
                                 gCTTCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTCCa              >  1:3008970/1‑69 (MQ=255)
                                 gCTCCCGCTGATTGGTCTGCCCTGCTCTGTCACCCATTTGCCATTACGTAAACAGGCTGCAAAACTCCa              >  1:146271/1‑69 (MQ=255)
                                   ttCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCTTTTGc                                          <  1:1032168/39‑1 (MQ=255)
                                     ccGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCAATACGTAAACAGGCGGCaa                     <  1:1661192/58‑1 (MQ=255)
                                              ggTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATc  <  1:271174/68‑1 (MQ=255)
                                                                |                                                 
ACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATC  >  minE/1202522‑1202635

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: