Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,283 C→T 11.6% E262E (GAG→GAA fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,2830CT11.6% 163.3 / 10.1 60E262E (GAG→GAAfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (23/30);  new base T (2/5);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTT  >  minE/1015213‑1015353
                                                                      |                                                                      
ccAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGc                                                                        >  1:1908612/1‑71 (MQ=255)
        aGATTTTGCGCTTCGATCAGCAATTCTTTTTACAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTccc                                                                >  1:3387117/1‑71 (MQ=255)
           ttttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTc                                                             >  1:1230565/1‑71 (MQ=255)
             ttGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTcc                                                            <  1:691656/70‑1 (MQ=255)
              tGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCg                                                           <  1:1863658/70‑1 (MQ=255)
              tGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCg                                                           <  1:3796645/70‑1 (MQ=255)
                cgcTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGtt                                                         <  1:3467956/70‑1 (MQ=255)
                 gcTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCg                                                           <  1:3005636/67‑1 (MQ=255)
                 gcTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTgg                                                       >  1:1797452/1‑71 (MQ=255)
                  cTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCa                                                                     >  1:2385587/1‑56 (MQ=255)
                  cTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCa                                                                     >  1:137271/1‑56 (MQ=255)
                  cTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCa                                                                     >  1:395392/1‑56 (MQ=255)
                  cTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGc                                                      <  1:279154/71‑1 (MQ=255)
                   ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTgg                                                       >  1:262112/1‑69 (MQ=255)
                     cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                                   >  1:1567985/1‑71 (MQ=255)
                     cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                                   >  1:1635787/1‑71 (MQ=255)
                     cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                                   >  1:638593/1‑71 (MQ=255)
                      gATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                                   <  1:1912878/70‑1 (MQ=255)
                      gATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                                   <  1:1938763/70‑1 (MQ=255)
                       aTCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTgg                                                       <  1:479704/65‑1 (MQ=255)
                         cagcagTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtcat                                               <  1:646300/71‑1 (MQ=255)
                         cAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcat                                               >  1:1757373/1‑71 (MQ=255)
                         cAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtca                                                >  1:321642/1‑70 (MQ=255)
                           gCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtcatcg                                             <  1:3297527/71‑2 (MQ=255)
                           gCAGTTCCTGTTCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtca                                                <  1:3419707/68‑1 (MQ=255)
                           gCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtcatca                                             <  1:1664489/71‑1 (MQ=255)
                            cAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTcatcat                                            <  1:290090/71‑1 (MQ=255)
                               ttCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCAtt                                                                   >  1:805861/1‑45 (MQ=255)
                                tCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCtt                                                              >  1:3218016/1‑49 (MQ=255)
                                tCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTtca                                                >  1:1262810/1‑63 (MQ=255)
                                 cTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAg                                          >  1:3740127/1‑68 (MQ=255)
                                  tttttCCAGTTCCACATCGCGGAACGCCTGACCTTGCTCATTCCCTTCCGtt                                                         <  1:3130134/52‑1 (MQ=255)
                                   ttttCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTAttt                                     >  1:2936579/1‑71 (MQ=255)
                                    tttCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGTCAGTTCATCATAGTAt                                       <  1:3444093/68‑1 (MQ=255)
                                    gttCCAGCTGGACATCGCGGAACGCCTGACCATGTTc                                                                      >  1:2890939/2‑37 (MQ=12)
                                     ttCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTc                                                             >  1:3805777/1‑45 (MQ=255)
                                     ttCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTc                                    <  1:3471748/70‑1 (MQ=255)
                                     ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTtcat                                               >  1:2544471/1‑59 (MQ=255)
                                     ttCCAGCTGGACATCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATagtg                                     <  1:1632628/69‑4 (MQ=255)
                                      tCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTcg                                   <  1:142764/70‑1 (MQ=255)
                                      tCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTcg                                   <  1:3540392/70‑1 (MQ=255)
                                       ccAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCGGTTg                                                        <  1:1109231/48‑1 (MQ=255)
                                       ccAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTg                                                        <  1:383607/48‑1 (MQ=255)
                                       ccAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTcg                                   <  1:1381048/69‑1 (MQ=255)
                                        cAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCgg                                <  1:1792455/71‑1 (MQ=255)
                                        cAGTTCCACATCGCGGAACGCCTGACCGTCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCgg                                <  1:3403505/70‑1 (MQ=255)
                                             ccACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGc                            >  1:63942/1‑70 (MQ=255)
                                                caTCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGt                        <  1:1568943/71‑1 (MQ=255)
                                                  tCGCGGAACGCCTGACCATGTTCGTTCCCTTCCGtc                                                         <  1:573673/36‑2 (MQ=25)
                                                       gAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCg                 >  1:1320581/1‑71 (MQ=255)
                                                         aCGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGaa                              >  1:1699660/1‑56 (MQ=255)
                                                          cGCCTGACCGTTCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGttt              <  1:2229574/71‑1 (MQ=255)
                                                          cGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTcgcg                                 <  1:2148689/52‑1 (MQ=255)
                                                          cGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGttt              >  1:3836725/1‑71 (MQ=255)
                                                          cGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGttt              <  1:3311719/71‑1 (MQ=255)
                                                           gCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGtttt             <  1:3181957/71‑1 (MQ=255)
                                                            ccTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAg                         <  1:105375/58‑1 (MQ=255)
                                                                 ccGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGtt       <  1:1105546/71‑1 (MQ=255)
                                                                    tGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGc                            <  1:1251879/47‑1 (MQ=255)
                                                                    tGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGt    <  1:561368/71‑1 (MQ=255)
                                                                     gCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGtt   >  1:3727133/1‑71 (MQ=255)
                                                                      cTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTcgcg                                 >  1:3797188/1‑40 (MQ=255)
                                                                      cTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTcgcg                                 >  1:1970032/1‑40 (MQ=255)
                                                                      cTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTcgcg                                 >  1:723388/1‑40 (MQ=255)
                                                                      cTCATTCCCTTCCGTTGGCAGGTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGttt  <  1:2197495/71‑1 (MQ=255)
                                                                      |                                                                      
CCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTT  >  minE/1015213‑1015353

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: