Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,243 A→G 14.4% L276L (TTG→CTG)  fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,2430AG14.4% 239.6 / 26.7 97L276L (TTG→CTG) fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (54/29);  new base G (11/3);  total (65/32)
Fisher's exact test for biased strand distribution p-value = 3.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 101 total aligned reads displayed.

GGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCA  >  minE/1015174‑1015310
                                                                     |                                                                   
ggATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAAt                                                                    >  1:2359578/1‑71 (MQ=255)
  aTGTCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTc                                                                  <  1:2320369/71‑1 (MQ=255)
   tGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCt                                                                 <  1:2910283/71‑1 (MQ=255)
     tCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTcctg                                                               >  1:2990150/1‑68 (MQ=255)
     tCGTGCGCGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTcctg                                                               >  1:5068033/1‑68 (MQ=255)
           gCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAg                                                         <  1:5613530/71‑1 (MQ=255)
            cGAAGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAgc                                                        >  1:2203378/1‑70 (MQ=255)
               aGTATTTACCGCCAAACTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAgctg                                                      <  1:4441711/70‑4 (MQ=255)
               aGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCtttt                                                              >  1:2533073/1‑62 (MQ=255)
               aGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCtttt                                                              >  1:4026848/1‑62 (MQ=255)
                gTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                    >  1:5691118/1‑71 (MQ=255)
                 tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTc                                                      <  1:2062810/68‑1 (MQ=255)
                 tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTc                                                      <  1:170607/68‑1 (MQ=255)
                 tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTc                                                      <  1:5335744/68‑1 (MQ=255)
                 tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCcac                                                   <  1:4015888/71‑1 (MQ=255)
                  aTTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCcac                                                   >  1:2587251/1‑70 (MQ=255)
                     tACCACCAAACTGCGCCCCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcg                                               >  1:4337732/1‑71 (MQ=255)
                     tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCttt                                                               >  1:2086504/1‑55 (MQ=255)
                     tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAg                                                         >  1:4687202/1‑61 (MQ=255)
                     tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcg                                               >  1:4518929/1‑71 (MQ=255)
                     tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcg                                               >  1:5206192/1‑71 (MQ=255)
                        caccaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCgg                                            <  1:4782366/71‑1 (MQ=255)
                          ccaAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGaa                                          >  1:4364390/1‑71 (MQ=255)
                           caAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAAc                                         >  1:2584735/1‑71 (MQ=255)
                            aaaCTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACg                                        >  1:3177950/1‑71 (MQ=255)
                             aaCTGCGCGCCAAGACCGAGTTTCTGGGCCTCTTCCAGCAGTTCCTGTTCCAGCTGGACATCGCGGAACGc                                       >  1:5115266/1‑71 (MQ=255)
                             aaCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCtttt                                                              <  1:3385621/48‑1 (MQ=255)
                             aaCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGtt                                                       >  1:2185789/1‑55 (MQ=255)
                             aaCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGc                                       >  1:3749115/1‑71 (MQ=255)
                              aCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTc                                                            <  1:2063972/49‑1 (MQ=255)
                              aCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                    <  1:1751104/57‑1 (MQ=255)
                              aCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGcc                                      >  1:310655/1‑71 (MQ=255)
                               cTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTc                                                            >  1:4284153/1‑48 (MQ=255)
                                  cgcACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTg                                    <  1:1220813/69‑1 (MQ=255)
                                    cACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGaa                                          >  1:1358784/1‑61 (MQ=255)
                                      cccAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                    <  1:2871839/49‑1 (MQ=255)
                                      cccAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                    <  1:4259719/49‑1 (MQ=255)
                                       ccAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                    <  1:4973735/48‑1 (MQ=255)
                                       ccAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGc                             >  1:5649973/1‑71 (MQ=255)
                                       ccAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGc                             >  1:4095523/1‑71 (MQ=255)
                                        cAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCt                            >  1:4954005/1‑71 (MQ=255)
                                        cAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCt                            <  1:5475389/71‑1 (MQ=255)
                                         aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                    >  1:1661203/1‑46 (MQ=255)
                                         aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                    >  1:1888053/1‑46 (MQ=255)
                                         aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                    >  1:1319329/1‑46 (MQ=255)
                                         aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCt                            >  1:3869904/1‑70 (MQ=255)
                                         aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCt                            >  1:5513014/1‑70 (MQ=255)
                                          gACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCa                          <  1:2054361/71‑1 (MQ=255)
                                            ccAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACAt                                                 >  1:1683929/1‑46 (MQ=255)
                                            ccAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCAtt                        <  1:3362463/71‑1 (MQ=255)
                                            ccAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCAtt                        >  1:7295/1‑71 (MQ=255)
                                             cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTcc                                                           >  1:3391593/1‑35 (MQ=255)
                                             cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTcc                                                           >  1:2762077/1‑35 (MQ=255)
                                             cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTcc                                                           >  1:389368/1‑35 (MQ=255)
                                             cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTcc                                                           >  1:4994845/1‑35 (MQ=255)
                                             cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTc                       >  1:4611896/1‑71 (MQ=255)
                                             cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTc                       >  1:1532803/1‑71 (MQ=255)
                                             cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCAAGTTCCACATcg                                               >  1:4392533/1‑47 (MQ=255)
                                                gATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcgcg                                             >  1:4825571/1‑46 (MQ=255)
                                                gATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcgc                                              >  1:2410770/1‑45 (MQ=255)
                                                gATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCt                    >  1:4712877/1‑71 (MQ=255)
                                                  ttttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTc                  <  1:1593021/71‑1 (MQ=255)
                                                   tttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                 >  1:3523973/1‑71 (MQ=255)
                                                    ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcgcg                                             >  1:4311236/1‑42 (MQ=255)
                                                    ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGAc                                  >  1:2753235/1‑53 (MQ=255)
                                                    ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGAc                                  >  1:3626270/1‑53 (MQ=255)
                                                    ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGAc                                  >  1:5560640/1‑53 (MQ=255)
                                                    ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                 >  1:2950151/1‑70 (MQ=255)
                                                    ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                 <  1:4880941/70‑1 (MQ=255)
                                                    ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCg                >  1:1366151/1‑71 (MQ=255)
                                                        gcTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCa                          <  1:538864/57‑1 (MQ=255)
                                                          ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCTTGCTCAt                         <  1:3681184/56‑1 (MQ=255)
                                                          ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCt                            <  1:1386262/53‑1 (MQ=255)
                                                          ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCAt                         <  1:1211112/56‑1 (MQ=255)
                                                          ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCAt                         <  1:4761850/56‑1 (MQ=255)
                                                          ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCAt                         <  1:3321426/56‑1 (MQ=255)
                                                          ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTgg            >  1:5237853/1‑69 (MQ=255)
                                                          ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTgg            >  1:5049742/1‑69 (MQ=255)
                                                            cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt        >  1:825667/1‑71 (MQ=255)
                                                            cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt        >  1:5597289/1‑71 (MQ=255)
                                                            cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt        >  1:1156936/1‑71 (MQ=255)
                                                            cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt        >  1:841212/1‑71 (MQ=255)
                                                             ttccagcagTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAg         <  1:2850847/66‑1 (MQ=255)
                                                              aTCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGa                                   >  1:3801326/1‑42 (MQ=255)
                                                              aTCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGa                                   >  1:3672237/1‑42 (MQ=255)
                                                               ccagcagTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTc                           >  1:659490/2‑49 (MQ=255)
                                                                cagcagTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtca     <  1:3529388/70‑1 (MQ=255)
                                                                 aGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGTCAGTtcat    <  1:462331/70‑1 (MQ=255)
                                                                 aGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcat    >  1:2260021/1‑70 (MQ=255)
                                                                 aGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcat    <  1:3075955/70‑1 (MQ=255)
                                                                 aGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcat    >  1:1069328/1‑70 (MQ=255)
                                                                 aGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCCTTGGCAGTtcat    <  1:3081069/70‑1 (MQ=255)
                                                                  gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                   >  1:890307/1‑54 (MQ=12)
                                                                  gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                   >  1:5273397/1‑54 (MQ=12)
                                                                  gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                   >  1:4980568/1‑54 (MQ=12)
                                                                  gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                   >  1:4227379/1‑54 (MQ=12)
                                                                  gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                   >  1:31893/1‑54 (MQ=12)
                                                                  gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                   >  1:1138379/1‑54 (MQ=12)
                                                                  gCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtcatcg  <  1:752292/71‑2 (MQ=255)
                                                                    aaTTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcatca  >  1:2610039/1‑69 (MQ=255)
                                                                     |                                                                   
GGATGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCA  >  minE/1015174‑1015310

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: