Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 108,755 C→T 48.0% Y333Y (TAC→TAT aceE → pyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE108,7550CT48.0% 1.0 / 43.0 52Y333Y (TAC→TATaceEpyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base C (11/16);  new base T (9/16);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

CTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAC  >  minE/108687‑108818
                                                                    |                                                               
cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc                                                                 >  1:3218466/1‑69 (MQ=255)
cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc                                                                 >  1:2417503/1‑69 (MQ=255)
cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc                                                                 >  1:2359416/1‑69 (MQ=255)
cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc                                                                 >  1:1301476/1‑69 (MQ=255)
cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc                                                                 >  1:3122267/1‑69 (MQ=255)
cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc                                                                 >  1:1369641/1‑69 (MQ=255)
cTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTAc                                                                 >  1:1843019/1‑69 (MQ=255)
     ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt                                                              <  1:2832595/67‑1 (MQ=255)
     ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt                                                              <  1:1043745/67‑1 (MQ=255)
     ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt                                                              <  1:449882/67‑1 (MQ=255)
     ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt                                                              <  1:2511369/67‑1 (MQ=255)
     ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt                                                              <  1:3128195/67‑1 (MQ=255)
     ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt                                                              <  1:1348829/67‑1 (MQ=255)
     ccAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATAGAAAGATGGTGCGTACGtt                                                              <  1:2948808/67‑1 (MQ=255)
               aaCGAAACCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACActtc                                                  <  1:766574/69‑1 (MQ=255)
               aaCGAAACCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACActtc                                                  <  1:2926386/69‑1 (MQ=255)
                aCGAAACCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACActtc                                                  <  1:379544/68‑1 (MQ=255)
                aCGAAACCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACActtc                                                  <  1:490651/68‑1 (MQ=255)
                    aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                              >  1:2318796/1‑68 (MQ=255)
                    aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                              >  1:3367725/1‑68 (MQ=255)
                    aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                              >  1:397446/1‑68 (MQ=255)
                    aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                              >  1:2637920/1‑68 (MQ=255)
                    aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                              >  1:1122547/1‑68 (MQ=255)
                    aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                              >  1:2133685/1‑68 (MQ=255)
                    aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                              >  1:177847/1‑68 (MQ=255)
                    aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                              >  1:1728159/1‑68 (MQ=255)
                    aaCCGTTGACGGCGATTACCAGACATTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                              >  1:778783/1‑68 (MQ=255)
                        gTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTaaa                                         <  1:1409904/69‑1 (MQ=255)
                              ggCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGtt                                                              <  1:820222/42‑1 (MQ=255)
                              ggCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATcc                                    >  1:67398/1‑68 (MQ=255)
                              ggCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATcc                                    >  1:3261032/1‑68 (MQ=255)
                                      ccAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTg                                                          <  1:2823902/38‑1 (MQ=255)
                                      ccAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTg                                                          <  1:730432/38‑1 (MQ=255)
                                      ccAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgc                           >  1:1247739/1‑69 (MQ=255)
                                      ccAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgc                           >  1:1979177/1‑69 (MQ=255)
                                               cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAACCCGCAGCACTGGt                  >  1:1192439/1‑69 (MQ=255)
                                               cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGt                  >  1:641985/1‑69 (MQ=255)
                                               cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGt                  >  1:2977258/1‑69 (MQ=255)
                                               cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGt                  >  1:3167529/1‑69 (MQ=255)
                                               cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGt                  >  1:3367919/1‑69 (MQ=255)
                                               cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGt                  >  1:1950930/1‑69 (MQ=255)
                                               cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGt                  >  1:1729074/1‑69 (MQ=255)
                                                  aTCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTg                <  1:2909664/68‑1 (MQ=255)
                                                  aTCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTg                <  1:1725705/68‑1 (MQ=255)
                                                  aTCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTg                <  1:878805/68‑1 (MQ=255)
                                                     gAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAga            <  1:2420895/69‑1 (MQ=255)
                                                     gAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAga            <  1:207365/69‑1 (MQ=255)
                                                               gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc  <  1:3316254/69‑1 (MQ=255)
                                                               gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc  <  1:2678852/69‑1 (MQ=255)
                                                               gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc  <  1:645039/69‑1 (MQ=255)
                                                               gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc  <  1:2324581/69‑1 (MQ=255)
                                                               gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc  <  1:1080736/69‑1 (MQ=255)
                                                               gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc  <  1:129927/69‑1 (MQ=255)
                                                               gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc  <  1:1888847/69‑1 (MQ=255)
                                                               gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc  <  1:2172456/69‑1 (MQ=255)
                                                               gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc  <  1:2453248/69‑1 (MQ=255)
                                                               gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc  <  1:2546544/69‑1 (MQ=255)
                                                               gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc  <  1:2722765/69‑1 (MQ=255)
                                                               gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc  <  1:2913315/69‑1 (MQ=255)
                                                                    cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCtcag                         <  1:2228583/41‑1 (MQ=255)
                                                                    cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgcag                         <  1:3208489/41‑1 (MQ=255)
                                                                    cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgcag                         <  1:2805167/41‑1 (MQ=255)
                                                                    cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgcag                         <  1:201780/41‑1 (MQ=255)
                                                                    cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgcag                         <  1:2812195/41‑1 (MQ=255)
                                                                    cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgcag                         <  1:1110167/41‑1 (MQ=255)
                                                                    cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgcag                         <  1:846090/41‑1 (MQ=255)
                                                                    cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAgactg         >  1:2877730/1‑57 (MQ=255)
                                                                    |                                                               
CTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAC  >  minE/108687‑108818

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: