Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,656 C→G 37.2% T471T (ACG→ACC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6560CG37.2% 40.4 / 49.0 51T471T (ACG→ACCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (13/19);  new base G (11/8);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATC  >  minE/1014593‑1014724
                                                               |                                                                    
cagtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACct                                                               >  1:1661580/4‑71 (MQ=255)
  gtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctg                                                              >  1:283608/2‑70 (MQ=255)
   tACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctgct                                                            >  1:3128585/1‑71 (MQ=255)
     cAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgctg                                                           <  1:3721172/70‑1 (MQ=255)
        cGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctggctg                                                       <  1:2730593/71‑1 (MQ=255)
        cGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctggct                                                        >  1:3006410/1‑70 (MQ=255)
          gCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCg                                                     >  1:1321545/1‑71 (MQ=255)
           ccGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCgg                                                    <  1:421606/71‑1 (MQ=255)
            cGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGt                                                   <  1:1790629/71‑1 (MQ=255)
             ggACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGtt                                                  <  1:367644/71‑1 (MQ=255)
             ggACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGt                                                   >  1:2382111/1‑70 (MQ=255)
               accgccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGttac                                                >  1:2005342/1‑69 (MQ=255)
               accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc                                                >  1:2014914/1‑71 (MQ=255)
               accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc                                                <  1:1255534/71‑1 (MQ=255)
               accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc                                                <  1:2117837/71‑1 (MQ=255)
                ccaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGt                                                                    >  1:3467307/1‑50 (MQ=255)
                  gccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCtt                                             <  1:1882317/70‑1 (MQ=255)
                   ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCttt                                            >  1:3632253/1‑71 (MQ=255)
                    cGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTc                                           <  1:3375926/71‑1 (MQ=255)
                    cGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTc                                           <  1:3852099/71‑1 (MQ=255)
                       tACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGcc                                        <  1:1266158/71‑1 (MQ=255)
                        aCTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTc                                                                   <  1:1152728/43‑1 (MQ=255)
                          tGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCtt                                             <  1:65776/63‑1 (MQ=255)
                            ccgaGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCAt                                   <  1:2384966/68‑1 (MQ=255)
                            ccgaGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCAt                                   <  1:1872678/68‑1 (MQ=255)
                            ccAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCAt                                   >  1:1803096/1‑71 (MQ=255)
                             cAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATg                                  <  1:208981/71‑1 (MQ=255)
                               aGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGAt                                <  1:113469/71‑1 (MQ=255)
                                ggTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGAt                                >  1:2152047/1‑70 (MQ=255)
                                   aGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCa                              <  1:1381217/69‑1 (MQ=255)
                                   aGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtg                            >  1:2552918/1‑70 (MQ=255)
                                    gAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtgc                           >  1:582794/1‑69 (MQ=255)
                                      aGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACtt                         <  1:2708418/71‑1 (MQ=255)
                                       gccgccGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCtgc                        >  1:2210878/1‑69 (MQ=255)
                                         cgccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                      >  1:3283334/1‑71 (MQ=255)
                                         cgccgTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCg                      >  1:520368/1‑71 (MQ=255)
                                            cgTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCAt                                   <  1:971827/55‑1 (MQ=255)
                                             gtgtTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                  >  1:3815333/1‑71 (MQ=255)
                                             gtgtTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                  <  1:2036719/71‑1 (MQ=255)
                                                ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                      >  1:434409/1‑64 (MQ=255)
                                                  ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                  >  1:2609980/1‑66 (MQ=255)
                                                    ttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTc                        >  1:3651580/1‑58 (MQ=255)
                                                     tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg           <  1:407881/70‑1 (MQ=255)
                                                      acacGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGc          <  1:3720889/70‑1 (MQ=255)
                                                      aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc          <  1:1571079/70‑1 (MQ=255)
                                                      aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCCCTTGc          <  1:2075309/70‑1 (MQ=255)
                                                          gCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGtt      >  1:3113426/1‑70 (MQ=255)
                                                          gCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGt       >  1:1543762/1‑69 (MQ=255)
                                                           cATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa    <  1:3804464/71‑1 (MQ=255)
                                                            aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGtt      <  1:3298396/68‑1 (MQ=255)
                                                             tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc  >  1:2074443/1‑71 (MQ=255)
                                                               |                                                                    
CAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATC  >  minE/1014593‑1014724

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: