Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,650 G→C 42.9% A473A (GCC→GCG fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6500GC42.9% 17.8 / 57.1 49A473A (GCC→GCGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (10/18);  new base C (13/8);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 8.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGC  >  minE/1014585‑1014716
                                                                 |                                                                  
tgcTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTc                                                               >  1:689321/3‑71 (MQ=255)
cccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                                                                 <  1:1169085/69‑1 (MQ=255)
  cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGcg                                                                 <  1:3625366/67‑2 (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:2117837/71‑1 (MQ=255)
                       accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc                                        <  1:1255534/71‑1 (MQ=255)
                       accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc                                        >  1:2014914/1‑71 (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: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: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)
                                                                 |                                                                  
CCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGC  >  minE/1014585‑1014716

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: