Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,280 G→A 11.9% 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,2800GA11.9% 157.6 / 10.5 59H263H (CAC→CATfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base A (2/5);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

CACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTT  >  minE/1015210‑1015348
                                                                      |                                                                    
cACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCg                                                                      <  1:2768556/71‑1 (MQ=255)
cACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCg                                                                      <  1:3950693/71‑1 (MQ=255)
   ccAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGc                                                                   >  1:1908612/1‑71 (MQ=255)
    cAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCg                                                                      <  1:2190500/67‑1 (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)
                        cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTg                                                                    >  1:2446395/1‑49 (MQ=255)
                        cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                              >  1:1635787/1‑71 (MQ=255)
                        cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                              >  1:1567985/1‑71 (MQ=255)
                        cGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                              >  1:638593/1‑71 (MQ=255)
                         gATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                              <  1:1938763/70‑1 (MQ=255)
                         gATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGt                                              <  1:1912878/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:3311719/71‑1 (MQ=255)
                                                             cGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGttt         >  1:3836725/1‑71 (MQ=255)
                                                              gCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGtttt        <  1:3181957/71‑1 (MQ=255)
                                                               ccTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAg                    <  1:105375/58‑1 (MQ=255)
                                                                    ccGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGtt  <  1:1105546/71‑1 (MQ=255)
                                                                      |                                                                    
CACCCAGACCAAGATTTTGCGCTTCGATCAGCAATTCTTTTTCCAGTTCCACATCGCGGAACGCCTGACCGTGCTCATTCCCTTCCGTTGGCAGTTCATCATAGTATTTCGCGGAAGCCAGTTTCACCGTTTTAAGGTT  >  minE/1015210‑1015348

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: