Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,199 A→G 15.1% G290G (GGT→GGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,1990AG15.1% 177.3 / 23.8 73G290G (GGT→GGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (34/28);  new base G (4/7);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

TACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCG  >  minE/1015129‑1015267
                                                                      |                                                                    
tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                                                      >  1:245509/1‑71 (MQ=255)
  ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                                                      >  1:401627/1‑69 (MQ=255)
  ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                                                      >  1:3443266/1‑69 (MQ=255)
  ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                                                      >  1:672198/1‑69 (MQ=255)
    gaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTAccacca                                                                   >  1:3647370/1‑70 (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:2118624/1‑71 (MQ=255)
          gCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTgcg                                                            >  1:2227199/1‑71 (MQ=255)
          gCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCaa                                                                  <  1:2857572/65‑1 (MQ=255)
               aTGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCAcc                                                        >  1:3522643/1‑70 (MQ=255)
                tGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCAcc                                                        >  1:3854277/1‑69 (MQ=255)
                tGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCa                                                      <  1:3799115/71‑1 (MQ=255)
                tGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGcc                                                        >  1:3839792/1‑69 (MQ=255)
                    ccGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGAcc                                                  <  1:3186068/71‑1 (MQ=255)
                      gTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGAc                                                   >  1:531074/1‑68 (MQ=255)
                        gACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGAcc                                                  >  1:2128846/1‑67 (MQ=255)
                         aCGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAt                                             >  1:320035/1‑71 (MQ=255)
                         aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                                                  <  1:2955778/66‑1 (MQ=255)
                         aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                                                  <  1:3892899/66‑1 (MQ=255)
                         aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                                                  <  1:3379144/66‑1 (MQ=255)
                         aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                                                  <  1:2315201/66‑1 (MQ=255)
                         aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                                                  <  1:811356/66‑1 (MQ=255)
                           gTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAttt                                           <  1:1067884/71‑1 (MQ=255)
                             ggCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTg                                         >  1:3628426/1‑71 (MQ=255)
                              gCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTg                                         <  1:2361308/70‑1 (MQ=255)
                              gCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTg                                         <  1:2654329/70‑1 (MQ=255)
                                aGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTgcgc                                      <  1:3627870/71‑1 (MQ=255)
                                 ggCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCt                                     <  1:1201057/71‑1 (MQ=255)
                                 ggCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCt                                     <  1:3157958/71‑1 (MQ=255)
                                 ggCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCt                                     <  1:2869084/71‑1 (MQ=255)
                                  gCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCt                                     <  1:3381379/70‑1 (MQ=255)
                                   cGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTc                                   <  1:2394728/71‑1 (MQ=255)
                                   cGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTc                                   >  1:1767833/1‑71 (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)
                                                          gAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATc                               >  1:2203206/1‑52 (MQ=255)
                                                           aaGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAgctg           <  1:1793654/71‑4 (MQ=255)
                                                             gTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTc                 <  1:3955348/63‑1 (MQ=255)
                                                             gTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTc           <  1:3830387/69‑1 (MQ=255)
                                                             gTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca         >  1:1588141/1‑71 (MQ=255)
                                                              tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTgcgc                                      >  1:1071623/1‑41 (MQ=255)
                                                              tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCg                                  <  1:1431807/45‑1 (MQ=255)
                                                              tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCa                              <  1:717190/49‑1 (MQ=255)
                                                              tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCcac        <  1:3416954/71‑1 (MQ=255)
                                                              tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCcac        <  1:3387526/71‑1 (MQ=255)
                                                               aTTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCttt                    <  1:964570/58‑1 (MQ=255)
                                                               aTTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCcaca       >  1:1579372/1‑71 (MQ=255)
                                                                tttACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGAt                                >  1:865288/1‑45 (MQ=255)
                                                                tttACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACAt      >  1:3581177/1‑71 (MQ=255)
                                                                tttACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACAt      <  1:773049/71‑1 (MQ=255)
                                                                tttACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACAt      >  1:1398539/1‑71 (MQ=255)
                                                                tttACCACCAAACTGCGCACCAAGACCAAGATTTTGCGCTTCGATCAGc                            >  1:2485430/1‑49 (MQ=255)
                                                                  tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTcc                >  1:3783762/1‑59 (MQ=255)
                                                                  tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCa               >  1:3472449/1‑60 (MQ=255)
                                                                  tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCa               >  1:952317/1‑60 (MQ=255)
                                                                  tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAg              >  1:3343906/1‑61 (MQ=255)
                                                                  tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcg    >  1:2168905/1‑71 (MQ=255)
                                                                  tACCACCAAACTGCGCACCCAGACCAAGATTTTGAGCTTCGATCAGCAATTCTTTTTCCa               >  1:254772/1‑60 (MQ=255)
                                                                  tACCACAAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCa                           >  1:116994/1‑48 (MQ=255)
                                                                    ccaccaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcgcg  >  1:2750842/1‑71 (MQ=255)
                                                                    ccaccaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcgc   <  1:2035775/70‑1 (MQ=255)
                                                                      accaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAg                             <  1:608130/42‑1 (MQ=255)
                                                                      accaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAg                             <  1:253724/42‑1 (MQ=255)
                                                                      accaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTc           >  1:3972340/1‑60 (MQ=255)
                                                                      |                                                                    
TACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCG  >  minE/1015129‑1015267

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: