Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,623 A→G 36.4% L482L (CTT→CTC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6230AG36.4% 44.7 / 49.5 55L482L (CTT→CTCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (14/21);  new base G (11/9);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGA  >  minE/1014560‑1014693
                                                               |                                                                      
ctcacaTTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                                  >  1:1360593/3‑70 (MQ=255)
aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                                  <  1:4134195/70‑1 (MQ=255)
aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                                  <  1:2641148/70‑1 (MQ=255)
aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                                  <  1:2520645/70‑1 (MQ=255)
aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                                  <  1:1980413/70‑1 (MQ=255)
       tCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCaa                                                                       >  1:5568517/1‑58 (MQ=255)
       tCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgccgcc                                                          >  1:924183/1‑71 (MQ=255)
          aGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtg                                                       <  1:3329862/71‑1 (MQ=255)
          aGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgt                                                        <  1:2919898/70‑1 (MQ=255)
            gCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgccgcc                                                          <  1:1879519/66‑1 (MQ=255)
               ctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGttttt                                                  >  1:4765620/1‑71 (MQ=255)
                ttAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCa                                                                        <  1:462459/48‑1 (MQ=255)
                  aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGa                                                                 >  1:1823041/1‑53 (MQ=255)
                    tACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAgg                                                                     <  1:3616730/47‑1 (MQ=255)
                       ttCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                                  <  1:162245/47‑1 (MQ=255)
                         cccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                                          <  1:5435483/69‑1 (MQ=255)
                         cccTGCGCCAATACAGCGGCCGGACCACCGATACAGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                                          <  1:3040774/69‑1 (MQ=255)
                           cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGcgt                                         <  1:2813926/68‑3 (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:2271464/70‑2 (MQ=255)
                              cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGta                                    <  1:2331422/70‑2 (MQ=255)
                              cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                                     <  1:4717182/69‑1 (MQ=255)
                              cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                                     <  1:3096114/69‑1 (MQ=255)
                              cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                                     <  1:2188830/69‑1 (MQ=255)
                              cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTc                                    <  1:1776041/70‑1 (MQ=255)
                                 cAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                >  1:2937934/1‑71 (MQ=255)
                                 cAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct                                >  1:5625313/1‑71 (MQ=255)
                                  agtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACct                                >  1:3344664/3‑70 (MQ=255)
                                   gtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctgc                              >  1:401206/2‑71 (MQ=255)
                                   aTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGccgccg                                                         >  1:3263390/1‑44 (MQ=255)
                                   aTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCAt                                         <  1:5582778/60‑1 (MQ=255)
                                    tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATcc                                       <  1:2153796/61‑1 (MQ=255)
                                    tACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgct                             >  1:477448/1‑71 (MQ=255)
                                     aCAGCGGCCGGACCACCGATACTGACAAGGTAGAAGCCGCCGTGTTTTTTAc                                               >  1:2034702/1‑52 (MQ=255)
                                        gccgCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTgctgact                         >  1:2804282/4‑71 (MQ=255)
                                         ccgCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTgctgactg                        >  1:1923016/3‑71 (MQ=255)
                                          ggCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTAc                                               <  1:3324028/47‑1 (MQ=255)
                                            ccGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCgg                     >  1:1332908/1‑71 (MQ=255)
                                              ggACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGtt                   <  1:3406118/71‑1 (MQ=255)
                                              ggACCACCGATACTGCCAAGGTAGAAGCCGCCGTGttttt                                                  <  1:2221431/40‑1 (MQ=255)
                                                 ccgccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTAcc                >  1:1428085/4‑71 (MQ=255)
                                                  caccGATACTGCCAAGGTAGAAGCCGCCGTGttttt                                                  <  1:5541609/36‑1 (MQ=255)
                                                   accGATACTGCCAAGGTAGAAGCCTCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGtt                   <  1:1485907/66‑1 (MQ=255)
                                                    ccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCttt             <  1:2124863/71‑1 (MQ=255)
                                                    ccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCa                                   >  1:5358990/1‑49 (MQ=255)
                                                     cGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACa                                              >  1:1968719/1‑37 (MQ=16)
                                                       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)
                                                              cAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGt                    >  1:3420216/1‑54 (MQ=255)
                                                              cAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATg   <  1:3314451/71‑1 (MQ=255)
                                                               aaGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGa  <  1:4642660/71‑1 (MQ=255)
                                                               |                                                                      
AACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGA  >  minE/1014560‑1014693

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: