Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,184 A→C 38.9% 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,1840AC38.9% 46.9 / 148.3 121A295A (GCT→GCGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (31/43);  new base C (20/27);  total (51/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 126 total aligned reads displayed.

CAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCC  >  minE/1015116‑1015253
                                                                    |                                                                     
cagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGc                                                                     <  1:44878/71‑1 (MQ=255)
 agGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGcgcg                                                                    <  1:1667475/71‑1 (MQ=255)
  ggAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGcgcg                                                                    >  1:5622372/1‑70 (MQ=255)
  ggAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGa                                                                   <  1:4194267/71‑1 (MQ=255)
    agaCGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAg                                                                 <  1:550342/71‑1 (MQ=255)
    agaCGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAg                                                                 <  1:2178286/71‑1 (MQ=255)
    agaCGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAg                                                                 <  1:1877428/71‑1 (MQ=255)
       cGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAt                                                              >  1:515647/1‑71 (MQ=255)
        gCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTAtt                                                             <  1:2853639/71‑1 (MQ=255)
         cccATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTgcgc                                                                     <  1:2705980/62‑1 (MQ=255)
           cATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTac                                                          >  1:1605888/1‑71 (MQ=255)
           cATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTac                                                          >  1:3480242/1‑71 (MQ=255)
            atgccgaccgGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTacc                                                         <  1:5183224/68‑1 (MQ=255)
             tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAttt                                                            <  1:5133979/67‑1 (MQ=255)
             tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAttt                                                            <  1:2331947/67‑1 (MQ=255)
             tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                                        >  1:761644/1‑71 (MQ=255)
             tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                                        >  1:4369366/1‑71 (MQ=255)
             tACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                                        >  1:3721887/1‑71 (MQ=255)
               ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAtt                                                             >  1:3575939/1‑64 (MQ=255)
               ccgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTac                                                          <  1:5179/67‑1 (MQ=255)
               ccgaccgGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTaccgc                                                       >  1:4794508/1‑68 (MQ=255)
               ccgaccgGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTaccgc                                                       >  1:3481542/1‑68 (MQ=255)
                cgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTAccacca                                                     <  1:766822/71‑1 (MQ=255)
                cgaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTAccacca                                                     <  1:3918109/71‑1 (MQ=255)
                cgaccgGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTAccgcca                                                     <  1:4751028/71‑1 (MQ=255)
                 gaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTacca                                                        >  1:3642065/1‑67 (MQ=255)
                 gaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTAccacca                                                     >  1:1116641/1‑70 (MQ=255)
                 gaccgGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTAccacca                                                     >  1:4505366/1‑70 (MQ=255)
                 gaccgGGCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTAccgcca                                                     <  1:3716781/70‑1 (MQ=255)
                      ggCAGGATGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTacc                                                         >  1:5363327/1‑61 (MQ=255)
                        cAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAttt                                                            <  1:868233/56‑1 (MQ=255)
                        cAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTgcg                                              >  1:1080147/1‑70 (MQ=255)
                         aGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCaa                                                    <  1:5593544/63‑1 (MQ=255)
                           gATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCAcc                                          >  1:3962549/1‑71 (MQ=255)
                           gATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCAcc                                          <  1:4296266/71‑1 (MQ=255)
                           gATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCAcc                                          <  1:5543269/71‑1 (MQ=255)
                             tGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTa                                                           >  1:1404639/1‑52 (MQ=255)
                             tGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTa                                                           >  1:3893736/1‑52 (MQ=255)
                             tGCGCCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGcc                                          >  1:3596755/1‑69 (MQ=255)
                               cgcCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCaaa                                                   >  1:4584156/1‑58 (MQ=255)
                               cgcCGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAg                                       <  1:3368764/70‑1 (MQ=255)
                                 ccGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTa                                                               <  1:3373173/44‑1 (MQ=255)
                                 ccGTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAcc                                    <  1:2376671/71‑1 (MQ=255)
                                   gTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGAc                                     >  1:4623788/1‑68 (MQ=255)
                                   gTGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGt                                                                <  1:3361534/41‑1 (MQ=19)
                                    tGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTg                                                <  1:2570077/56‑1 (MQ=255)
                                    tGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGAc                                     <  1:502631/67‑1 (MQ=255)
                                    tGACGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAg                                 >  1:4065708/1‑71 (MQ=255)
                                      aCGTGGCAGGCGAATCACGCGGATGTCGTGAGCGaa                                                                  >  1:2317324/1‑36 (MQ=255)
                                      aCGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAAc                                                  >  1:1413253/1‑52 (MQ=255)
                                      aCGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAt                               >  1:2604271/1‑71 (MQ=255)
                                      aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAgtt                               <  1:5281815/71‑3 (MQ=255)
                                      aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAgtt                               <  1:4418975/71‑3 (MQ=255)
                                      aCGTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAgt                                <  1:2940748/70‑2 (MQ=255)
                                       cGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCaa                                                    <  1:4795824/49‑1 (MQ=255)
                                       cGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCaa                                                    <  1:2038483/49‑1 (MQ=255)
                                        gTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAttt                             <  1:2163767/71‑1 (MQ=255)
                                        gTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAttt                             <  1:1041444/71‑1 (MQ=255)
                                        gTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTac                                                          >  1:5358333/1‑42 (MQ=19)
                                        gTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTa                                                           >  1:5089089/1‑41 (MQ=19)
                                        gTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTAccgcca                                                     <  1:3620845/47‑1 (MQ=255)
                                        gTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTAccgcca                                                     <  1:2529111/47‑1 (MQ=255)
                                        gTGGCAGACGGATAACGCGAATGTCGTGCGCGAAGTATTTAccgcca                                                     <  1:2535727/47‑1 (MQ=255)
                                         tGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGt                                                                <  1:3063156/35‑1 (MQ=255)
                                         tGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAtt                                                             <  1:404026/38‑1 (MQ=255)
                                         tGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTAt                                                              <  1:653327/37‑1 (MQ=255)
                                         tGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAtttt                            <  1:3543896/71‑1 (MQ=255)
                                         tGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAtttt                            <  1:2860717/71‑1 (MQ=255)
                                               gCGAATCACGCGGATGTCGTGAGCGAAGTATTTAccacca                                                     <  1:2215510/40‑1 (MQ=255)
                                               gCGAATCACGCGGATGTCGTGAGCGAAGTATTTAccacca                                                     <  1:2625634/40‑1 (MQ=255)
                                               gCGAATCACGCGGATGTCGTGAGCGAAGTATTTAccacca                                                     <  1:1577095/40‑1 (MQ=255)
                                               gCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAtt                              <  1:2447383/63‑1 (MQ=255)
                                               gCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAtt                              <  1:2112087/63‑1 (MQ=255)
                                                cGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAtttt                            <  1:2560277/64‑1 (MQ=255)
                                                 gAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAAc                                                  >  1:79475/1‑41 (MQ=255)
                                                 gAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCg                    <  1:1013409/71‑1 (MQ=255)
                                                 gAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCg                    <  1:2654855/71‑1 (MQ=255)
                                                  aaTCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAg                                       >  1:5381085/1‑51 (MQ=255)
                                                   aTCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGAt                  >  1:1184460/1‑71 (MQ=255)
                                                    tCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGAt                  <  1:4131335/70‑1 (MQ=255)
                                                    tCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATc                 <  1:2829265/71‑1 (MQ=255)
                                                    tCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATc                 <  1:2457002/71‑1 (MQ=255)
                                                      aCGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCAc                                           >  1:2862767/1‑43 (MQ=255)
                                                      aCGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAg               >  1:187400/1‑71 (MQ=255)
                                                      aCGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAg               >  1:3904013/1‑71 (MQ=255)
                                                       cgcgGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGc              <  1:3186520/71‑1 (MQ=255)
                                                       cgcgAATTTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGtttca                            >  1:4392976/1‑55 (MQ=255)
                                                       cgcgAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGttt                              <  1:5358661/55‑1 (MQ=255)
                                                       cgcgAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGttt                              <  1:1347161/55‑1 (MQ=255)
                                                       cgcgAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGttt                              <  1:828371/55‑1 (MQ=255)
                                                       cgcgAATGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTggg                         >  1:3676836/1‑58 (MQ=255)
                                                        gcgGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGc              >  1:2286505/1‑70 (MQ=255)
                                                            aTGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTc         <  1:2320369/71‑1 (MQ=255)
                                                            aTGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGAt                               <  1:1610945/49‑1 (MQ=255)
                                                             tGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCt        <  1:2910283/71‑1 (MQ=255)
                                                               tCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGtttct                            >  1:199845/1‑47 (MQ=19)
                                                               tCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTcctg      >  1:5068033/1‑68 (MQ=255)
                                                               tCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTcctg      >  1:2990150/1‑68 (MQ=255)
                                                               tCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGAt                  <  1:3395182/59‑1 (MQ=255)
                                                                    aGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTcc  <  1:89427/70‑1 (MQ=255)
                                                                    |                                                                     
CAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCC  >  minE/1015116‑1015253

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