Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,617 A→G 34.5% S484S (AGT→AGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6170AG34.5% 52.0 / 46.7 55S484S (AGT→AGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (13/23);  new base G (10/9);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGC  >  minE/1014548‑1014685
                                                                     |                                                                    
ttCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATAc                                                                     <  1:607906/71‑1 (MQ=255)
    ggATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGcc                                                                 <  1:3217333/71‑1 (MQ=255)
            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)
                             aaaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATa                                                                      >  1:4431457/2‑41 (MQ=255)
                              aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATa                                                                      >  1:4257964/1‑40 (MQ=255)
                              aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATa                                                                      >  1:1266360/1‑40 (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:3096114/69‑1 (MQ=255)
                                          cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                             <  1:2188830/69‑1 (MQ=255)
                                          cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt                             <  1:4717182/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)
                                               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)
                                                                     |                                                                    
TTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGC  >  minE/1014548‑1014685

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: