Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,336,510 A→G 49.4% Y137Y (TAT→TAC gmd ← GDP‑D‑mannose dehydratase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,336,5100AG49.4% ‑6.1 / 71.9 77Y137Y (TAT→TACgmdGDP‑D‑mannose dehydratase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (20/19);  new base G (13/25);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAGGAAGCGGATCGCC  >  minE/1336444‑1336576
                                                                  |                                                                  
gACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                  <  1:1062748/69‑1 (MQ=255)
gACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                  <  1:696065/69‑1 (MQ=255)
gACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                  <  1:3367114/69‑1 (MQ=255)
gACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                  <  1:3033870/69‑1 (MQ=255)
gACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                  <  1:2987160/69‑1 (MQ=255)
gACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                  <  1:2721312/69‑1 (MQ=255)
gACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                  <  1:2647213/69‑1 (MQ=255)
gACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                  <  1:2622958/69‑1 (MQ=255)
gACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATa                                                                  <  1:1369635/69‑1 (MQ=255)
         cGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCaga                                                         >  1:640465/1‑69 (MQ=255)
         cGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCaga                                                         >  1:1793300/1‑69 (MQ=255)
         cGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCaga                                                         >  1:1976396/1‑69 (MQ=255)
         cGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCag                                                          >  1:1594681/1‑68 (MQ=255)
                cgcgGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGaa                                                  >  1:1510018/1‑69 (MQ=255)
                cgcgGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGaa                                                  >  1:401164/1‑69 (MQ=255)
                cgcgGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGaa                                                  >  1:843180/1‑69 (MQ=255)
                   gggTAGAACGGCGTGGTCTC‑‑TCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:3087840/69‑1 (MQ=255)
                     gTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:98572/69‑1 (MQ=255)
                     gTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:883562/69‑1 (MQ=255)
                     gTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:908440/69‑1 (MQ=255)
                     gTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:561376/69‑1 (MQ=255)
                     gTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:2141339/69‑1 (MQ=255)
                     gTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:111895/69‑1 (MQ=255)
                     gTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:1598256/69‑1 (MQ=255)
                     gTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:154838/69‑1 (MQ=255)
                      tAGAACGGCGTTGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:2314412/68‑1 (MQ=255)
                      tAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:1372668/68‑1 (MQ=255)
                      tAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:709343/68‑1 (MQ=255)
                      tAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:1333200/68‑1 (MQ=255)
                      tAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:324092/68‑1 (MQ=255)
                      tAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:2535005/68‑1 (MQ=255)
                       aGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:1821269/67‑1 (MQ=255)
                                         tCTGCGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTg                                             <  1:2338226/49‑1 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCa                      >  1:3207854/1‑68 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAg                     >  1:2610666/1‑69 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAg                     >  1:3394482/1‑69 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAg                     >  1:1386301/1‑69 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAg                     >  1:3413038/1‑69 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAg                     >  1:343573/1‑69 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAg                     >  1:618168/1‑69 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAg                     >  1:818721/1‑69 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAg                     >  1:2098260/1‑69 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAg                     >  1:1268126/1‑69 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAg                     >  1:916488/1‑69 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAg                     >  1:541717/1‑69 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTATTTTCCAg                     >  1:2670077/1‑69 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGa                                   <  1:604137/55‑1 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGa                                   <  1:621964/55‑1 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGa                                   <  1:522220/55‑1 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGa                                   <  1:2718904/55‑1 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGa                                   <  1:1066147/55‑1 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGa                                   <  1:3275044/55‑1 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGa                                   <  1:2973951/55‑1 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGa                                   <  1:262926/55‑1 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGa                                   <  1:2617225/55‑1 (MQ=255)
                                             cGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGa                                   <  1:1804268/55‑1 (MQ=255)
                                                    tCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCcagacaga                >  1:143084/1‑67 (MQ=255)
                                                    tCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGa                >  1:3292932/1‑67 (MQ=255)
                                                    tCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGa                >  1:3264417/1‑67 (MQ=255)
                                                    tCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGa                >  1:3150901/1‑67 (MQ=255)
                                                    tCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGa                >  1:2296820/1‑67 (MQ=255)
                                                    tCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGa                >  1:2096032/1‑67 (MQ=255)
                                                    tCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGa                >  1:18335/1‑67 (MQ=255)
                                                    tCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGa                >  1:699259/1‑67 (MQ=255)
                                                            cAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAGACCGAGGAAGCGGAt      <  1:712933/69‑1 (MQ=255)
                                                            cAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAGACCGAGGAAGCGGAt      <  1:3027138/69‑1 (MQ=255)
                                                            cAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAGACCGAGGAAGCGGAt      <  1:709020/69‑1 (MQ=255)
                                                            cAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAGACCGAGGAAGCGGAt      <  1:1047191/69‑1 (MQ=255)
                                                            cAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAGACCGAGGAAGCGGAt      <  1:450914/69‑1 (MQ=255)
                                                            cAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAGACCGAGGAAGCGGAt      <  1:3231928/69‑1 (MQ=255)
                                                            cAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAGACCGAGGAAGCGGAt      <  1:1235791/69‑1 (MQ=255)
                                                            cAGACCGTACAGTTCAGAAGTAGAAGCCTGATAGAAACGGGTTTTCTTTTCCAGACCGAGGAAGCGGAt      <  1:2981961/69‑1 (MQ=255)
                                                                 cATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAGGAAGCGGATCGcc  >  1:734163/1‑68 (MQ=255)
                                                                 cATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAGGAAGCGGATCGcc  >  1:2621648/1‑68 (MQ=255)
                                                                 cATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAGGAAGCGGATCGcc  >  1:252337/1‑68 (MQ=255)
                                                                 cATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAGGAAGCGGATCGcc  >  1:2244345/1‑68 (MQ=255)
                                                                 cATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAGGAAGCGGATCGcc  >  1:1374276/1‑68 (MQ=255)
                                                                  |                                                                  
GACCGCATACGGAGATCGCGGGTAGAACGGCGTGGTCTCTTTCTGCGGAATTTCCTGCACCAGACCATACAGTTCAGAGGTGGAAGCCTGATAGAAACGAGTTTTCTTTTCCAGACCGAGGAAGCGGATCGCC  >  minE/1336444‑1336576

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: