Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,680,438 A→C 33.0% L1512L (CTT→CTG yfhM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,680,4380AC33.0% 73.8 / 69.2 91L1512L (CTT→CTGyfhMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (46/15);  new base C (8/22);  total (54/37)
Fisher's exact test for biased strand distribution p-value = 1.22e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 107 total aligned reads displayed.

ACACCAGCACCAGATCGCCGCTACGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGC  >  minE/1680377‑1680506
                                                             |                                                                    
acacCAGCACCAGATCGCCGCTACGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGAatat                                                               >  1:2230593/1‑69 (MQ=255)
acacCAGCACCAGATCGCCGCTACGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGAatat                                                               >  1:1360146/1‑69 (MQ=255)
         ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGaa                                                                  >  1:869739/1‑57 (MQ=255)
         ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAatat                                                               >  1:2361200/1‑60 (MQ=255)
         ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGa                                                      >  1:2647119/1‑69 (MQ=255)
         ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGa                                                      >  1:2361780/1‑69 (MQ=255)
         ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGa                                                      >  1:2080480/1‑69 (MQ=255)
         ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGa                                                      >  1:1916803/1‑69 (MQ=255)
         ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGa                                                      >  1:368050/1‑69 (MQ=255)
         ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGa                                                      >  1:615176/1‑69 (MQ=255)
         ccagATCGCCGCTACGTAACGAGTCCAGCGACTTGCTCTTACCATCAGTACCAAGAATATGAAGCTCGa                                                      >  1:1291055/1‑69 (MQ=255)
                      aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                                  >  1:1067007/1‑44 (MQ=255)
                      aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                                  >  1:2992885/1‑44 (MQ=255)
                      aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                                  >  1:232692/1‑44 (MQ=255)
                      aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                                  >  1:228796/1‑44 (MQ=255)
                      aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                                  >  1:85235/1‑44 (MQ=255)
                      aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                                  >  1:786924/1‑44 (MQ=255)
                      aCGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGaa                                                                  >  1:1332502/1‑44 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:1475999/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:2634297/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:2977438/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:1168551/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:2354138/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:3168181/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:2191917/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:209842/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:3257791/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:2041503/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:431392/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:1544247/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:701273/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:630421/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:1606064/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:761137/69‑1 (MQ=255)
                          aaCGAGTCCAGCGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCACCACATTGCTc                                     <  1:2069956/69‑1 (MQ=255)
                                tCCAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTc                                     >  1:1783113/1‑63 (MQ=255)
                                tCCAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCgg                                >  1:193906/1‑68 (MQ=255)
                                tCCAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCgg                                >  1:3340585/1‑68 (MQ=255)
                                  cAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTaa                             <  1:2477789/69‑1 (MQ=255)
                                     cGATTTGCTCTTACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:1829940/58‑1 (MQ=255)
                                         ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:2419497/69‑1 (MQ=255)
                                         ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:2445363/69‑1 (MQ=255)
                                         ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:1417996/69‑1 (MQ=255)
                                         ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:2978992/69‑1 (MQ=255)
                                         ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:2801221/69‑1 (MQ=255)
                                         ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:266390/69‑1 (MQ=255)
                                         ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:1119687/69‑1 (MQ=255)
                                         ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:1493398/69‑1 (MQ=255)
                                         ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:126743/69‑1 (MQ=255)
                                         ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:1305825/69‑1 (MQ=255)
                                         ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:3254433/69‑1 (MQ=255)
                                         ttGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:3416577/69‑1 (MQ=255)
                                          tGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcg                      <  1:1487190/68‑1 (MQ=255)
                                                tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAgg                          >  1:786730/1‑58 (MQ=255)
                                                tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg               >  1:3134919/1‑69 (MQ=255)
                                                tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg               >  1:861807/1‑69 (MQ=255)
                                                tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg               >  1:541742/1‑69 (MQ=255)
                                                tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg               >  1:518988/1‑69 (MQ=255)
                                                tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg               >  1:2401994/1‑69 (MQ=255)
                                                tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg               >  1:2580015/1‑69 (MQ=255)
                                                tACCGTCAGTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTCGCCGGTAAAGGtgcggattgcg               >  1:2021227/1‑69 (MQ=255)
                                                 aCCATCAGTACCAATAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAgg                          <  1:3285114/57‑1 (MQ=255)
                                                 aCCATCAGTACCAAGAATATGCCGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcggattgcg               >  1:1790379/1‑68 (MQ=255)
                                                 aCCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGc                                   >  1:2876206/1‑48 (MQ=255)
                                                 aCCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcggattgcg               >  1:2842619/1‑68 (MQ=255)
                                                 aCCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcggattgcg               >  1:2956425/1‑68 (MQ=255)
                                                 aCCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcggattgcg               >  1:3429185/1‑68 (MQ=255)
                                                 aCCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGtgcggattgcg               >  1:2628371/1‑68 (MQ=255)
                                                 aCCATCAGTACCAAGAAAATGACGCTCGATTTGCAGCACATAGTTCGCCGGTAAAGGtgcggattgcg               >  1:1343807/1‑68 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATaa          >  1:1219209/1‑68 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAc         >  1:1722404/1‑69 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAc         >  1:174454/1‑69 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAc         >  1:194149/1‑69 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAc         >  1:1392279/1‑69 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAc         >  1:1185925/1‑69 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAc         >  1:2373438/1‑69 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAc         >  1:2517755/1‑69 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAc         >  1:299739/1‑69 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAc         >  1:2377492/1‑69 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAc         >  1:1002770/1‑69 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAc         >  1:2467901/1‑69 (MQ=255)
                                                      cAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAc         >  1:2814725/1‑69 (MQ=255)
                                                       aGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAgg                          >  1:1233534/1‑51 (MQ=255)
                                                       aGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAAcc        >  1:2956466/1‑69 (MQ=255)
                                                        gTACCCAGAATATGACGCTCAATTTGCAGCACATTGCTc                                     <  1:2071604/39‑1 (MQ=37)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:1284607/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:69363/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:1582317/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:933476/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:625569/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:1869846/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:1909310/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:3353364/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:2665908/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:2047792/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:321050/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:2339611/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:3150730/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:3047002/69‑1 (MQ=255)
                                                             aaGAATATGACGCTCGATTTGCAGCACATTGTCCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGc  <  1:1574756/69‑1 (MQ=255)
                                                             |                                                                    
ACACCAGCACCAGATCGCCGCTACGTAACGAGTCCAGCGATTTGCTCTTACCATCAGTACCAAGAATATGACGCTCGATTTGCAGCACATTGTTCGCCGGTAAAGGTGCGGATTGCGGATAACCGCTGGC  >  minE/1680377‑1680506

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