Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,620 G→A 35.2% G483G (GGC→GGT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6200GA35.2% 53.0 / 44.6 54G483G (GGC→GGTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (13/22);  new base A (10/9);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGC  >  minE/1014552‑1014685
                                                                    |                                                                 
ggATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGcc                                                                 <  1:3217333/71‑1 (MQ=255)
        ctcacaTTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                          >  1:1360593/3‑70 (MQ=255)
        aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                          <  1:2641148/70‑1 (MQ=255)
        aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                          <  1:2520645/70‑1 (MQ=255)
        aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                          <  1:4134195/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)
                                                                    |                                                                 
GGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGC  >  minE/1014552‑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: