Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 348,416 A→G 37.1% L238L (CTA→CTG)  gcl → glyoxylate carboligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE348,4160AG37.1% 37.7 / 43.4 70L238L (CTA→CTG) gclglyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base A (25/19);  new base G (21/5);  total (46/24)
Fisher's exact test for biased strand distribution p-value = 6.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCA  >  minE/348348‑348475
                                                                    |                                                           
aatGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTa                                                             >  1:3133398/1‑69 (MQ=255)
aatGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTa                                                             >  1:3140977/1‑69 (MQ=255)
aatGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTa                                                             >  1:1781248/1‑69 (MQ=255)
aatGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTa                                                             >  1:175969/1‑69 (MQ=255)
aatGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTa                                                             >  1:1230141/1‑69 (MQ=255)
aatGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTa                                                             >  1:794159/1‑69 (MQ=255)
aatGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTa                                                             >  1:1700880/1‑69 (MQ=255)
aatGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTa                                                             >  1:551008/1‑69 (MQ=255)
          cAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATgggctggg                                                   >  1:1125972/1‑69 (MQ=255)
          cAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATgggctggg                                                   >  1:1798060/1‑69 (MQ=255)
          cAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATgggctggg                                                   >  1:2127907/1‑69 (MQ=255)
          cAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATgggctggg                                                   >  1:738456/1‑69 (MQ=255)
               gCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGt                                              >  1:242812/1‑69 (MQ=255)
               gCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGt                                              >  1:453756/1‑69 (MQ=255)
               gCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGt                                              >  1:765871/1‑69 (MQ=255)
               gCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGt                                              >  1:951221/1‑69 (MQ=255)
               gCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGt                                              >  1:951973/1‑69 (MQ=255)
               gCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGt                                              >  1:2702807/1‑69 (MQ=255)
               gCACTCTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGt                                              >  1:2920399/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGt                                              >  1:1168259/1‑68 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:1724019/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:2610891/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:2910793/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:221073/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:2611029/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:2058274/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:204396/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:272599/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:1982662/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:1948823/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:3200473/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:1671971/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:1578300/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:1547480/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:1521453/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:499272/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:629572/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:1075409/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:1017264/1‑69 (MQ=255)
                cACTGTTACAACAGTTTGCTGAACTGAACAGCATTCCGGTGATCCCAACGCTGATGGGCTGGGGCTGTa                                             >  1:752855/1‑69 (MQ=255)
                           caGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGAt                                  <  1:1063744/69‑1 (MQ=255)
                           caGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGAt                                  <  1:3118590/69‑1 (MQ=255)
                           caGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGAt                                  <  1:3047315/69‑1 (MQ=255)
                           caGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGAt                                  <  1:2899930/69‑1 (MQ=255)
                           caGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGAt                                  <  1:1291906/69‑1 (MQ=255)
                           caGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGAt                                  <  1:1303594/69‑1 (MQ=255)
                           caGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGAt                                  <  1:2339839/69‑1 (MQ=255)
                           caGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGAt                                  <  1:1496511/69‑1 (MQ=255)
                            aGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGAt                                  <  1:2075801/68‑1 (MQ=255)
                             gTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGAt                                  <  1:454305/67‑1 (MQ=255)
                               ttGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATg                              >  1:2728333/1‑69 (MQ=255)
                               ttGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATg                              >  1:91891/1‑69 (MQ=255)
                               ttGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATg                              >  1:168819/1‑69 (MQ=255)
                               ttGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATg                              >  1:2947453/1‑69 (MQ=255)
                               ttGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATg                              >  1:2639673/1‑69 (MQ=255)
                                tGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGa                             >  1:1350476/1‑69 (MQ=255)
                                tGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGa                             >  1:603657/1‑69 (MQ=255)
                                       cTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATg                      >  1:681713/1‑69 (MQ=255)
                                           ccAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACt                          <  1:2126627/61‑1 (MQ=255)
                                           ccAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACt                          <  1:928631/61‑1 (MQ=255)
                                           ccAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACt                          <  1:1823917/61‑1 (MQ=255)
                                           ccAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACt                          <  1:1265920/61‑1 (MQ=255)
                                           ccAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACt                          <  1:1282968/61‑1 (MQ=255)
                                           ccAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACt                          <  1:2750998/61‑1 (MQ=255)
                                           ccAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACt                          <  1:2791026/61‑1 (MQ=255)
                                              gCGTTCCGGTGATCCTAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACt                          <  1:957125/58‑1 (MQ=255)
                                                       tGATCCCAACGCTGATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTgggcc      <  1:1149675/69‑3 (MQ=255)
                                                       tGATCCCAACGCTGATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTgggcc      <  1:2985321/69‑3 (MQ=255)
                                                        gATCCCAACGCTGATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTgggcc      <  1:2216154/68‑3 (MQ=255)
                                                           cccAACGCTGATGGGCTGGGGCTGTATCCCGGACGATCATGAACTAATGGCCGGGATGGTgggcc      <  1:465374/65‑3 (MQ=255)
                                                           cccAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCa  >  1:1487412/1‑69 (MQ=255)
                                                           cccAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCa  >  1:402697/1‑69 (MQ=255)
                                                           cccAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCa  >  1:3426247/1‑69 (MQ=255)
                                                           cccAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCa  >  1:2685980/1‑69 (MQ=255)
                                                           cccAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCa  >  1:1951252/1‑69 (MQ=255)
                                                           cccAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCa  >  1:1206713/1‑69 (MQ=255)
                                                            ccAACGCTGATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTgggcc      <  1:1315190/64‑3 (MQ=255)
                                                                  cTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACt                          <  1:80134/38‑1 (MQ=255)
                                                                    |                                                           
AATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCA  >  minE/348348‑348475

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: ‑