Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,653 A→G 12.1% 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,6530AG12.1% 81.4 / 5.0 33D472D (GAT→GACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (13/16);  new base G (2/2);  total (15/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGA  >  minE/1014586‑1014722
                                                                   |                                                                     
ccTGCGCCAATACAGCGGCCGGCCCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATcc                                                                    >  1:343947/1‑71 (MQ=255)
         gtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctgc                                                           >  1:2303198/2‑71 (MQ=255)
          tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTc                                                                 <  1:1881849/64‑1 (MQ=255)
           aCAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgctg                                                         >  1:4007781/1‑71 (MQ=255)
           aCAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgctg                                                         >  1:2401664/1‑71 (MQ=255)
              gCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctggct                                                      >  1:1782672/1‑71 (MQ=255)
               cGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTTCAGGCATCCGTCACCTgctggct                                                      >  1:1105988/1‑70 (MQ=255)
                ggCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGc                                                    >  1:3350080/1‑71 (MQ=255)
                  ccGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATc                                                                     >  1:3629039/1‑52 (MQ=255)
                          ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTc                                                                     >  1:446543/1‑44 (MQ=12)
                           cGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGt                                                                  <  1:708229/46‑1 (MQ=255)
                             aTACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGc                                       <  1:1980606/71‑1 (MQ=255)
                                cTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgc                                                           >  1:2080275/1‑48 (MQ=255)
                                   ccgaGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCt                                            <  1:2282755/57‑1 (MQ=255)
                                     aaGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGa                               >  1:586954/1‑71 (MQ=255)
                                      aGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGTCTGCGGTTGCCTTTCGCCAGCATGAt                              <  1:2841483/71‑1 (MQ=255)
                                      aGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGc                                   <  1:2502461/66‑1 (MQ=255)
                                      aGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGc                                   <  1:3700317/66‑1 (MQ=255)
                                       ggTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgct                                                          <  1:2800780/42‑1 (MQ=255)
                                           gAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATa                          <  1:562015/70‑1 (MQ=255)
                                                 gccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGc                   >  1:3228494/1‑71 (MQ=255)
                                                    gtgtTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGAt                              <  1:227388/57‑1 (MQ=255)
                                                       ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGa                               <  1:1671168/53‑1 (MQ=255)
                                                       ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgcg             >  1:4009513/1‑71 (MQ=255)
                                                       ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgc              <  1:1063912/70‑1 (MQ=255)
                                                         ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATc                             >  1:2770714/1‑53 (MQ=255)
                                                         ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgcgc            <  1:4126777/70‑1 (MQ=255)
                                                            tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg         <  1:1973623/70‑1 (MQ=255)
                                                            tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg         <  1:583684/70‑1 (MQ=255)
                                                            tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg         <  1:2135436/70‑1 (MQ=255)
                                                             aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc        <  1:426451/70‑1 (MQ=255)
                                                                 gCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGa                               <  1:2636776/43‑1 (MQ=255)
                                                                  cATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa  >  1:1971411/1‑71 (MQ=255)
                                                                   |                                                                     
CCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGA  >  minE/1014586‑1014722

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: