Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,517 G→C 22.9% A184A (GCC→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,5170GC22.9% 230.3 / 63.0 122A184A (GCC→GCGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (44/50);  new base C (13/15);  total (57/65)
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.66e-01

Only 100 of 125 total aligned reads displayed.

GTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGT  >  minE/1015447‑1015585
                                                                      |                                                                    
gTTTCCCCAGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTg                                                                      >  1:1338241/1‑71 (MQ=255)
   tCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc                                                                    <  1:3608082/70‑1 (MQ=255)
   tCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc                                                                    <  1:3885495/70‑1 (MQ=255)
     cccGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGa                                                                  <  1:1886545/70‑1 (MQ=255)
     cccGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGa                                                                  <  1:862705/70‑1 (MQ=255)
        ggCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACa                                                              <  1:668128/71‑1 (MQ=255)
              aGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTgg                                                                     >  1:247607/1‑58 (MQ=255)
              aGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGaa                                                        >  1:5405961/1‑71 (MQ=255)
              aGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGaa                                                        >  1:1553698/1‑71 (MQ=255)
              aGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGaa                                                        >  1:3129290/1‑71 (MQ=255)
               gTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAt                                                       <  1:3288018/71‑1 (MQ=255)
                  gggCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTg                                                    >  1:4443382/3‑71 (MQ=255)
                  gggCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTg                                                    <  1:289407/69‑1 (MQ=255)
                  aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTgg                                                                     >  1:2317734/1‑54 (MQ=255)
                  aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGa                                                                  >  1:4311400/1‑57 (MQ=255)
                  aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    >  1:5656091/1‑71 (MQ=255)
                  aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    >  1:2861294/1‑71 (MQ=255)
                  aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    <  1:2210864/71‑1 (MQ=255)
                   cGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGt                                                   <  1:523433/71‑1 (MQ=255)
                    gCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTa                                                  <  1:1073069/71‑1 (MQ=255)
                     cTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTac                                                 >  1:279029/1‑70 (MQ=255)
                     cTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTac                                                 >  1:4437263/1‑70 (MQ=255)
                     cTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTAt                                                 <  1:2561958/71‑1 (MQ=255)
                      tttGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTg                                                    <  1:2037477/67‑1 (MQ=255)
                        tGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACtc                                               <  1:1782442/70‑1 (MQ=255)
                        tGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtc                                               <  1:1586922/70‑1 (MQ=255)
                        tGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGAGATACAGAGGAATTTGTATtc                                               <  1:1905209/70‑1 (MQ=255)
                         ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc                                                                    <  1:4533583/48‑1 (MQ=255)
                         ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc                                                                    <  1:288495/48‑1 (MQ=255)
                         ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgt                                             >  1:5297930/1‑71 (MQ=255)
                         ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgt                                             >  1:3400818/1‑71 (MQ=255)
                         ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcg                                              >  1:1547230/1‑70 (MQ=255)
                          gTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgt                                             <  1:4363897/70‑1 (MQ=255)
                           tttCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGAt                                                                 <  1:4619618/49‑1 (MQ=255)
                           tttCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGAt                                                                 <  1:5346723/49‑1 (MQ=255)
                           tttCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGaa                                                        >  1:5251844/1‑58 (MQ=255)
                             tCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGc                                                                    <  1:3300920/44‑1 (MQ=9)
                             tCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGc                                          >  1:2245741/1‑70 (MQ=255)
                             tCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGc                                          <  1:3758264/70‑1 (MQ=255)
                             tCCTGGTAGAGATACGTTTTGTTGGCAGAGCCCCGCCTTTCGCAACGCAAAGGaa                                                        <  1:1819799/55‑1 (MQ=255)
                             tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGa                                                                  >  1:921680/1‑46 (MQ=255)
                             tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGa                                                                  >  1:5035740/1‑46 (MQ=255)
                             tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgtc                                            <  1:3839973/68‑1 (MQ=255)
                             tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGc                                          >  1:2995995/1‑70 (MQ=255)
                             tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGc                                          <  1:4641737/70‑1 (MQ=255)
                             tCCTGATAGAGATACGTCTTGTTAGCCGAACCACCACCTTTGGCGa                                                                  >  1:638196/1‑46 (MQ=255)
                              ccTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCa                                        >  1:3948516/1‑71 (MQ=255)
                               cTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGt                                                   <  1:4192742/59‑1 (MQ=255)
                               cTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGt                                                   <  1:262458/59‑1 (MQ=255)
                                tggtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                                       >  1:3978994/4‑70 (MQ=255)
                                tggtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                                       >  1:1893027/4‑70 (MQ=255)
                                 gATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    <  1:1344612/56‑1 (MQ=255)
                                 gATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    <  1:1502902/56‑1 (MQ=255)
                                  gtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCTATc                                      <  1:2453718/68‑1 (MQ=255)
                                  gtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATc                                      <  1:184854/68‑1 (MQ=255)
                                  aTAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATc                                      <  1:1020268/69‑1 (MQ=255)
                                  aTAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATc                                      <  1:4116209/69‑1 (MQ=255)
                                  aTAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa                                    <  1:3230301/71‑1 (MQ=255)
                                   tAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                                       >  1:1283257/1‑67 (MQ=255)
                                        aTACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCa                               >  1:2072584/1‑70 (MQ=255)
                                         tACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTa                                                  <  1:896968/50‑1 (MQ=255)
                                         tACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCAATCAACGGCAt                              >  1:2007269/1‑70 (MQ=255)
                                          aCGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACtcatcg                                           >  1:791211/1‑56 (MQ=255)
                                          aCGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATcta                                    >  1:2418078/1‑61 (MQ=255)
                                          aCGTCTTGTTTGCCGAACCACCACCTTTGGCGATACaga                                                            >  1:1900631/1‑39 (MQ=255)
                                          aCGTCTTGTTTGCCGAACCACCACCTTTGGCGATACaga                                                            >  1:2055945/1‑39 (MQ=255)
                                           cGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTaca                           <  1:4390952/71‑3 (MQ=255)
                                            gTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    <  1:5427927/45‑1 (MQ=255)
                                            gTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    <  1:1648357/45‑1 (MQ=255)
                                              cTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAt                                                       <  1:1160647/40‑1 (MQ=255)
                                               ttGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTAt                                                 <  1:1750528/45‑1 (MQ=255)
                                               ttGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAga                         <  1:1196073/69‑1 (MQ=255)
                                                 gTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa                                    >  1:4382572/1‑56 (MQ=255)
                                                 gTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa                                    >  1:4162793/1‑56 (MQ=255)
                                                  tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                                     <  1:3911620/38‑1 (MQ=255)
                                                  tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                                     <  1:1592012/38‑1 (MQ=255)
                                                  tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                                     <  1:4485693/38‑1 (MQ=255)
                                                  tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                                     <  1:1429455/38‑1 (MQ=255)
                                                  tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                                     <  1:3977140/38‑1 (MQ=255)
                                                  tttGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcgat                    >  1:360859/1‑71 (MQ=255)
                                                  tttGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcgat                    >  1:4978529/1‑71 (MQ=255)
                                                  tttGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcga                     >  1:453718/1‑70 (MQ=255)
                                                    tGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGt                                                   >  1:3242740/1‑38 (MQ=255)
                                                       cGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTgcg               <  1:4692434/71‑1 (MQ=255)
                                                         aaCCACCACCTTTGGCGATACAGAGGAATTTGTATtcg                                              >  1:3157211/1‑38 (MQ=255)
                                                          accaccaccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTg            <  1:725818/71‑1 (MQ=255)
                                                           ccaccaccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTgg           <  1:4156516/71‑1 (MQ=255)
                                                                accTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCt                  <  1:2720264/59‑1 (MQ=255)
                                                                 ccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGc                                >  1:2795932/1‑44 (MQ=255)
                                                                 ccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                             <  1:4458752/47‑1 (MQ=255)
                                                                 ccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                             <  1:4584255/47‑1 (MQ=255)
                                                                 ccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGAtt     <  1:4524809/71‑1 (MQ=255)
                                                                 ccTTTCGCGATACAGAGGAATTTGTATTCGTCGCCa                                        >  1:1318878/1‑36 (MQ=255)
                                                                  cTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgat                        >  1:2622213/1‑51 (MQ=255)
                                                                  cTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTg    <  1:3775936/71‑1 (MQ=255)
                                                                   tttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTgg   >  1:2072133/1‑71 (MQ=255)
                                                                    ttGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa                                    >  1:2488213/1‑37 (MQ=255)
                                                                    ttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAga                         >  1:4013716/1‑48 (MQ=255)
                                                                    ttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAg                          <  1:4629740/47‑1 (MQ=255)
                                                                    ttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGt  >  1:226183/1‑71 (MQ=255)
                                                                      |                                                                    
GTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGT  >  minE/1015447‑1015585

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: