Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,247 T→C 14.7% K274K (AAA→AAG) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,2470TC14.7% 232.8 / 27.3 95K274K (AAA→AAG) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (52/29);  new base C (10/4);  total (62/33)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTAT  >  minE/1015177‑1015316
                                                                      |                                                                     
tGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCt                                                                       <  1:2910283/71‑1 (MQ=255)
       aGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTcc                                                                 <  1:89427/70‑1 (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:4026848/1‑62 (MQ=255)
            aGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCtttt                                                                    >  1:2533073/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:5475389/71‑1 (MQ=255)
                                     cAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCt                                  >  1:4954005/1‑71 (MQ=255)
                                      aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                          >  1:1888053/1‑46 (MQ=255)
                                      aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                          >  1:1319329/1‑46 (MQ=255)
                                      aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                          >  1:1661203/1‑46 (MQ=255)
                                      aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCt                                  >  1:5513014/1‑70 (MQ=255)
                                      aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCt                                  >  1:3869904/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:4994845/1‑35 (MQ=255)
                                          cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTcc                                                                 >  1:389368/1‑35 (MQ=255)
                                          cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTcc                                                                 >  1:3391593/1‑35 (MQ=255)
                                          cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTcc                                                                 >  1:2762077/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:5560640/1‑53 (MQ=255)
                                                 ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGAc                                        >  1:2753235/1‑53 (MQ=255)
                                                 ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGAc                                        >  1:3626270/1‑53 (MQ=255)
                                                 ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                       <  1:4880941/70‑1 (MQ=255)
                                                 ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                       >  1:2950151/1‑70 (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:4761850/56‑1 (MQ=255)
                                                       ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCAt                               <  1:3321426/56‑1 (MQ=255)
                                                       ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCAt                               <  1:1211112/56‑1 (MQ=255)
                                                       ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTgg                  >  1:5049742/1‑69 (MQ=255)
                                                       ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTgg                  >  1:5237853/1‑69 (MQ=255)
                                                         cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt              >  1:841212/1‑71 (MQ=255)
                                                         cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt              >  1:825667/1‑71 (MQ=255)
                                                         cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt              >  1:5597289/1‑71 (MQ=255)
                                                         cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt              >  1:1156936/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:1069328/1‑70 (MQ=255)
                                                              aGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcat          <  1:3075955/70‑1 (MQ=255)
                                                              aGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCCTTGGCAGTtcat          <  1:3081069/70‑1 (MQ=255)
                                                               gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                         >  1:31893/1‑54 (MQ=12)
                                                               gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                         >  1:1138379/1‑54 (MQ=12)
                                                               gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                         >  1:890307/1‑54 (MQ=12)
                                                               gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                         >  1:5273397/1‑54 (MQ=12)
                                                               gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                         >  1:4227379/1‑54 (MQ=12)
                                                               gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                         >  1:4980568/1‑54 (MQ=12)
                                                               gCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtcatcg        <  1:752292/71‑2 (MQ=255)
                                                                 aaTTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcatca        >  1:2610039/1‑69 (MQ=255)
                                                                   ttCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCtt                         <  1:3602112/50‑1 (MQ=255)
                                                                   ttCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCtt                         <  1:5444176/50‑1 (MQ=255)
                                                                     cTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTAt  >  1:750885/1‑71 (MQ=255)
                                                                      |                                                                     
TGTCGTGAGCGAAGTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTAT  >  minE/1015177‑1015316

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: