Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,298 T→C 11.5% P257P (CCA→CCG fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,2980TC11.5% 157.1 / 11.0 61P257P (CCA→CCGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (22/32);  new base C (3/4);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

CGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCAC  >  minE/1015229‑1015368
                                                                     |                                                                      
cgcTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGtt                                                                        <  1:3467956/70‑1 (MQ=255)
 gcTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTgg                                                                      >  1:1797452/1‑71 (MQ=255)
  cTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGc                                                                     <  1:279154/71‑1 (MQ=255)
   ttCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTgg                                                                      >  1:262112/1‑69 (MQ=255)
     cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                                                  >  1:638593/1‑71 (MQ=255)
     cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                                                  >  1:1635787/1‑71 (MQ=255)
     cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                                                  >  1:1567985/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)
                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)
                     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)
                                       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:1970032/1‑40 (MQ=255)
                                                      cTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTcgcg                                                >  1:3797188/1‑40 (MQ=255)
                                                      cTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTcgcg                                                >  1:723388/1‑40 (MQ=255)
                                                      cTCATTCCCTTCCGTTGGCAGGTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGttt                 <  1:2197495/71‑1 (MQ=255)
                                                       tcgttCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTc                >  1:1332389/4‑71 (MQ=255)
                                                       tcgttCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTc                >  1:3909374/4‑71 (MQ=255)
                                                       tCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGttt                                     <  1:3012355/50‑1 (MQ=255)
                                                         aTTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTc                <  1:167899/69‑1 (MQ=255)
                                                          ttCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTg              >  1:2860197/1‑70 (MQ=255)
                                                            cccTTCCGTTGGCAGTTCATCATAGTATTTCGCGGaa                                             >  1:1215210/1‑37 (MQ=255)
                                                             ccTTCCGTTGGCAGGTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGttt                 >  1:2532486/1‑64 (MQ=255)
                                                               ttCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAg           <  1:2010506/68‑1 (MQ=255)
                                                               ttCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAg        >  1:3420221/1‑71 (MQ=255)
                                                                 ccGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGt       <  1:2809015/70‑1 (MQ=255)
                                                                  cGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGtt      <  1:361565/70‑1 (MQ=255)
                                                                  cGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTc     <  1:3370886/71‑1 (MQ=255)
                                                                   gTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTcc    <  1:1300484/71‑1 (MQ=255)
                                                                     tGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGc             >  1:779978/1‑60 (MQ=255)
                                                                     tGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTccac  >  1:536662/1‑71 (MQ=255)
                                                                     |                                                                      
CGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTTCGTTTCTGCAGAAGTTCCAC  >  minE/1015229‑1015368

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: