Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,256 T→C 21.5% E271E (GAA→GAG) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,2560TC21.5% 182.7 / 44.5 93E271E (GAA→GAG) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (47/26);  new base C (10/10);  total (57/36)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

GTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAG  >  minE/1015190‑1015326
                                                                  |                                                                      
gTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                                    >  1:5691118/1‑71 (MQ=255)
 tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTc                                                                      <  1:5335744/68‑1 (MQ=255)
 tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTc                                                                      <  1:2062810/68‑1 (MQ=255)
 tATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTc                                                                      <  1:170607/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)
     tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcg                                                               >  1:5206192/1‑71 (MQ=255)
     tACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATcg                                                               >  1:4518929/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)
             aaCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGtt                                                                       >  1:2185789/1‑55 (MQ=255)
             aaCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGc                                                       >  1:3749115/1‑71 (MQ=255)
              aCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                                    <  1:1751104/57‑1 (MQ=255)
              aCTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGcc                                                      >  1:310655/1‑71 (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:4095523/1‑71 (MQ=255)
                       ccAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGc                                             >  1:5649973/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:1661203/1‑46 (MQ=255)
                         aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCca                                                                    >  1:1319329/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)
                             cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTc                                       >  1:1532803/1‑71 (MQ=255)
                             cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTc                                       >  1:4611896/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:3626270/1‑53 (MQ=255)
                                    ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGAc                                                  >  1:2753235/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:3321426/56‑1 (MQ=255)
                                          ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCAt                                         <  1:4761850/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:1156936/1‑71 (MQ=255)
                                            cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                        >  1:5597289/1‑71 (MQ=255)
                                            cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                        >  1:841212/1‑71 (MQ=255)
                                            cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                        >  1:825667/1‑71 (MQ=255)
                                             ttccagcagTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAg                         <  1:2850847/66‑1 (MQ=255)
                                              aTCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGa                                                   >  1:3672237/1‑42 (MQ=255)
                                              aTCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGa                                                   >  1:3801326/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:1069328/1‑70 (MQ=255)
                                                 aGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcat                    >  1:2260021/1‑70 (MQ=255)
                                                 aGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcat                    <  1:3075955/70‑1 (MQ=255)
                                                 aGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCCTTGGCAGTtcat                    <  1:3081069/70‑1 (MQ=255)
                                                  gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                                   >  1:4227379/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:4980568/1‑54 (MQ=12)
                                                  gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                                   >  1:31893/1‑54 (MQ=12)
                                                  gCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtcatcg                  <  1:752292/71‑2 (MQ=255)
                                                    aaTTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcatca                  >  1:2610039/1‑69 (MQ=255)
                                                      ttCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCtt                                   <  1:5444176/50‑1 (MQ=255)
                                                      ttCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCtt                                   <  1:3602112/50‑1 (MQ=255)
                                                        cTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTAt            >  1:750885/1‑71 (MQ=255)
                                                          ttttCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTAtt           >  1:3862745/1‑70 (MQ=255)
                                                            ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTTATCGTAATagt           <  1:17031/68‑3 (MQ=255)
                                                            ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTaa               >  1:3190261/1‑63 (MQ=255)
                                                            ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATagtg          <  1:5569910/69‑4 (MQ=255)
                                                            ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATagtg          <  1:4404743/69‑4 (MQ=255)
                                                            ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATagtg          <  1:2284907/69‑4 (MQ=255)
                                                            ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATagt           <  1:4736138/68‑3 (MQ=255)
                                                            ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATagt           <  1:4515649/68‑3 (MQ=255)
                                                            ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATagt           <  1:4292728/68‑3 (MQ=255)
                                                             tCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTcgc       >  1:983451/1‑71 (MQ=255)
                                                                  ttCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                                 <  1:4043971/40‑1 (MQ=255)
                                                                  ttCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                                 <  1:34462/40‑1 (MQ=255)
                                                                  ttCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                                 <  1:3129600/40‑1 (MQ=255)
                                                                  ttCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                                 <  1:1152013/40‑1 (MQ=255)
                                                                  ttCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAg  >  1:4799669/1‑71 (MQ=255)
                                                                  |                                                                      
GTATTTACCACCAAACTGCGCACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAG  >  minE/1015190‑1015326

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: