Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,653 A→G 41.2% D472D (GAT→GAC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6530AG41.2% 23.5 / 56.8 51D472D (GAT→GACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (12/18);  new base G (13/8);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGA  >  minE/1014585‑1014722
                                                                    |                                                                     
tgcTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTc                                                                     >  1:689321/3‑71 (MQ=255)
cccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                                                                       <  1:1169085/69‑1 (MQ=255)
  cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcgg                                                                   >  1:3057193/1‑69 (MQ=255)
        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:2117837/71‑1 (MQ=255)
                       accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc                                              <  1:1255534/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:3852099/71‑1 (MQ=255)
                            cGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTc                                         <  1:3375926/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)
                                                                    |                                                                     
CCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGA  >  minE/1014585‑1014722

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: