Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,843,402 A→G 24.4% T193T (ACT→ACC frdA ← fumarate reductase (anaerobic) catalytic and NAD/flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,843,4020AG24.4% 133.6 / 39.8 110T193T (ACT→ACCfrdAfumarate reductase (anaerobic) catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (31/52);  new base G (14/13);  total (45/65)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 110 total aligned reads displayed.

CCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCT  >  minE/2843337‑2843465
                                                                  |                                                               
cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTa                                                              >  1:1027206/1‑69 (MQ=255)
cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTa                                                              >  1:2728213/1‑69 (MQ=255)
cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTa                                                              >  1:257956/1‑69 (MQ=255)
cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTa                                                              >  1:255225/1‑69 (MQ=255)
cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTa                                                              >  1:2459908/1‑69 (MQ=255)
cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTa                                                              >  1:2339178/1‑69 (MQ=255)
cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTa                                                              >  1:3191121/1‑69 (MQ=255)
cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTa                                                              >  1:3239020/1‑69 (MQ=255)
cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTa                                                              >  1:3380228/1‑69 (MQ=255)
cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTa                                                              >  1:1524127/1‑69 (MQ=255)
cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTa                                                              >  1:827466/1‑69 (MQ=255)
cccataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTa                                                              >  1:1151028/1‑69 (MQ=255)
  cataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGc                                                            <  1:2851774/69‑1 (MQ=255)
  cataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGc                                                            <  1:2396848/69‑1 (MQ=255)
   ataccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGc                                                            <  1:182064/68‑1 (MQ=255)
      ccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATa                                                        >  1:404646/1‑69 (MQ=255)
      ccGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATa                                                        >  1:2153840/1‑69 (MQ=255)
        gTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCGGTAGCCATAAc                                                      <  1:367216/69‑1 (MQ=255)
        gTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCGGTAGCCATAAc                                                      <  1:2913374/69‑1 (MQ=255)
        gTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCGGTAGCCATAAc                                                      <  1:1844976/69‑1 (MQ=255)
        gTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCGGTAGCCATAAc                                                      <  1:3324624/69‑1 (MQ=255)
         tCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCGGTAGCCATAAc                                                      <  1:1987050/68‑1 (MQ=255)
         tCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCGGTAGCCATAAc                                                      <  1:3002947/68‑1 (MQ=255)
         tCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCGGTAGCCATAAc                                                      <  1:151747/68‑1 (MQ=255)
         tCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACg                                                     <  1:2249737/69‑1 (MQ=255)
           aCCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGAc                                                   <  1:2948493/69‑1 (MQ=255)
             cGGAAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCGGTAGCCATAAc                                                      <  1:3192666/64‑1 (MQ=255)
              ggTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCGGTAGCCATAAc                                                      <  1:653283/63‑1 (MQ=255)
                tAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCGGTAGCCATAAc                                                      <  1:153655/61‑1 (MQ=255)
                    gATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                           <  1:2728394/68‑1 (MQ=255)
                    gATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                           <  1:750158/68‑1 (MQ=255)
                    gATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                           <  1:1925564/68‑1 (MQ=255)
                    gATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                           <  1:1944830/68‑1 (MQ=255)
                    gATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                           <  1:1427526/68‑1 (MQ=255)
                     attccgccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCTTTAGCa                                         >  1:455898/4‑69 (MQ=255)
                     aTGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGc                                          >  1:1791022/1‑68 (MQ=255)
                     aTGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                         >  1:1643800/1‑69 (MQ=255)
                     aTGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                         >  1:3306973/1‑69 (MQ=255)
                     aTGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                         >  1:2852451/1‑69 (MQ=255)
                     aTGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                         >  1:2337413/1‑69 (MQ=255)
                     aTGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                         >  1:1084972/1‑69 (MQ=255)
                     aTGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                         >  1:2347999/1‑69 (MQ=255)
                     aTGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                         >  1:2577727/1‑69 (MQ=255)
                     aTGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCa                                         >  1:3372847/1‑69 (MQ=255)
                       gccgccGTTGGTGTTGTAACGATAACCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACg                                       >  1:1986393/1‑69 (MQ=255)
                       gccgccGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCAc                                        >  1:3234724/1‑68 (MQ=255)
                        ccgccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCGGTAGCCATAAc                                                      <  1:2612263/53‑1 (MQ=255)
                        ccgccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACgg                                      <  1:3419797/69‑1 (MQ=255)
                        ccgccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACgg                                      <  1:2863537/69‑1 (MQ=255)
                        ccgccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACgg                                      <  1:1839929/69‑1 (MQ=255)
                        ccgccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACgg                                      <  1:2237780/69‑1 (MQ=255)
                        ccgccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACgg                                      <  1:2165756/69‑1 (MQ=255)
                          gccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGa                                     <  1:432101/68‑1 (MQ=255)
                          gccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGa                                     <  1:3408959/68‑1 (MQ=255)
                           ccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATc                                   <  1:892806/69‑1 (MQ=255)
                           ccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATc                                   <  1:2620058/69‑1 (MQ=255)
                           ccgTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATc                                   <  1:1462154/69‑1 (MQ=255)
                              ttGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGc                                <  1:1071698/69‑1 (MQ=255)
                              ttGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGc                                <  1:2668616/69‑1 (MQ=255)
                              ttGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGc                                <  1:1673657/69‑1 (MQ=255)
                                 gtgtTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACg                                                     <  1:2783658/45‑1 (MQ=255)
                                 gtgtTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACg                                                     <  1:731534/45‑1 (MQ=255)
                                 gtgtTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACg                                                     <  1:1455206/45‑1 (MQ=255)
                                 gtgtTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACg                                                     <  1:1461940/45‑1 (MQ=255)
                                 gtgtTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACg                                                     <  1:1655983/45‑1 (MQ=255)
                                    ttgtAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACgg                                      <  1:2105846/57‑1 (MQ=255)
                                      gtAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGc                                <  1:2688185/61‑1 (MQ=255)
                                       tAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGc                                <  1:1300803/60‑1 (MQ=255)
                                         aCGATAAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                           >  1:1093251/1‑48 (MQ=255)
                                         aCGATAAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGc                      >  1:907771/1‑69 (MQ=255)
                                         aCGATAAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGc                      >  1:1816521/1‑69 (MQ=255)
                                         aCGATAAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGc                      >  1:1938434/1‑69 (MQ=255)
                                          cGATAAACGCGACCCGCACCGCCGGTAGCCATAAc                                                      <  1:1237576/35‑1 (MQ=255)
                                          cGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGa                                     <  1:2273410/52‑1 (MQ=255)
                                           gATAAACGCGACCCGCACCGCCAGTAGCCATAACGAc                                                   <  1:2614319/37‑1 (MQ=255)
                                           gATAAACGCGACCCGCACCGCCAGTAGCCATAACGAc                                                   <  1:505846/37‑1 (MQ=255)
                                                aaCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                           <  1:2323828/41‑1 (MQ=255)
                                                 aCGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                           <  1:3072662/40‑1 (MQ=255)
                                                  cgcgACCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                           <  1:1767565/39‑1 (MQ=255)
                                                  cgcgACCCGCACCGCCAGTAGCCATAACGACCGCGTTAg                                           <  1:1222083/39‑1 (MQ=255)
                                                   gcgACCCGCACCGCCGGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat            >  1:1681523/1‑69 (MQ=255)
                                                   gcgACCCGCACCGCCGGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat            >  1:3062752/1‑69 (MQ=255)
                                                   gcgACCCGCACCGCCGGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat            >  1:2855949/1‑69 (MQ=255)
                                                   gcgACCCGCACCGCCGGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat            >  1:1139926/1‑69 (MQ=255)
                                                   gcgACCCGCACCGCCGGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat            >  1:2567989/1‑69 (MQ=255)
                                                   gcgACCCGCACCGCCGGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat            >  1:302979/1‑69 (MQ=255)
                                                   gcgACCCGCACCGCCGGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat            >  1:1984131/1‑69 (MQ=255)
                                                   gcgACCCGCACCGCCGGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat            >  1:35625/1‑69 (MQ=255)
                                                   gcgACCCGCACCGCCGGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat            >  1:3400158/1‑69 (MQ=255)
                                                   gcgACCCGCACCGCCGGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCcatcat            >  1:3007561/1‑69 (MQ=255)
                                                      aCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGtt         <  1:2933020/69‑1 (MQ=255)
                                                      aCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGtt         <  1:2068950/69‑1 (MQ=255)
                                                      aCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGtt         <  1:511444/69‑1 (MQ=255)
                                                      aCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGtt         <  1:650428/69‑1 (MQ=255)
                                                      aCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGtt         <  1:1607260/69‑1 (MQ=255)
                                                      aCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGtt         <  1:72948/69‑1 (MQ=255)
                                                      aCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGtt         <  1:1144069/69‑1 (MQ=255)
                                                      aCCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGtt         <  1:1670852/69‑1 (MQ=255)
                                                       cccGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGtt         <  1:3150201/68‑1 (MQ=255)
                                                              cGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCt  <  1:2986784/68‑1 (MQ=255)
                                                                  |                                                               
CCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCAGTAGCCATAACGACCGCGTTAGCACGGATCTGCACCAGCGTGCCTTCCATCATGTTCATTGCT  >  minE/2843337‑2843465

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: