Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,508 A→G 26.9% G187G (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,5080AG26.9% 167.6 / 69.6 110G187G (GGT→GGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (37/43);  new base G (13/16);  total (50/60)
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

Only 100 of 121 total aligned reads displayed.

AATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGG  >  minE/1015438‑1015576
                                                                      |                                                                    
aaTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAacca                                                                      >  1:3940664/1‑71 (MQ=255)
 aTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccac                                                                     <  1:4047212/71‑1 (MQ=255)
 aTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccac                                                                     <  1:3980521/71‑1 (MQ=255)
   ttttCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccacc                                                                    <  1:3273283/70‑1 (MQ=255)
    tttCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAgccgccgc                                                                  >  1:5472185/1‑69 (MQ=255)
       cAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCt                                                                <  1:3567305/70‑1 (MQ=255)
        aGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCtt                                                               <  1:1191864/70‑1 (MQ=255)
         gTTTCCCCAGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTg                                                             >  1:1338241/1‑71 (MQ=255)
            tCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc                                                           <  1:3885495/70‑1 (MQ=255)
            tCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc                                                           >  1:1573830/1‑70 (MQ=255)
              cccGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGa                                                         <  1:862705/70‑1 (MQ=255)
              cccGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGa                                                         <  1:1886545/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:3129290/1‑71 (MQ=255)
                        gTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAt                                              <  1:3288018/71‑1 (MQ=255)
                         tAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccaccac                                                                  <  1:629610/50‑1 (MQ=255)
                         tAACGCTTTGGGTTCCTGATAGAGATACGTCTTGTTTGCCGAaccaccac                                                                  <  1:1187300/50‑1 (MQ=255)
                          aaCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccaccac                                                                  <  1:1484955/49‑1 (MQ=255)
                           gggCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTg                                           >  1:3761804/3‑71 (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:2210864/71‑1 (MQ=255)
                           aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                           >  1:2861294/1‑71 (MQ=255)
                            cGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGt                                          <  1:523433/71‑1 (MQ=255)
                             gCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTa                                         <  1:1073069/71‑1 (MQ=255)
                             gCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTa                                         <  1:2198266/71‑1 (MQ=255)
                              cTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCttt                                                              >  1:1305058/1‑49 (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:4594302/70‑1 (MQ=255)
                                 tGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACtc                                      <  1:1782442/70‑1 (MQ=255)
                                 tGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACtc                                      <  1:1224949/70‑1 (MQ=255)
                                 tGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCtt                                                               <  1:5695635/45‑1 (MQ=255)
                                 tGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtc                                      <  1:1586922/70‑1 (MQ=255)
                                 tGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGAGATACAGAGGAATTTGTATtc                                      <  1:1905209/70‑1 (MQ=255)
                                  ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc                                                           <  1:288495/48‑1 (MQ=255)
                                  ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc                                                           <  1:4533583/48‑1 (MQ=255)
                                  ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgt                                    >  1:422688/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)
                                      tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATAc                                                      >  1:3922263/1‑49 (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)
                                          gATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                           <  1:1502902/56‑1 (MQ=255)
                                          gATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                           <  1:1344612/56‑1 (MQ=255)
                                           gtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCTATc                             <  1:2453718/68‑1 (MQ=255)
                                           gtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATc                             <  1:184854/68‑1 (MQ=255)
                                           aTAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATc                             <  1:4116209/69‑1 (MQ=255)
                                           aTAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATc                             <  1:1020268/69‑1 (MQ=255)
                                           aTAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa                           <  1:3230301/71‑1 (MQ=255)
                                            tAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                              >  1:1283257/1‑67 (MQ=255)
                                            tAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAAc                          >  1:1370218/1‑71 (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)
                                                   aCGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACtcatcg                                  >  1:2748907/1‑56 (MQ=255)
                                                   aCGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATcta                           >  1:2418078/1‑61 (MQ=255)
                                                   aCGTCTTGTTTGCCGAACCACCACCTTTGGCGATACaga                                                   >  1:1900631/1‑39 (MQ=255)
                                                    cGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTaca                  <  1:4390952/71‑3 (MQ=255)
                                                     gTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                           <  1:5538074/45‑1 (MQ=255)
                                                     gTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                           <  1:1648357/45‑1 (MQ=255)
                                                     gTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                           <  1:5427927/45‑1 (MQ=255)
                                                       cTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAt                                              <  1:1160647/40‑1 (MQ=255)
                                                        ttGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTAt                                        <  1:1750528/45‑1 (MQ=255)
                                                        ttGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                    >  1:5291368/1‑65 (MQ=255)
                                                          gTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                           <  1:2019427/40‑1 (MQ=255)
                                                          gTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa                           >  1:4382572/1‑56 (MQ=255)
                                                          gTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCaa                           >  1:4162793/1‑56 (MQ=255)
                                                           tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                            <  1:3977140/38‑1 (MQ=255)
                                                           tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                            <  1:1592012/38‑1 (MQ=255)
                                                           tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                            <  1:1429455/38‑1 (MQ=255)
                                                           tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                            <  1:4485693/38‑1 (MQ=255)
                                                           tttGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcgat           >  1:4978529/1‑71 (MQ=255)
                                                           tttGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcgat           >  1:360859/1‑71 (MQ=255)
                                                             tGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGt                                          >  1:2529448/1‑38 (MQ=255)
                                                             tGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGt                                          >  1:4518594/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)
                                                                      |                                                                    
AATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGG  >  minE/1015438‑1015576

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: