Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,509,579 G→A 30.9% A42T (GCA→ACA)  ackA → acetate kinase A and propionate kinase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,509,5790GA30.9% 79.8 / 42.6 55A42T (GCA→ACA) ackAacetate kinase A and propionate kinase 2
Reads supporting (aligned to +/- strand):  ref base G (16/22);  new base A (5/12);  total (21/34)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCATCGATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAGGTGCAGGCGCCGCTCA  >  minE/1509513‑1509643
                                                                  |                                                                
atcatcGATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACg                                                              <  1:74760/71‑1 (MQ=255)
 tcatcGATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGt                                                             <  1:1997282/71‑1 (MQ=255)
     cGATGCAGTAAATGGTGAAGATTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATc                                                          >  1:185271/1‑70 (MQ=255)
     cGATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATc                                                          >  1:503030/1‑70 (MQ=255)
      gATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaa                                                        <  1:1972831/71‑1 (MQ=255)
      gATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaa                                                        <  1:1282013/71‑1 (MQ=255)
      gATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaa                                                        <  1:1284660/71‑1 (MQ=255)
      gATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCa                                                         >  1:737569/1‑70 (MQ=255)
      gATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCa                                                         <  1:2061452/70‑1 (MQ=255)
      gATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCa                                                         <  1:851385/70‑1 (MQ=255)
            gTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCaaatgga                                                   >  1:1342000/1‑70 (MQ=255)
            gTAAATGGTGAAGAGTACCTTTCTGGTTTAGACGAATGTTTCCACCTGCCCGAAGCACGTATCa                                                         <  1:1343242/64‑1 (MQ=255)
               aaTGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaaaa                                                <  1:1982101/70‑1 (MQ=255)
               aaTGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCaaatggaaaat                                               >  1:1354607/1‑71 (MQ=255)
               aaTGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCaaatggaaaat                                               >  1:156113/1‑71 (MQ=255)
                aTGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaaaatg                                              <  1:1095843/71‑1 (MQ=255)
                aTGGGGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaaaatg                                              <  1:467874/71‑1 (MQ=255)
                  ggTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaaaatgga                                            >  1:42599/1‑71 (MQ=255)
                  ggTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaaaatgga                                            >  1:1893927/1‑71 (MQ=255)
                   gTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaaaatgg                                             >  1:303430/1‑69 (MQ=255)
                      aaGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatgg                                                    <  1:393136/59‑1 (MQ=255)
                        gagTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCaaatggaa                                                  >  1:474362/1‑59 (MQ=255)
                        gagTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCaa                                      <  1:2200792/71‑1 (MQ=255)
                         agTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCaa                                      >  1:1305914/1‑70 (MQ=255)
                          gTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATa                                    <  1:388835/71‑1 (MQ=255)
                           tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCa                                                         <  1:869092/49‑1 (MQ=255)
                           tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCa                                                         <  1:10670/49‑1 (MQ=255)
                           tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCa                                                         <  1:1111094/49‑1 (MQ=255)
                           tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATa                                    >  1:1060067/1‑70 (MQ=255)
                           tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATa                                    >  1:1431783/1‑70 (MQ=255)
                           tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCa                                                         <  1:1693450/49‑1 (MQ=255)
                           tACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCa                                                         <  1:232851/49‑1 (MQ=255)
                             ccTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATaaa                                  <  1:50033/70‑1 (MQ=255)
                             ccTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATaaa                                  <  1:306416/70‑1 (MQ=255)
                                 tCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCaaatggaaaa                                                >  1:41597/1‑52 (MQ=255)
                                  cTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGa                             <  1:686138/70‑1 (MQ=255)
                                  cTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGa                             <  1:1073730/70‑1 (MQ=255)
                                  cTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGa                             <  1:234917/70‑1 (MQ=255)
                                  cTGGTTTAGCCGAATGTTTCCACCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATAAACAGGa                             <  1:1227913/70‑1 (MQ=255)
                                     gTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATaaa                                  >  1:232628/1‑62 (MQ=255)
                                          gCCGAATGTTTCCACCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGAGGCttt                    <  1:370385/71‑1 (MQ=255)
                                             gAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAgg                 >  1:28689/1‑71 (MQ=255)
                                             gAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAgg                 >  1:26102/1‑71 (MQ=255)
                                               aTGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAgg                 >  1:1359644/1‑69 (MQ=255)
                                               aTGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAg                  >  1:1637691/1‑68 (MQ=255)
                                                  tttCCACCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCttt                    <  1:1972934/63‑1 (MQ=255)
                                                       aCCTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATa                                    <  1:1933316/42‑1 (MQ=255)
                                                        ccTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCttt                    <  1:2148543/57‑1 (MQ=255)
                                                        ccTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCttt                    <  1:498338/57‑1 (MQ=255)
                                                         cTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATaaa                                  >  1:1586534/1‑42 (MQ=255)
                                                         cTGCCCGAAACACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAGGt                <  1:1178269/60‑1 (MQ=255)
                                                            cccGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAGGTGCAGGCGCCGCt    <  1:1993433/69‑1 (MQ=255)
                                                            cccGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAGGTGCAGGCGCCGCTca  <  1:63066/71‑1 (MQ=255)
                                                             ccGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCttt                    <  1:1806343/52‑1 (MQ=255)
                                                                aaGCACGTATCAAATGGAAAATGGACGGCAATaaa                                  >  1:291289/1‑35 (MQ=255)
                                                                  gCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTa                   <  1:1231779/48‑1 (MQ=255)
                                                                  |                                                                
ATCATCGATGCAGTAAATGGTGAAGAGTACCTTTCTGGTTTAGCCGAATGTTTCCACCTGCCCGAAGCACGTATCAAATGGAAAATGGACGGCAATAAACAGGAAGCGGCTTTAGGTGCAGGCGCCGCTCA  >  minE/1509513‑1509643

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: