Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,780,768 A→G 22.3% A581A (GCA→GCG metH → homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,780,7680AG22.3% 86.6 / 18.7 67A581A (GCA→GCGmetHhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent
Reads supporting (aligned to +/- strand):  ref base A (17/35);  new base G (4/11);  total (21/46)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATGGATATGGGGATC  >  minE/2780709‑2780813
                                                           |                                             
gtttCTAACGTGTCGTTCTCGTTTCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCGGTGTTCctc                                      <  1:202547/66‑1 (MQ=255)
gtttCTAACGTGTCGTTCTCGTTTCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCGGTGTTCctc                                      <  1:119092/66‑1 (MQ=255)
gtttCTAACGTGTCGTTCTCGTTTCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCGGTGTTCctc                                      <  1:908776/66‑1 (MQ=255)
gtttCTAACGTGTCGTTCTCGTTTCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCGGTGTTCctc                                      <  1:2324775/66‑1 (MQ=255)
gtttCTAACGTGTCGTTCTCGTTTCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCGGTGTTCctc                                      <  1:3311004/66‑1 (MQ=255)
gtttCTAACGTGTCGTTCTCGTTTCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCGGTGTTCctc                                      <  1:865621/66‑1 (MQ=255)
gtttCTAACGTGTCGTTCTCGTTTCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCGGTGTTCctc                                      <  1:2212005/66‑1 (MQ=255)
gtttCTAACGTGTCGTTCTCGTTTCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCGGTGTTCctc                                      <  1:1685656/66‑1 (MQ=255)
gtttCTAACGTGTCGTTCTCGTTTCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCGGTGTTCctc                                      <  1:1706552/66‑1 (MQ=255)
gtttCTAACGTGTCGTTCTCGTTTCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCGGTGTTCctc                                      <  1:3129154/66‑1 (MQ=255)
gtttCTAACGTGTCGTTCTCGTTTCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCGGTGTTCctc                                      <  1:3417390/66‑1 (MQ=255)
gTATCTAACGTTTCTTTCTCGTTCCTTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc                                      <  1:3135273/69‑1 (MQ=255)
gTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc                                      <  1:3258733/69‑1 (MQ=255)
gTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc                                      <  1:1996381/69‑1 (MQ=255)
gTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc                                      <  1:68335/69‑1 (MQ=255)
gTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc                                      <  1:2771015/69‑1 (MQ=255)
gTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc                                      <  1:939471/69‑1 (MQ=255)
gTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCctc                                      <  1:992675/69‑1 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:2651895/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   <  1:110202/69‑1 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:2542307/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:3286608/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   <  1:3382963/69‑1 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:3408681/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:798480/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:868650/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:893037/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   <  1:93868/69‑1 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   <  1:240757/69‑1 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:1070948/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:1146821/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:1434707/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:1585296/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:1798058/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   <  1:2006356/69‑1 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:2210375/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   <  1:2459054/69‑1 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   <  1:2281706/69‑1 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:2256021/1‑69 (MQ=255)
   tCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGAGAAGCCATTCACGCAGTGTTCCTctac                                   >  1:2835969/1‑69 (MQ=255)
    cTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctat                                   <  1:2943791/68‑2 (MQ=255)
     tAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTctacta                                 >  1:2232641/1‑69 (MQ=255)
            tctttctCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTAtt                          >  1:3395649/1‑69 (MQ=255)
                  tcGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGc                       <  1:3221088/66‑1 (MQ=255)
                  tcGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGc                       <  1:1190182/66‑1 (MQ=255)
                  tcGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGc                       <  1:2676751/66‑1 (MQ=255)
                  tcGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGc                       <  1:500822/66‑1 (MQ=255)
                  tcGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGc                       <  1:361009/66‑1 (MQ=255)
                  tcGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGc                       <  1:3281100/66‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:1884117/68‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:2951187/68‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:986982/68‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:2585918/68‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:2805546/68‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:920740/68‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:1255847/68‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:1315900/68‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:3012249/68‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:146032/68‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:2753135/68‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:492799/68‑1 (MQ=255)
                         gTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATg              <  1:1939421/68‑1 (MQ=255)
                          tGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATGGa            <  1:3430446/69‑1 (MQ=255)
                                    ccGGTGCGCGAAGCCATTCACGCGGTGTTCCTCTACTACGCTATTCGCAATGGCATGGATATGGGGAt   >  1:1971567/1‑68 (MQ=255)
                                    ccGGTGCGCGAAGCCATTCACGCGGTGTTCCTCTACTACGCTATTCGCAATGGCATGGATATGGGGATc  >  1:832483/1‑69 (MQ=255)
                                    ccGGTGCGCGAAGCCATTCACGCGGTGTTCCTCTACTACGCTATTCGCAATGGCATGGATATGGGGATc  >  1:2037697/1‑69 (MQ=255)
                                    ccGGTGCGCGAAGCCATTCACGCGGTGTTCCTCTACTACGCTATTCGCAATGGCATGGATATGGGGATc  >  1:1836625/1‑69 (MQ=255)
                                                           |                                             
GTATCTAACGTTTCTTTCTCGTTCCGTGGCAACGATCCGGTGCGCGAAGCCATTCACGCAGTGTTCCTCTACTACGCTATTCGCAATGGCATGGATATGGGGATC  >  minE/2780709‑2780813

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: