Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 108,761 T→C 42.6% R335R (CGT→CGC 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,7610TC42.6% 14.1 / 45.2 61R335R (CGT→CGCaceEpyruvate dehydrogenase, decarboxylase component E1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base T (12/23);  new base C (9/17);  total (21/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGCAGA  >  minE/108702‑108825
                                                           |                                                                
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:2637920/1‑68 (MQ=255)
     aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                                     >  1:3367725/1‑68 (MQ=255)
     aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                                     >  1:2318796/1‑68 (MQ=255)
     aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                                     >  1:397446/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)
     aaCCGTTGACGGCGATTACCAGACCTTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                                     >  1:1122547/1‑68 (MQ=255)
     aaCCGTTGACGGCGATTACCAGACATTCAAATCGAAAGATGGTGCGTATGTTCGCGAACACTTCTTCg                                                     >  1:778783/1‑68 (MQ=255)
         gTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTaaa                                                <  1:1409904/69‑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:1979177/1‑69 (MQ=255)
                       ccAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgc                                  >  1:1247739/1‑69 (MQ=255)
                                cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAACCCGCAGCACTGGt                         >  1:1192439/1‑69 (MQ=255)
                                cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGt                         >  1:1729074/1‑69 (MQ=255)
                                cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGt                         >  1:3367919/1‑69 (MQ=255)
                                cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGt                         >  1:3167529/1‑69 (MQ=255)
                                cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGt                         >  1:2977258/1‑69 (MQ=255)
                                cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGt                         >  1:641985/1‑69 (MQ=255)
                                cAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGt                         >  1:1950930/1‑69 (MQ=255)
                                   aTCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTg                       <  1:2909664/68‑1 (MQ=255)
                                   aTCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTg                       <  1:878805/68‑1 (MQ=255)
                                   aTCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTg                       <  1:1725705/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:645039/69‑1 (MQ=255)
                                                gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc         <  1:2913315/69‑1 (MQ=255)
                                                gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc         <  1:2722765/69‑1 (MQ=255)
                                                gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc         <  1:129927/69‑1 (MQ=255)
                                                gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc         <  1:1080736/69‑1 (MQ=255)
                                                gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc         <  1:1888847/69‑1 (MQ=255)
                                                gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc         <  1:2678852/69‑1 (MQ=255)
                                                gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc         <  1:2172456/69‑1 (MQ=255)
                                                gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc         <  1:2546544/69‑1 (MQ=255)
                                                gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc         <  1:2453248/69‑1 (MQ=255)
                                                gCGTATGTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc         <  1:2324581/69‑1 (MQ=255)
                                                     tgTTCGCGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGAc         <  1:2746511/63‑1 (MQ=255)
                                                     cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCtcag                                <  1:2228583/41‑1 (MQ=255)
                                                     cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgcag                                <  1:1110167/41‑1 (MQ=255)
                                                     cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgcag                                <  1:846090/41‑1 (MQ=255)
                                                     cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgcag                                <  1:201780/41‑1 (MQ=255)
                                                     cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgcag                                <  1:3208489/41‑1 (MQ=255)
                                                     cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgcag                                <  1:2805167/41‑1 (MQ=255)
                                                     cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCgcag                                <  1:2812195/41‑1 (MQ=255)
                                                     cGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAgactg                >  1:2877730/1‑57 (MQ=255)
                                                       ttCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGCAGa  <  1:3332940/69‑1 (MQ=255)
                                                       ttCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGCAGa  <  1:2568045/69‑1 (MQ=255)
                                                       ttCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGCAGa  <  1:1943984/69‑1 (MQ=255)
                                                       ttCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGCAGa  <  1:1480567/69‑1 (MQ=255)
                                                       ttCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGCAGa  <  1:1134691/69‑1 (MQ=255)
                                                       ttCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGCAGa  <  1:848241/69‑1 (MQ=255)
                                                       ttCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGCAGa  <  1:851164/69‑1 (MQ=255)
                                                       ttCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGCAGa  <  1:852172/69‑1 (MQ=255)
                                                           |                                                                
AACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGCAGA  >  minE/108702‑108825

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: