Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,520 G→A 22.6% I183I (ATC→ATT) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,5200GA22.6% 216.5 / 58.1 115I183I (ATC→ATT) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (44/45);  new base A (12/14);  total (56/59)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 119 total aligned reads displayed.

CCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGT  >  minE/1015452‑1015585
                                                                    |                                                                 
cccGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGa                                                                  <  1:1886545/70‑1 (MQ=255)
cccGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGa                                                                  <  1:862705/70‑1 (MQ=255)
         aGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGaa                                                        >  1:3129290/1‑71 (MQ=255)
         aGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGaa                                                        >  1:5405961/1‑71 (MQ=255)
          gTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAt                                                       <  1:3288018/71‑1 (MQ=255)
             gggCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTg                                                    >  1:3761804/3‑71 (MQ=255)
             gggCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTg                                                    <  1:289407/69‑1 (MQ=255)
             gggCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTg                                                    >  1:4443382/3‑71 (MQ=255)
             aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGa                                                                  >  1:4311400/1‑57 (MQ=255)
             aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    >  1:2861294/1‑71 (MQ=255)
             aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    <  1:2210864/71‑1 (MQ=255)
             aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    >  1:5656091/1‑71 (MQ=255)
              cGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGt                                                   <  1:523433/71‑1 (MQ=255)
               gCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTa                                                  <  1:2198266/71‑1 (MQ=255)
               gCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTa                                                  <  1:1073069/71‑1 (MQ=255)
                cTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTac                                                 >  1:4437263/1‑70 (MQ=255)
                cTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTac                                                 >  1:279029/1‑70 (MQ=255)
                cTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTAt                                                 <  1:2561958/71‑1 (MQ=255)
                 tttGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTg                                                    <  1:2037477/67‑1 (MQ=255)
                 tttGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTg                                                    <  1:4939290/67‑1 (MQ=255)
                   tGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACtc                                               <  1:4594302/70‑1 (MQ=255)
                   tGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACtc                                               <  1:1224949/70‑1 (MQ=255)
                   tGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtc                                               <  1:1586922/70‑1 (MQ=255)
                    ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgt                                             >  1:3400818/1‑71 (MQ=255)
                    ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgt                                             >  1:422688/1‑71 (MQ=255)
                    ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgt                                             >  1:5297930/1‑71 (MQ=255)
                    ggTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcg                                              >  1:1547230/1‑70 (MQ=255)
                     gTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgt                                             <  1:4363897/70‑1 (MQ=255)
                      tttCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGaa                                                        >  1:5251844/1‑58 (MQ=255)
                        tCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGc                                          >  1:2245741/1‑70 (MQ=255)
                        tCCTGGTAGAGATACGTTTTGTTGGCAGAGCCCCGCCTTTCGCAACGCAAAGGaa                                                        <  1:1819799/55‑1 (MQ=255)
                        tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATAc                                                               >  1:3922263/1‑49 (MQ=255)
                        tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATtcgtc                                            <  1:3839973/68‑1 (MQ=255)
                        tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGc                                          <  1:4641737/70‑1 (MQ=255)
                        tCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGc                                          >  1:2995995/1‑70 (MQ=255)
                        tCCTGATAGAGATACGTCTTGTTAGCCGAACCACCACCTTTGGCGa                                                                  >  1:638196/1‑46 (MQ=255)
                          cTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGt                                                   <  1:262458/59‑1 (MQ=255)
                          cTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGt                                                   <  1:4192742/59‑1 (MQ=255)
                           tggtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                                       >  1:3978994/4‑70 (MQ=255)
                           tggtAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                                       >  1:1893027/4‑70 (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:1750586/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)
                             aTAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTAGTCGCCATc                                      <  1:4432689/69‑1 (MQ=255)
                              tAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCAt                                       >  1:1283257/1‑67 (MQ=255)
                              tAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAAc                                   >  1:1370218/1‑71 (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:2055945/1‑39 (MQ=255)
                                     aCGTCTTGTTTGCCGAACCACCACCTTTGGCGATACaga                                                            >  1:1900631/1‑39 (MQ=255)
                                      cGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTaca                           <  1:4390952/71‑3 (MQ=255)
                                       gTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    <  1:1648357/45‑1 (MQ=255)
                                       gTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    <  1:5427927/45‑1 (MQ=255)
                                       gTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    <  1:5538074/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)
                                          ttGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAga                         <  1:1196073/69‑1 (MQ=255)
                                            gTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg                                                    <  1:2019427/40‑1 (MQ=255)
                                             tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                                     <  1:3977140/38‑1 (MQ=255)
                                             tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                                     <  1:1592012/38‑1 (MQ=255)
                                             tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                                     <  1:3911620/38‑1 (MQ=255)
                                             tttGCCGAACCACCACCTTTGGCGATACAGAGGAAttt                                                     <  1:1429455/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:4518594/1‑38 (MQ=255)
                                               tGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGt                                                   >  1:2529448/1‑38 (MQ=255)
                                               tGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGt                                                   >  1:3242740/1‑38 (MQ=255)
                                                  cGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTgcg               <  1:4692434/71‑1 (MQ=255)
                                                    aaCCACCACCTTTGGCGATACAGAGGAATTTGTATtcg                                              >  1:3157211/1‑38 (MQ=255)
                                                    aaCCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgatcgatc                   <  1:743428/65‑1 (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:4584255/47‑1 (MQ=255)
                                                            ccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                             <  1:4458752/47‑1 (MQ=255)
                                                            ccTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGAtt     <  1:4524809/71‑1 (MQ=255)
                                                             cTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAgat                        >  1:3159724/1‑51 (MQ=255)
                                                             cTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTg    <  1:3775936/71‑1 (MQ=255)
                                                              tttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTgg   >  1:2072133/1‑71 (MQ=255)
                                                              tttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTgg   >  1:3325390/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:4099552/1‑71 (MQ=255)
                                                               ttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGt  >  1:1891320/1‑71 (MQ=255)
                                                               ttGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGt  >  1:226183/1‑71 (MQ=255)
                                                                    gATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                             <  1:5315601/39‑1 (MQ=255)
                                                                    gATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                             <  1:4176970/39‑1 (MQ=255)
                                                                    gATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                             <  1:4971077/39‑1 (MQ=255)
                                                                    gATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATa                             <  1:2108055/39‑1 (MQ=255)
                                                                    |                                                                 
CCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGT  >  minE/1015452‑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: