Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,622,786 G→A 23.2% V620V (GTC→GTT maeB ← fused malic enzyme predicted oxidoreductase and predicted phosphotransacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,622,7860GA23.2% 120.4 / 27.7 90V620V (GTC→GTTmaeBfused malic enzyme predicted oxidoreductase and predicted phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base G (26/43);  new base A (14/7);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 2.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

CAGTCAGAAGAACCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGT  >  minE/1622721‑1622854
                                                                 |                                                                    
cAGTCAGAAGAACCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACa                                                                   <  1:1047569/69‑1 (MQ=255)
cAGTCAGAAGAACCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACa                                                                   <  1:1586941/69‑1 (MQ=255)
cAGTCAGAAGAACCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACa                                                                   <  1:677762/69‑1 (MQ=255)
cAGTCAGAAGAACCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACa                                                                   <  1:2575243/69‑1 (MQ=255)
cAGTCAGAAGAACCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACa                                                                   <  1:1340626/69‑1 (MQ=255)
cAGTCAGAAGAACCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACa                                                                   <  1:2055830/69‑1 (MQ=255)
cAGTCAGAAGAACCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACa                                                                   <  1:283970/69‑1 (MQ=255)
cAGTCAGAAGAACCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACa                                                                   <  1:3307588/69‑1 (MQ=255)
 aGTCAGAAGAACCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACa                                                                   <  1:1028001/68‑1 (MQ=255)
       aagaaCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAATCGACGGACAGTTTCTg                                                            >  1:1416831/1‑69 (MQ=255)
       aagaaCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCt                                                             >  1:3346797/1‑68 (MQ=255)
       aagaaCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCt                                                             >  1:1237461/1‑68 (MQ=255)
       aagaaCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCt                                                             >  1:913946/1‑68 (MQ=255)
       aagaaCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTg                                                            >  1:496521/1‑69 (MQ=255)
       aagaaCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTg                                                            >  1:3264058/1‑69 (MQ=255)
       aagaaCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTg                                                            >  1:2062090/1‑69 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTc                                                              >  1:2371016/1‑66 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:1912248/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:1710183/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:2108726/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:1755462/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:1966508/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:1476204/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:1369787/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:2794306/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:1338927/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:1234684/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:1210821/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:582224/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAGAGCCACACGCGGCTCAATACCAAAACGACGAACAGTTTCt                                                             >  1:712193/1‑67 (MQ=255)
        agaaCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTg                                                            >  1:250222/1‑68 (MQ=255)
           aCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTgccgc                                                        <  1:2688252/69‑1 (MQ=255)
           aCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTgccgc                                                        <  1:2446741/69‑1 (MQ=255)
           aCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTgccgc                                                        <  1:2442884/69‑1 (MQ=255)
           aCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTgccgc                                                        <  1:2806175/69‑1 (MQ=255)
           aCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTgccgc                                                        <  1:133346/69‑1 (MQ=255)
           aCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTgccgc                                                        <  1:2552385/69‑1 (MQ=255)
           aCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTgccgc                                                        <  1:1965640/69‑1 (MQ=255)
           aCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTgccgc                                                        <  1:3057996/69‑1 (MQ=255)
           aCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTgccgc                                                        <  1:2114340/69‑1 (MQ=255)
           aCCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTgccgc                                                        <  1:72926/69‑1 (MQ=255)
             cAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGAGGGACa                                                                   <  1:2465229/56‑1 (MQ=255)
                    ggAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAATGACGGa                                                                     <  1:436739/47‑1 (MQ=255)
                    ggAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGa                                                                     <  1:1588996/47‑1 (MQ=255)
                    ggAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGa                                                                     <  1:529473/47‑1 (MQ=255)
                    ggAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGa                                                                     <  1:876986/47‑1 (MQ=255)
                    ggAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGa                                                                     <  1:2500164/47‑1 (MQ=255)
                    ggAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGa                                                                     <  1:1344125/47‑1 (MQ=255)
                    ggAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGa                                                                     <  1:2756628/47‑1 (MQ=255)
                     gAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGa                                                                     <  1:804240/46‑1 (MQ=255)
                                  aaGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGct                                 >  1:2906951/1‑69 (MQ=255)
                                  aaGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGc                                  >  1:981131/1‑68 (MQ=255)
                                  aaGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTGCCGCCATAAAGGTGATCTCCGCCAGct                                 >  1:2951846/1‑69 (MQ=255)
                                   aGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTgccgc                                                        <  1:2866256/45‑1 (MQ=255)
                                        cacgcgGCTCAATACCAAAACGACGAACAGTTTCTGCCGCCATCAAGGTGAtct                                          <  1:3358009/52‑1 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc              >  1:2693684/1‑68 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc              >  1:1349014/1‑68 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc              >  1:1066055/1‑68 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc              >  1:739402/1‑68 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc              >  1:1140680/1‑68 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc              >  1:2816531/1‑68 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc              >  1:2839139/1‑68 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc              >  1:1481388/1‑68 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc              >  1:1149543/1‑68 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc              >  1:3389829/1‑68 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcatc              >  1:2225783/1‑68 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcat               >  1:1814415/1‑67 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcat               >  1:1621865/1‑67 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcat               >  1:2503416/1‑67 (MQ=255)
                                                      ccAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTtcat               >  1:2349717/1‑67 (MQ=255)
                                                        aaaaCGACGAACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCAGGCTCATCAtt           <  1:899314/69‑1 (MQ=255)
                                                        aaaaCGACGAACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCAGGCTCATCAtt           <  1:882526/69‑1 (MQ=255)
                                                        aaaaCGACGAACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCAGGCTCATCAtt           <  1:2011391/69‑1 (MQ=255)
                                                        aaaaCGACGAACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCAGGCTCATCAtt           <  1:1318082/69‑1 (MQ=255)
                                                        aaaaCGACGAACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCAGGCTCATCAtt           <  1:2039153/69‑1 (MQ=255)
                                                        aaaaCGACGAACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCAGGCTCATCAtt           <  1:2396684/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:2706509/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:980800/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:871230/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:865923/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:1019936/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:1624093/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:287371/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:1283077/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:2173916/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:3025529/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:3317568/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:1327998/69‑1 (MQ=255)
                                                            cgacgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACa       <  1:3302927/69‑1 (MQ=255)
                                                              acgGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACata     <  1:602429/69‑1 (MQ=255)
                                                                 gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt  <  1:439527/69‑1 (MQ=255)
                                                                 gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt  <  1:1124038/69‑1 (MQ=255)
                                                                 gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt  <  1:391823/69‑1 (MQ=255)
                                                                 gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt  <  1:391619/69‑1 (MQ=255)
                                                                 gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt  <  1:2246157/69‑1 (MQ=255)
                                                                 gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt  <  1:3385806/69‑1 (MQ=255)
                                                                 gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt  <  1:3337888/69‑1 (MQ=255)
                                                                 gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt  <  1:1333329/69‑1 (MQ=255)
                                                                 gACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGt  <  1:3181425/69‑1 (MQ=255)
                                                                 |                                                                    
CAGTCAGAAGAACCAAAGTTGGAGTGCGACAACAAAGCAACGCGCGGCTCAATACCAAAACGACGGACAGTTTCTGCCGCCATCAAGGTGATCTCCGCCAGCTCTTCTGCATCCGGTTCATCATTAACATATGT  >  minE/1622721‑1622854

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: