Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,280 G→A 27.1% H263H (CAC→CAT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,2800GA27.1% 165.0 / 62.7 100H263H (CAC→CATfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (43/29);  new base A (14/13);  total (57/43)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTC  >  minE/1015213‑1015349
                                                                   |                                                                     
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)
  aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCt                                                                   >  1:5513014/1‑70 (MQ=255)
  aGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCt                                                                   >  1:3869904/1‑70 (MQ=255)
   gACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCa                                                                 <  1:2054361/71‑1 (MQ=255)
     ccAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCAtt                                                               <  1:3362463/71‑1 (MQ=255)
     ccAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCAtt                                                               >  1:7295/1‑71 (MQ=255)
      cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTc                                                              >  1:4611896/1‑71 (MQ=255)
      cAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTc                                                              >  1:1532803/1‑71 (MQ=255)
         gATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCt                                                           >  1:4712877/1‑71 (MQ=255)
           ttttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTc                                                         <  1:1593021/71‑1 (MQ=255)
            tttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                                                        >  1:3523973/1‑71 (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:5237853/1‑69 (MQ=255)
                   ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTgg                                                   >  1:5049742/1‑69 (MQ=255)
                     cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                               >  1:5597289/1‑71 (MQ=255)
                     cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                               >  1:825667/1‑71 (MQ=255)
                     cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                               >  1:841212/1‑71 (MQ=255)
                     cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                               >  1:1156936/1‑71 (MQ=255)
                      ttccagcagTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAg                                                <  1:2850847/66‑1 (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:31893/1‑54 (MQ=12)
                           gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                                                          >  1:4980568/1‑54 (MQ=12)
                           gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                                                          >  1:1138379/1‑54 (MQ=12)
                           gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                                                          >  1:5273397/1‑54 (MQ=12)
                           gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                                                          >  1:890307/1‑54 (MQ=12)
                           gCAGTTCCTGTTCCAGCTGGATATCGCGGAACGCCTGACCATGTTCGTTCCCtt                                                          >  1:4227379/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:2284907/69‑4 (MQ=255)
                                     ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATagtg                                 <  1:4404743/69‑4 (MQ=255)
                                     ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATagtg                                 <  1:5569910/69‑4 (MQ=255)
                                     ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATagt                                  <  1:4736138/68‑3 (MQ=255)
                                     ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATagt                                  <  1:4292728/68‑3 (MQ=255)
                                     ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATagt                                  <  1:4515649/68‑3 (MQ=255)
                                      tCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTcgc                              >  1:983451/1‑71 (MQ=255)
                                           ttCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                                                        <  1:1152013/40‑1 (MQ=255)
                                           ttCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                                                        <  1:3129600/40‑1 (MQ=255)
                                           ttCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                                                        <  1:34462/40‑1 (MQ=255)
                                           ttCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                                                        <  1:4043971/40‑1 (MQ=255)
                                           ttCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAg                         >  1:4799669/1‑71 (MQ=255)
                                             ccACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAg                                                >  1:3122701/1‑46 (MQ=255)
                                             ccACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcatca                                         >  1:5571181/1‑53 (MQ=255)
                                             ccACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGc                        >  1:404647/1‑70 (MQ=255)
                                             ccACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGc                        >  1:3446307/1‑70 (MQ=255)
                                                 aTCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAg                                      >  1:4489761/1‑52 (MQ=255)
                                                 aTCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGtt                   >  1:4543378/1‑71 (MQ=255)
                                                 aTCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtc                                             <  1:5339815/45‑1 (MQ=255)
                                                 aTCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtc                                             <  1:3609211/45‑1 (MQ=255)
                                                  tCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGt                                                      >  1:1574021/1‑35 (MQ=255)
                                                  tCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcat                                           <  1:4574071/46‑1 (MQ=255)
                                                  tCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGttt                  <  1:694299/71‑1 (MQ=255)
                                                  tCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGttt                  <  1:3077542/71‑1 (MQ=255)
                                                  tCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGttt                  <  1:1834510/71‑1 (MQ=255)
                                                  tCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGttt                  >  1:4235186/1‑71 (MQ=255)
                                                   cgcgGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGtt                                              >  1:476484/1‑42 (MQ=255)
                                                   cgcgGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGtt                                              >  1:524866/1‑42 (MQ=255)
                                                   cgcgGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGtt                                              >  1:751471/1‑42 (MQ=255)
                                                   cgcgGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGtt                                              >  1:4184905/1‑42 (MQ=255)
                                                   cgcgGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGt                                               >  1:4101945/1‑41 (MQ=25)
                                                   cgcgGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtcatc                                          >  1:2080178/1‑46 (MQ=255)
                                                   cgcgGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtcat                                           >  1:5415063/1‑45 (MQ=255)
                                                    gcgGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGaa                          >  1:2008294/1‑61 (MQ=255)
                                                    gcgGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTa                                       >  1:5257797/1‑48 (MQ=255)
                                                    gcgGAACGCCTGACCATGTTCGTTACCTTCCGTCGGc                                                  >  1:5326464/1‑37 (MQ=16)
                                                    gcgGAACGCCCGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTa                                       >  1:3393111/1‑48 (MQ=19)
                                                     cgGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCAc               <  1:2302416/71‑1 (MQ=255)
                                                     cgGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtcat                                           <  1:5137597/43‑1 (MQ=25)
                                                      ggAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAg                                      <  1:2423036/47‑1 (MQ=255)
                                                        aaCGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTa                                    >  1:2587117/1‑47 (MQ=255)
                                                        aaCGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTa                                    >  1:1589301/1‑47 (MQ=255)
                                                        aaCGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTa                                    >  1:5109643/1‑47 (MQ=255)
                                                           gCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGtttt         <  1:5389410/71‑1 (MQ=255)
                                                            ccTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGtttt         >  1:2418531/1‑70 (MQ=255)
                                                            ccTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGtttt         <  1:4389633/70‑1 (MQ=255)
                                                            ccTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGtttt         >  1:1930855/1‑70 (MQ=255)
                                                                aCCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGt    <  1:3751543/71‑1 (MQ=255)
                                                                  cGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTc  <  1:1029168/71‑1 (MQ=255)
                                                                   |                                                                     
CCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTC  >  minE/1015213‑1015349

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: