Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,184 A→C 17.4% A295A (GCT→GCG fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,1840AC17.4% 118.7 / 35.2 69A295A (GCT→GCGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base C (5/7);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

GACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTC  >  minE/1015121‑1015252
                                                               |                                                                    
gaCGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAg                                                                >  1:1770329/1‑70 (MQ=255)
  cGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAt                                                             >  1:3319399/1‑71 (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)
        tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                                       >  1:245509/1‑71 (MQ=255)
         accgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCg                                                                   >  1:946243/1‑58 (MQ=255)
          ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAtt                                                            >  1:1802242/1‑64 (MQ=255)
          ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                                       >  1:3443266/1‑69 (MQ=255)
          ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                                       >  1:401627/1‑69 (MQ=255)
          ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                                       >  1:672198/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)
                 ggCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTgc                                              >  1:1363658/1‑71 (MQ=255)
                  gCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTgcg                                             >  1:2227199/1‑71 (MQ=255)
                  gCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTgcg                                             >  1:2118624/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:3379144/66‑1 (MQ=255)
                                 aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                                   <  1:2955778/66‑1 (MQ=255)
                                 aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                                   <  1:2315201/66‑1 (MQ=255)
                                 aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                                   <  1:3892899/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:2869084/71‑1 (MQ=255)
                                         ggCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCt                      <  1:3157958/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)
                                              aTCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGAt                 <  1:1137028/71‑1 (MQ=255)
                                               tCACGCGGATGTCGTGAGCGAAGTATTTACCACCaaa                                                  >  1:18305/1‑37 (MQ=255)
                                               tCACGCGGATGTCGTGAGCGAAGTATTTACCACCaaa                                                  >  1:1743315/1‑37 (MQ=255)
                                               tCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTg                                               >  1:2131783/1‑40 (MQ=255)
                                                cACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCa                                           <  1:2163469/43‑1 (MQ=255)
                                                 aCGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATc                >  1:1280967/1‑69 (MQ=255)
                                                 aCGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCa               >  1:1857560/1‑70 (MQ=255)
                                                    cgGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCaa           <  1:4100535/71‑1 (MQ=255)
                                                      gATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAAt          <  1:2250479/70‑1 (MQ=255)
                                                         gTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAttt                            <  1:3256401/49‑1 (MQ=255)
                                                         gTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAttt                            <  1:2915426/49‑1 (MQ=255)
                                                          tCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTcctg     >  1:1487542/1‑68 (MQ=255)
                                                          tCGTGAGCGAAGTATTTACCACAAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCttt     <  1:1919404/71‑1 (MQ=255)
                                                             tGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCaa                                 <  1:2455580/40‑1 (MQ=255)
                                                              gAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCaa           >  1:1163824/1‑61 (MQ=255)
                                                              gAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTc  >  1:1704624/1‑70 (MQ=255)
                                                               aGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTc                    >  1:3123302/1‑51 (MQ=255)
                                                               |                                                                    
GACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTC  >  minE/1015121‑1015252

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: