Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,130 A→G 24.2% G313G (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,1300AG24.2% 111.4 / 35.9 62G313G (GGT→GGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base G (7/8);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCA  >  minE/1015064‑1015199
                                                                  |                                                                     
gATCCAGATACCTTCGCGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGGGCAGGAGACGCCCATgccg                                                                    <  1:676450/70‑1 (MQ=255)
 aTCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgac                                                                  >  1:2805520/1‑71 (MQ=255)
   ccAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgaccg                                                                >  1:1803164/1‑71 (MQ=255)
      gATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCggg                                                              >  1:21813/1‑70 (MQ=255)
        tCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAg                                                           <  1:2251038/71‑1 (MQ=255)
              cgcGGTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATgcg                                                     >  1:2523653/3‑71 (MQ=255)
               aCGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgcgc                                                    >  1:2693676/1‑71 (MQ=255)
               aCGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc                                                      <  1:366663/69‑1 (MQ=255)
               aCGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc                                                      <  1:864029/69‑1 (MQ=255)
                  gTTGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                                                 >  1:767302/1‑71 (MQ=255)
                  gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGt                                                 <  1:3655157/71‑1 (MQ=255)
                  gTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGt                                                 <  1:2306050/71‑1 (MQ=255)
                   ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                                                 >  1:163156/1‑70 (MQ=255)
                   ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                                                 >  1:2901124/1‑70 (MQ=255)
                   ttGATTTTCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGt                                                 >  1:1126405/1‑70 (MQ=255)
                   ttGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGt                                                 >  1:1146590/1‑70 (MQ=255)
                    tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTg                                                >  1:626491/1‑70 (MQ=255)
                    tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                                               >  1:2695710/1‑71 (MQ=255)
                    tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                                               >  1:2511486/1‑71 (MQ=255)
                    tGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                                               >  1:2700594/1‑71 (MQ=255)
                       ttttCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCggg                                                              <  1:3164221/51‑1 (MQ=255)
                       ttttCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCggg                                                              <  1:3168649/51‑1 (MQ=255)
                       tCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGt                                            >  1:1859979/1‑71 (MQ=255)
                        tttCGCTTTAATGTTACGGTCAGCGGAGCAGGAGACGCCCATGCCGACCggg                                                              <  1:2070230/51‑1 (MQ=255)
                        cTTCGCTTTTATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGa                                               <  1:173596/67‑1 (MQ=255)
                        cTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg                                           >  1:2779622/1‑71 (MQ=255)
                        cTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg                                           >  1:2658490/1‑71 (MQ=255)
                         ttCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATgc                                                      <  1:2952554/59‑1 (MQ=255)
                         ttCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTgg                                          >  1:4109445/1‑71 (MQ=255)
                         ttCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTgg                                          >  1:2984807/1‑71 (MQ=255)
                           cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACg                                             <  1:529128/66‑1 (MQ=255)
                           cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACg                                             <  1:3364572/66‑1 (MQ=255)
                           cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg                                           <  1:771055/68‑1 (MQ=255)
                           cGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg                                           <  1:2939146/68‑1 (MQ=255)
                           cGCTTTGATATTACGGGCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg                                           <  1:2031427/68‑1 (MQ=255)
                            gCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCa                                        <  1:209572/70‑1 (MQ=255)
                                tGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCg                                    >  1:835920/1‑70 (MQ=255)
                                tGATATTACGGGCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCg                                    <  1:1151465/70‑1 (MQ=255)
                                 gATATTACGGTCAGCAGAGCAGGAGACGCCCATaccgac                                                                  >  1:1671867/1‑39 (MQ=255)
                                 gATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGa                                   <  1:4014703/70‑1 (MQ=255)
                                     ttACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAg                                                           <  1:3241916/42‑1 (MQ=255)
                                     ttACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGa                                                         >  1:3993592/1‑44 (MQ=255)
                                     ttACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCa                               <  1:3981594/70‑1 (MQ=255)
                                       aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGc                                     >  1:3357711/1‑62 (MQ=255)
                                       aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGc                                     >  1:1534337/1‑62 (MQ=255)
                                       aCGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgc                            >  1:911659/1‑71 (MQ=255)
                                        cGGTCAGCGGAGCAGGAGACGCCCATGCCGACCGGGCAGGATGCGCCGTGACGTGGCAGACGGATAAcgcg                           <  1:2736628/71‑1 (MQ=255)
                                         ggTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCAcgcg                           <  1:2175118/70‑1 (MQ=255)
                                             agcagagcagGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTg                                           <  1:965196/50‑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)
                                                                 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)
                                                                  |                                                                     
GATCCAGATCCCCTGACGGTTGATCTTCGCTTTGATATTACGGTCAGCAGAGCAGGAGACGCCCATACCGACCGGGCAGGATGCGCCGTGACGTGGCAGGCGAATCACGCGGATGTCGTGAGCGAAGTATTTACCA  >  minE/1015064‑1015199

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: