Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,668 G→A 39.1% S467S (AGC→AGT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6680GA39.1% 32.5 / 45.2 46S467S (AGC→AGTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (12/16);  new base A (8/10);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCC  >  minE/1014601‑1014737
                                                                   |                                                                     
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:1255534/71‑1 (MQ=255)
       accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc                                                             <  1:2117837/71‑1 (MQ=255)
       accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc                                                             >  1:2014914/1‑71 (MQ=255)
          gccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCtt                                                          <  1:1882317/70‑1 (MQ=255)
           ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCttt                                                         >  1:3632253/1‑71 (MQ=255)
            cGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTc                                                        <  1:3852099/71‑1 (MQ=255)
            cGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTc                                                        <  1:3375926/71‑1 (MQ=255)
               tACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGcc                                                     <  1:1266158/71‑1 (MQ=255)
                  tGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCtt                                                          <  1:65776/63‑1 (MQ=255)
                    ccgaGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCAt                                                <  1:1872678/68‑1 (MQ=255)
                    ccgaGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCAt                                                <  1:2384966/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:2036719/71‑1 (MQ=255)
                                     gtgtTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                               >  1:3815333/1‑71 (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)
                                                         taaCCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCAcg           <  1:3545738/69‑2 (MQ=255)
                                                           aCCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTa         <  1:928210/71‑1 (MQ=255)
                                                             ctgctgGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCtt                         >  1:3947847/1‑53 (MQ=255)
                                                               gctggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGaa     >  1:345416/1‑71 (MQ=255)
                                                                ctggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGa      <  1:3414534/69‑1 (MQ=255)
                                                                ctggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGa      <  1:3387460/69‑1 (MQ=255)
                                                                ctggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGa      <  1:3302231/69‑1 (MQ=255)
                                                                   gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGa             <  1:2970784/59‑1 (MQ=255)
                                                                   gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatcc  <  1:664677/70‑1 (MQ=255)
                                                                   |                                                                     
CGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCC  >  minE/1014601‑1014737

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: