Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,166 A→G 26.4% I301I (ATT→ATC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,1660AG26.4% 117.0 / 49.8 72I301I (ATT→ATCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base G (6/13);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTC  >  minE/1015097‑1015234
                                                                     |                                                                    
gATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGa                                                                      <  1:4014703/70‑1 (MQ=255)
    ttACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCa                                                                  <  1:3981594/70‑1 (MQ=255)
      aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgc                                                               >  1:911659/1‑71 (MQ=255)
       cGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAAcgcg                                                              <  1:2736628/71‑1 (MQ=255)
        ggTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgcg                                                              <  1:2175118/70‑1 (MQ=255)
             gcagagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTc                                                        <  1:2543808/71‑1 (MQ=255)
              cggagcagGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCg                                                       >  1:1212214/3‑71 (MQ=255)
              caaagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCg                                                       >  1:3985808/4‑71 (MQ=255)
                 agcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGt                                                      >  1:1554758/1‑69 (MQ=255)
                      gagaCGCCCATGCCGACCGGGCAGGATGCCCCGTGACGTGGCAGACGGATAACGCGAATg                                                          <  1:3274067/60‑1 (MQ=255)
                       agaCGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATg                                                          <  1:14321/59‑1 (MQ=255)
                       agaCGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCActcg                                                              <  1:2696733/55‑1 (MQ=255)
                        gaCGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAg                                              >  1:1770329/1‑70 (MQ=255)
                          cGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAt                                           >  1:3319399/1‑71 (MQ=255)
                           gCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATg                                                          <  1:638669/55‑1 (MQ=255)
                            cccATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTAttt                                         <  1:1809316/71‑1 (MQ=255)
                              cATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGaa                                               >  1:655367/1‑63 (MQ=255)
                              cATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTac                                       <  1:1319905/71‑1 (MQ=255)
                              cATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTac                                       <  1:2589799/71‑1 (MQ=255)
                               aTACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacc                                      <  1:2001476/71‑1 (MQ=255)
                                tgccgaccgGGCAGGATGCGCCGTGACGTGGCAGACGGAt                                                                    <  1:3496351/38‑1 (MQ=37)
                                tgccgaccgGGCAGGATGCGCCGTGACGTGGCAGACGGAt                                                                    <  1:1277961/38‑1 (MQ=37)
                                tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgcg                                                              <  1:242371/46‑1 (MQ=255)
                                tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                     >  1:245509/1‑71 (MQ=255)
                                 accgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcg                                                                >  1:1644411/1‑43 (MQ=255)
                                 accgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCg                                                 >  1:946243/1‑58 (MQ=255)
                                  ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAc                                                                 <  1:2034086/41‑1 (MQ=255)
                                  ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGa                                                            >  1:2877680/1‑46 (MQ=255)
                                  ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAtt                                          >  1:1802242/1‑64 (MQ=255)
                                  ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                     >  1:3443266/1‑69 (MQ=255)
                                  ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                     >  1:672198/1‑69 (MQ=255)
                                  ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                     >  1:401627/1‑69 (MQ=255)
                                  ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTa                                        >  1:2288802/1‑66 (MQ=255)
                                    gaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTAccacca                                  >  1:3647370/1‑70 (MQ=255)
                                       cgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTa                                            <  1:2279056/57‑1 (MQ=255)
                                       cgGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACt                              >  1:2804835/1‑71 (MQ=255)
                                       cgGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGAGCGAAGTATTTACCGCCAAACt                              >  1:423132/1‑71 (MQ=255)
                                         ggCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCgg                                                             >  1:2236292/1‑38 (MQ=255)
                                         ggCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTgc                            >  1:1363658/1‑71 (MQ=255)
                                          gCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTgcg                           >  1:2118624/1‑71 (MQ=255)
                                          gCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTgcg                           >  1:2227199/1‑71 (MQ=255)
                                          gCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCaa                                 <  1:2857572/65‑1 (MQ=255)
                                              gATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCg                                                 <  1:1907417/45‑1 (MQ=255)
                                               aTGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCAcc                       >  1:3522643/1‑70 (MQ=255)
                                                tGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAttt                                         <  1:3428559/51‑1 (MQ=255)
                                                tGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCAcc                       >  1:3854277/1‑69 (MQ=255)
                                                tGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCa                     <  1:3799115/71‑1 (MQ=255)
                                                tGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAg                                              >  1:2012674/1‑46 (MQ=255)
                                                tGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGcc                       >  1:3839792/1‑69 (MQ=255)
                                                    ccGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGAcc                 <  1:3186068/71‑1 (MQ=255)
                                                      gTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGAc                  >  1:531074/1‑68 (MQ=255)
                                                       tGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTa                                            <  1:1768903/41‑1 (MQ=255)
                                                       tGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTa                                            <  1:1357955/41‑1 (MQ=255)
                                                        gACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGAcc                 >  1:2128846/1‑67 (MQ=255)
                                                         aCGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAt            >  1:320035/1‑71 (MQ=255)
                                                         aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                 <  1:811356/66‑1 (MQ=255)
                                                         aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                 <  1:2315201/66‑1 (MQ=255)
                                                         aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                 <  1:3892899/66‑1 (MQ=255)
                                                         aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                 <  1:3379144/66‑1 (MQ=255)
                                                         aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                 <  1:2955778/66‑1 (MQ=255)
                                                           gTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAttt          <  1:1067884/71‑1 (MQ=255)
                                                             ggCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTg        >  1:3628426/1‑71 (MQ=255)
                                                              gCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTg        <  1:2654329/70‑1 (MQ=255)
                                                              gCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTg        <  1:2361308/70‑1 (MQ=255)
                                                                aGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTgcgc     <  1:3627870/71‑1 (MQ=255)
                                                                 ggCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCt    <  1:3157958/71‑1 (MQ=255)
                                                                 ggCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCt    <  1:2869084/71‑1 (MQ=255)
                                                                 ggCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCt    <  1:1201057/71‑1 (MQ=255)
                                                                  gCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCt    <  1:3381379/70‑1 (MQ=255)
                                                                   cGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTc  >  1:1767833/1‑71 (MQ=255)
                                                                   cGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTc  <  1:2394728/71‑1 (MQ=255)
                                                                   cGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCACACCCAGACCAAGATTTTGCGCTTc  <  1:360509/71‑1 (MQ=255)
                                                                    gAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAtt           >  1:2044610/1‑61 (MQ=255)
                                                                     |                                                                    
GATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTC  >  minE/1015097‑1015234

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: