Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,511,011 A→G 25.8% A93A (GCA→GCG pta → phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,511,0110AG25.8% 122.7 / 41.5 108A93A (GCA→GCGptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (46/34);  new base G (16/12);  total (62/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 108 total aligned reads displayed.

TTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTCCCGACACGTAAG  >  minE/1510956‑1511077
                                                       |                                                                  
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:190050/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:1644580/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:3410361/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:1542038/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:3186633/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:1694638/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:3186175/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:3152212/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:1807001/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:1324457/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:1960553/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:809158/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:264249/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:123137/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:953628/68‑1 (MQ=255)
ttGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:1005656/68‑1 (MQ=255)
 tGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGCGATCGTCGCAAACTACCACGCTaa                                                      <  1:395753/69‑1 (MQ=255)
 tGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTaa                                                      <  1:1174661/69‑1 (MQ=255)
 tGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTaa                                                      <  1:2538194/69‑1 (MQ=255)
       tCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAacac                                                   >  1:165034/1‑66 (MQ=255)
       tCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAacac                                                   >  1:1798695/1‑66 (MQ=255)
       tCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAacac                                                   >  1:2029940/1‑66 (MQ=255)
       tCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAacac                                                   >  1:2064778/1‑66 (MQ=255)
       tCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAacac                                                   >  1:2095673/1‑66 (MQ=255)
       tCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAacac                                                   >  1:809800/1‑66 (MQ=255)
       tCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAacac                                                   >  1:527491/1‑66 (MQ=255)
       tCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAacac                                                   >  1:1229980/1‑66 (MQ=255)
       tCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAacac                                                   >  1:1095460/1‑66 (MQ=255)
       tCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAacac                                                   >  1:1252486/1‑66 (MQ=255)
        cTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCaaa                                               >  1:602089/1‑69 (MQ=255)
             ttCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGAc                                            <  1:2430032/67‑1 (MQ=255)
                           aaGATGTTCTGATGGAAGAGATCGTCGCAAACTACCACGCt                                                        <  1:768710/41‑1 (MQ=255)
                                        ggAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTgg                <  1:1351389/68‑1 (MQ=255)
                                        ggAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTgg                <  1:603363/68‑1 (MQ=255)
                                         gAAGAGATCGTCGCAAACTAGCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:1482475/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCCCGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:140436/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCg                                  >  1:3043505/1‑49 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGt               >  1:2571456/1‑68 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:2489172/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:2284279/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:964412/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:1080465/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:2206731/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:2633229/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:269440/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:2708927/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:3019804/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:582646/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:3376560/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:1797725/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:1873249/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:2092852/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:144344/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTc              >  1:2075573/1‑69 (MQ=255)
                                         gAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTAAGGTCTGGTc              >  1:1276053/1‑68 (MQ=255)
                                           agagATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGaa                                      >  1:1661393/1‑43 (MQ=255)
                                           agagATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTccc            >  1:827994/1‑69 (MQ=255)
                                           agagATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTccc            >  1:2412582/1‑69 (MQ=255)
                                           agagATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTcc             >  1:3306607/1‑68 (MQ=255)
                                             agaTCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTCCCGa          <  1:2445728/69‑1 (MQ=255)
                                             agaTCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTCCCGa          <  1:2449690/69‑1 (MQ=255)
                                              gaTCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTCCCGa          <  1:1371957/68‑1 (MQ=255)
                                                 cgtcgCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGt               >  1:3068598/1‑60 (MQ=255)
                                                 cgtcgCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGt               >  1:3316774/1‑60 (MQ=255)
                                                 cgtcgCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGt               >  1:3339707/1‑60 (MQ=255)
                                                 cgtcgCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGt               >  1:1608273/1‑60 (MQ=255)
                                                 cgtcgCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGt               >  1:1370869/1‑60 (MQ=255)
                                                 cgtcgCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTCCCGacac       >  1:1354242/1‑68 (MQ=255)
                                                 cgtcgCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTCCCGacac       >  1:1146426/1‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:576260/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:1188278/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:3236378/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:1724227/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:2273226/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:1787647/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:2114072/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:1740482/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:2447502/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:2014646/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:2505301/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:2064293/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:2118269/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:2117152/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGacac       >  1:2084731/3‑68 (MQ=255)
                                                 catcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGaca        >  1:2302412/3‑67 (MQ=255)
                                                  atcgCGAACTACCACGCTAACACCAAAGACGCTGAATTCGTTCTGGTGGAAGGTCTGGTCCCGACACGt     <  1:2520232/68‑1 (MQ=255)
                                                  atcgCGAACTACCACGCTAACACCAAAGACGCTGAATTCGTTCTGGTGGAAGGTCTGGTCCCGACACGt     <  1:3357000/68‑1 (MQ=255)
                                                  atcgCGAACTACCACGCTAACACCAAAGACGCTGAATTCGTTCTGGTGGAAGGTCTGGTCCCGACACGt     <  1:1684731/68‑1 (MQ=255)
                                                  atcgCGAACTACCACGCTAACACCAAAGACGCTGAATTCGTTCTGGTGGAAGGTCTGGTCCCGACACGt     <  1:2394390/68‑1 (MQ=255)
                                                  atcgCGAACTACCACGCTAACACCAAAGACGCTGAATTCGTTCTGGTGGAAGGTCTGGTCCCGACACGt     <  1:1279188/68‑1 (MQ=255)
                                                  atcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGACACGt     <  1:3323674/68‑1 (MQ=255)
                                                  atcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGACACGt     <  1:1893342/68‑1 (MQ=255)
                                                  atcgCGAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTGGAAGGTCTGGTCCCGACACGt     <  1:636408/68‑1 (MQ=255)
                                                  atcgCGAACTACCACGCTAAAACCAAAGACGCTGAATTCGTTCTGGTGGAAGGTCTGGTCCCGACACGt     <  1:2551028/68‑1 (MQ=255)
                                                   tcgCGAACTACCACGCTAACACCAAAGACGCTGAATTCGTTCTGGTGGAAGGTCTGGTCCCGACACGt     <  1:1111903/68‑1 (MQ=255)
                                                   tcgCGAACTACCACGCTAACACCAAAGACGCTGAATTCGTTCTGGTGGAAGGTCTGGTCCCGACACGt     <  1:125454/68‑1 (MQ=255)
                                                   tcgCGAACTACCACGCTAACACCAAAGACGCTGAATTCGTTCTGGTGGAAGGTCTGGTCCCGACACGt     <  1:3345336/68‑1 (MQ=255)
                                                     gCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTCCCGACTCGTAAg  <  1:437936/69‑1 (MQ=255)
                                                     gCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTCCCGACACGTAAg  <  1:2293796/69‑1 (MQ=255)
                                                     gCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTCCCGACACGTAAg  <  1:2720206/69‑1 (MQ=255)
                                                       |                                                                  
TTGAAGGTCTGCTTTCCAGCAATCAGAAAGATGTGCTGATGGAAGAGATCGTCGCAAACTACCACGCTAACACCAAAGACGCTGAAGTCGTTCTGGTTGAAGGTCTGGTCCCGACACGTAAG  >  minE/1510956‑1511077

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: