Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,649,084 G→A 36.6% V498I (GTT→ATT)  glgB → 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,649,0840GA36.6% 71.1 / 66.7 109V498I (GTT→ATT) glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base G (42/27);  new base A (16/24);  total (58/51)
Fisher's exact test for biased strand distribution p-value = 4.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 118 total aligned reads displayed.

TGGGCTTCTGGTACAAGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTACACTGAAAACTTCG  >  minE/2649016‑2649150
                                                                    |                                                                  
tGGGCTTCTGGTACAAGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCgg                                                                    >  1:2972635/1‑69 (MQ=255)
tGGGCTTCTGGTACAAGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCgg                                                                    >  1:375324/1‑69 (MQ=255)
tGGGCTTCTGGTACAAGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCgg                                                                    <  1:754798/69‑1 (MQ=255)
tGGGCTTCTGGTACAAGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCgg                                                                    >  1:3304981/1‑69 (MQ=255)
tGGGCTTCTGGTACAAGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCgg                                                                    >  1:3386218/1‑69 (MQ=255)
tGGGCTTCTGGTACAAGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCgg                                                                    >  1:1541771/1‑69 (MQ=255)
   gCTTCTGGTACAAGTGGAACCTCGGCTGGATGCATGCCACCCTGGACTACATGAAGCTCGACCCGGttt                                                                 >  1:2487968/1‑69 (MQ=255)
   gCTTCTGGTACAAGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGttt                                                                 >  1:2772190/1‑69 (MQ=255)
            aCAAGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAg                                                         <  1:1406826/68‑1 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACCCCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTatc                                                     >  1:2070044/1‑69 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCgg                                                                    >  1:3090772/1‑54 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCg                                                             >  1:212545/1‑61 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGt                                                        >  1:1951262/1‑66 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGt                                                        >  1:2659572/1‑66 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGt                                                        >  1:2803999/1‑66 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGt                                                        >  1:2655599/1‑66 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGt                                                        >  1:87133/1‑66 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGt                                                        >  1:1025251/1‑66 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTatc                                                     >  1:2903591/1‑69 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTatc                                                     >  1:698182/1‑69 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTatc                                                     >  1:727626/1‑69 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTat                                                      <  1:206622/68‑1 (MQ=255)
               aGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTat                                                      <  1:2199516/68‑1 (MQ=255)
                  ggAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTatcatc                                                  >  1:1670853/1‑69 (MQ=255)
                  ggAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTatcatc                                                  >  1:2788351/1‑69 (MQ=255)
                  ggAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTatcatc                                                  >  1:179437/1‑69 (MQ=255)
                        tCGGCTGGATGCATGACCCCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGCTa                                            >  1:636147/1‑69 (MQ=255)
                        tCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATa                                            >  1:340576/1‑69 (MQ=255)
                        tCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATa                                            >  1:2964352/1‑69 (MQ=255)
                        tCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATa                                            >  1:1678370/1‑69 (MQ=255)
                        tCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATa                                            >  1:1313851/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTg                                       >  1:971896/1‑68 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:9814/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:3408114/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:3417824/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:463062/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:46442/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:565166/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:63031/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:963739/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:105939/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:1810558/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:1954395/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:1378126/1‑69 (MQ=255)
                              ggATGCATGACACCCTGGACTACATGAAGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGa                                      >  1:1657572/1‑69 (MQ=255)
                                     tGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCgg                               <  1:1913207/69‑1 (MQ=255)
                                     tGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCgg                               <  1:3136644/69‑1 (MQ=255)
                                     tGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCgg                               <  1:3160742/69‑1 (MQ=255)
                                     tGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCgg                               <  1:3195489/69‑1 (MQ=255)
                                      gACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCgg                               <  1:3293246/68‑1 (MQ=255)
                                       acacCCTGGACTACATGAAGCTCTACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGa                             >  1:1808519/1‑69 (MQ=255)
                                       acacCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCggg                              >  1:742877/1‑68 (MQ=255)
                                       acacCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCggg                              >  1:2404042/1‑68 (MQ=255)
                                       acacCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGa                             >  1:729966/1‑69 (MQ=255)
                                       acacCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGa                             >  1:2365554/1‑69 (MQ=255)
                                       acacCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATAATCACGATAAACTGACCTTCggg                              >  1:3312495/1‑68 (MQ=255)
                                                 cTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTAc                     <  1:1085295/67‑1 (MQ=255)
                                                 cTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTAc                     <  1:100879/67‑1 (MQ=255)
                                                 cTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTAc                     <  1:2130686/67‑1 (MQ=255)
                                                 cTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACa                    >  1:1215160/1‑68 (MQ=255)
                                                    cATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATtctc                        <  1:713484/61‑1 (MQ=255)
                                                    cATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATtctc                        <  1:1413567/61‑1 (MQ=255)
                                                    cATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATtctc                        <  1:2583987/61‑1 (MQ=255)
                                                    cATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATtctc                        <  1:2146677/61‑1 (MQ=255)
                                                    cATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATtctc                        <  1:3075850/61‑1 (MQ=255)
                                                    cATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTa                >  1:3186770/1‑69 (MQ=255)
                                                    cATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTa                >  1:2676881/1‑69 (MQ=255)
                                                    cATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTa                >  1:1856261/1‑69 (MQ=255)
                                                    cATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTa                >  1:63419/1‑69 (MQ=255)
                                                    cATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTa                >  1:253629/1‑69 (MQ=255)
                                                    cATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTa                >  1:1657735/1‑69 (MQ=255)
                                                     aTGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATa                                            >  1:424105/1‑40 (MQ=255)
                                                        aaGCTCTACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:941852/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:127253/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:234052/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:624916/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:1302546/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:127798/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:671302/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:1066328/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:2360292/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:1593221/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:2363514/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:1540802/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:2904281/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:1934525/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:1963122/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:3160810/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:2014989/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:3399767/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:1674260/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:2883809/69‑1 (MQ=255)
                                                        aaGCTCGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:2533891/69‑1 (MQ=255)
                                                            tCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTACACTGaaa        <  1:1724398/69‑1 (MQ=255)
                                                            tCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTACACTGaaa        <  1:3241997/69‑1 (MQ=255)
                                                            tCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTACACTGaaa        <  1:350367/69‑1 (MQ=255)
                                                            tCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTACACTGaaa        <  1:2843199/69‑1 (MQ=255)
                                                            tCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTACACTGaaa        <  1:1817988/69‑1 (MQ=255)
                                                             cGACCCGATTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATGCTCTACAACTACACt            <  1:1019850/64‑1 (MQ=255)
                                                                  cGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTACACTGAAAACTTCg  <  1:2486501/69‑1 (MQ=255)
                                                                    |                                                                  
TGGGCTTCTGGTACAAGTGGAACCTCGGCTGGATGCATGACACCCTGGACTACATGAAGCTCGACCCGGTTTATCGTCAGTATCATCACGATAAACTGACCTTCGGGATTCTCTACAACTACACTGAAAACTTCG  >  minE/2649016‑2649150

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: