Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,653 A→G 40.0% 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,6530AG40.0% 33.5 / 60.7 60D472D (GAT→GACfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (15/21);  new base G (11/13);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAT  >  minE/1014585‑1014723
                                                                    |                                                                      
cccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                                                                        <  1:5435483/69‑1 (MQ=255)
cccTGCGCCAATACAGCGGCCGGACCACCGATACAGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                                                                        <  1:3040774/69‑1 (MQ=255)
  cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcgg                                                                    >  1:5648506/1‑69 (MQ=255)
  cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcg                                                                     >  1:5141058/1‑69 (MQ=255)
  cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTc                                                                      <  1:227697/69‑1 (MQ=255)
  cTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                                                                        >  1:5498553/1‑67 (MQ=255)
  cTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCg                                                                    <  1:406601/71‑1 (MQ=255)
     cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGtaa                                                                 >  1:4676832/1‑69 (MQ=255)
     cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGta                                                                  <  1:2331422/70‑2 (MQ=255)
     cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGta                                                                  <  1:2271464/70‑2 (MQ=255)
     cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                                                                   <  1:3096114/69‑1 (MQ=255)
     cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                                                                   <  1:4717182/69‑1 (MQ=255)
     cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                                                                   <  1:2188830/69‑1 (MQ=255)
     cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTc                                                                  <  1:1776041/70‑1 (MQ=255)
        cAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                                              >  1:5625313/1‑71 (MQ=255)
        cAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                                              >  1:2937934/1‑71 (MQ=255)
         agtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACct                                                              >  1:3344664/3‑70 (MQ=255)
          gtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctgc                                                            >  1:401206/2‑71 (MQ=255)
          aTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCAt                                                                       <  1:5582778/60‑1 (MQ=255)
           tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATcc                                                                     <  1:2153796/61‑1 (MQ=255)
           tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgct                                                           >  1:477448/1‑71 (MQ=255)
               gccgCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTgctgact                                                       >  1:2804282/4‑71 (MQ=255)
                ccgCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTgctgactg                                                      >  1:1923016/3‑71 (MQ=255)
                   ccGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCgg                                                   >  1:1332908/1‑71 (MQ=255)
                     ggACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGtt                                                 <  1:3406118/71‑1 (MQ=255)
                        ccgccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTAcc                                              >  1:1428085/4‑71 (MQ=255)
                          accGATACTGCCAAGGTAGAAGCCTCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGtt                                                 <  1:1485907/66‑1 (MQ=255)
                           ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCttt                                           <  1:2124863/71‑1 (MQ=255)
                           ccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCa                                                                 >  1:5358990/1‑49 (MQ=255)
                              aTACTGCCGAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGc                                        >  1:4718648/1‑71 (MQ=255)
                               tACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCg                                                    <  1:584932/58‑1 (MQ=255)
                                   gCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATcc                                                                     >  1:1468642/1‑37 (MQ=255)
                                   gCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCAc                                                                >  1:1089860/1‑42 (MQ=255)
                                    ccgaGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCAt                                  >  1:4711726/4‑71 (MQ=255)
                                     cAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGt                                                  >  1:3420216/1‑54 (MQ=255)
                                     cAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATg                                 <  1:3314451/71‑1 (MQ=255)
                                      aaGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGa                                <  1:4642660/71‑1 (MQ=255)
                                         gTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGc                                        >  1:9102/1‑60 (MQ=255)
                                         gTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCa                             <  1:2046162/71‑1 (MQ=255)
                                          tAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGt                                                  <  1:3664711/49‑1 (MQ=255)
                                          tAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGt                                                  <  1:836799/49‑1 (MQ=255)
                                          tAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCa                             >  1:2626776/1‑70 (MQ=255)
                                          tAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCAt                            <  1:887462/71‑1 (MQ=255)
                                          tAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAt                            <  1:2955299/71‑1 (MQ=255)
                                           aGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCAt                            >  1:2383907/1‑70 (MQ=255)
                                           aGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtg                           >  1:875607/1‑70 (MQ=255)
                                            gAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtg                           <  1:2329419/70‑2 (MQ=255)
                                             aaGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtgc                          <  1:3357362/70‑3 (MQ=255)
                                             aaGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCAtg                           <  1:3328683/69‑2 (MQ=255)
                                             aaGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCt                         <  1:530231/71‑1 (MQ=255)
                                                 cgccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                     <  1:5149872/71‑1 (MQ=255)
                                                 cgccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                     <  1:4932291/71‑1 (MQ=255)
                                                 cgccgTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCg                     >  1:4858543/1‑71 (MQ=255)
                                                      tgtTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCAt                            <  1:4868250/59‑1 (MQ=255)
                                                      tgtTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCAt                            <  1:4861554/59‑1 (MQ=255)
                                                        ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCAt                            <  1:4304837/57‑1 (MQ=255)
                                                         tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgcgc             >  1:3348219/1‑71 (MQ=255)
                                                          ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTc                                          <  1:328727/41‑1 (MQ=255)
                                                             tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg          <  1:2412482/70‑1 (MQ=255)
                                                             tACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg          >  1:118167/1‑70 (MQ=255)
                                                               cAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc         <  1:3500963/69‑1 (MQ=255)
                                                               cAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCa        <  1:1257073/70‑1 (MQ=255)
                                                                   cATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa   >  1:1866316/1‑71 (MQ=255)
                                                                    aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt  >  1:785037/1‑71 (MQ=255)
                                                                    |                                                                      
CCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAT  >  minE/1014585‑1014723

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