Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,211 A→G 10.5% G286G (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,2110AG10.5% 234.3 / 17.3 86G286G (GGT→GGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (38/39);  new base G (2/7);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

ATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCG  >  minE/1015144‑1015280
                                                                   |                                                                     
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:811356/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)
          aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                                                               <  1:2315201/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)
                         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:3387526/71‑1 (MQ=255)
                                               tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCcac                     <  1:3416954/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:1398539/1‑71 (MQ=255)
                                                 tttACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACAt                   <  1:773049/71‑1 (MQ=255)
                                                 tttACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACAt                   >  1:3581177/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:253724/42‑1 (MQ=255)
                                                       accaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAg                                          <  1:608130/42‑1 (MQ=255)
                                                       accaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTc                        >  1:3972340/1‑60 (MQ=255)
                                                        ccaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGAt                                             <  1:389891/38‑1 (MQ=255)
                                                        ccaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTc                              <  1:3769897/53‑1 (MQ=255)
                                                         caAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGc                                         <  1:2669928/41‑1 (MQ=255)
                                                         caAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAAc           <  1:1353792/71‑1 (MQ=255)
                                                          aaaCTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGaa            <  1:3727336/69‑1 (MQ=255)
                                                          aaaCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCtttta                               >  1:4038818/1‑49 (MQ=255)
                                                          aaaCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCa                            >  1:2288710/1‑53 (MQ=255)
                                                          aaaCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACg          >  1:689281/1‑71 (MQ=255)
                                                           aaCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCTGAACg          <  1:1013152/70‑1 (MQ=255)
                                                             cTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCt       <  1:3510793/71‑1 (MQ=255)
                                                             cTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCt       >  1:671241/1‑71 (MQ=255)
                                                               gcgcgcCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGa     >  1:2178208/1‑71 (MQ=255)
                                                               gcgcACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                      >  1:4035986/1‑54 (MQ=255)
                                                                cgcACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCttttt                               >  1:1249604/1‑44 (MQ=255)
                                                                cgcACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGAc    <  1:1670164/71‑1 (MQ=255)
                                                                cgcACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGAc    <  1:2262588/71‑1 (MQ=255)
                                                                 gcACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGAcc   >  1:832801/1‑71 (MQ=255)
                                                                  cACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCg  >  1:4093193/1‑71 (MQ=255)
                                                                  cACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCg  <  1:1430437/71‑1 (MQ=255)
                                                                   |                                                                     
ATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCG  >  minE/1015144‑1015280

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: