Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,545,006 A→G 100% V239A (GTC→GCC)  yfcY ← beta‑ketoacyl‑CoA thiolase, anaerobic, subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,545,0060AG77.2% 87.5 / 16.5 65V239A (GTC→GCC) yfcYbeta‑ketoacyl‑CoA thiolase, anaerobic, subunit
Reads supporting (aligned to +/- strand):  ref base A (14/1);  new base G (21/29);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 3.77e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAC  >  minE/1544968‑1545071
                                      |                                                                 
tagCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATAAAg                                      <  1:1586503/66‑1 (MQ=255)
cagCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATAAAg                                      <  1:3327845/66‑1 (MQ=255)
cagCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATAAAg                                      <  1:825358/66‑1 (MQ=255)
cagCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATAAAg                                      <  1:1239470/66‑1 (MQ=255)
cagCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATAAAg                                      <  1:1824042/66‑1 (MQ=255)
 agCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATAAAg                                      <  1:419728/66‑1 (MQ=255)
 agCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATAAAg                                      <  1:3230011/66‑1 (MQ=255)
 agCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATAAAg                                      <  1:1532572/66‑1 (MQ=255)
 agCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATAAAg                                      <  1:456914/66‑1 (MQ=255)
 agCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATAAAg                                      <  1:787273/66‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:2628046/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:1089250/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:2638511/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:2513837/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:2651812/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:2660997/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:2019242/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:2225949/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:2945077/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:2147860/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:3289810/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:3392774/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:1541741/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:1431450/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:1329351/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:579281/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:592911/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGGCAAGCGGTTGTTTATAag                                                 <  1:1193045/41‑1 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAATGCAGTCATCACCtcttc                     >  1:1609968/1‑69 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCtcttc                     >  1:2728335/1‑69 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCtcttc                     >  1:255651/1‑69 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCtcttc                     >  1:1123486/1‑69 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCtcttc                     >  1:1778167/1‑69 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCtcttc                     >  1:1219751/1‑69 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCtcttc                     >  1:1325864/1‑69 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCtcttc                     >  1:1243246/1‑69 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCtcttc                     >  1:1427316/1‑69 (MQ=255)
                aCCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGAAAAAGGCAGTCATCACct                         >  1:1283110/1‑65 (MQ=255)
                      aaTATTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATAAAg                                      <  1:2471115/46‑1 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:3401424/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:3045909/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:3068037/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:586443/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:1718502/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:582892/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:1583551/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:2816421/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:1323342/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:578396/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:2980623/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:1885758/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:1982646/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:2005045/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:2555241/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:2473736/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:945948/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:2101940/1‑43 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATaa                                        >  1:654927/1‑42 (MQ=255)
                        tatTGTTGTCTTCGGCAAGCGGTTGTTTATAAGGAGGGATa                                         >  1:575781/1‑41 (MQ=255)
                        tatTGTTGTCTTCGGAAAGCGGTTGTTTATAAGGAGGGATaaa                                       >  1:1877570/1‑43 (MQ=255)
                               gtCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCg        <  1:1835101/67‑1 (MQ=255)
                                  ttCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTg    >  1:2100979/1‑68 (MQ=255)
                                  ttCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTg    >  1:678349/1‑68 (MQ=255)
                                  ttCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTg    >  1:3277292/1‑68 (MQ=255)
                                  ttCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTg    >  1:223259/1‑68 (MQ=255)
                                      aCAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAc  <  1:561988/66‑1 (MQ=255)
                                      aCAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAc  <  1:1664592/66‑1 (MQ=255)
                                      aCAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAc  <  1:1733085/66‑1 (MQ=255)
                                      aCAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAc  <  1:2776013/66‑1 (MQ=255)
                                      aCAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAc  <  1:1943639/66‑1 (MQ=255)
                                      aCAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAc  <  1:201249/66‑1 (MQ=255)
                                      aCAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAc  <  1:854338/66‑1 (MQ=255)
                                      aCAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAc  <  1:870954/66‑1 (MQ=255)
                                      aCAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAc  <  1:2020807/66‑1 (MQ=255)
                                      aCAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAc  <  1:960748/66‑1 (MQ=255)
                                      aCAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAc  <  1:967227/66‑1 (MQ=255)
                                      aCAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAc  <  1:971239/66‑1 (MQ=255)
                                      |                                                                 
CGGCAAGCGAGGAATTACCGCGAATATTGTTGTCTTCGACAAGCGGTTGTTTATAAGGAGGGATAAAGGCAGTCATCACCTCTTCTTTGAGTTTTCCGTCTGAC  >  minE/1544968‑1545071

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