Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 901,719 A→G 18.6% K522E (AAA→GAA)  ycjM → predicted glucosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE901,7190AG18.6% 170.9 / 26.7 112K522E (AAA→GAA) ycjMpredicted glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (22/69);  new base G (13/8);  total (35/77)
Fisher's exact test for biased strand distribution p-value = 1.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

Only 100 of 115 total aligned reads displayed.

TAATGCGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTGAAGTTGATATATT  >  minE/901657‑901786
                                                              |                                                                   
tAATGCGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGTATTGTTTAATGTCAGTAAAAata                                                               <  1:2200324/69‑1 (MQ=255)
tAATGCGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTTTTTAATGTCAGTAAAAata                                                               <  1:1928486/69‑1 (MQ=255)
tAATGCGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAata                                                               <  1:899964/69‑1 (MQ=255)
tAATGCGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAata                                                               <  1:2191229/69‑1 (MQ=255)
tAATGCGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAata                                                               <  1:2352742/69‑1 (MQ=255)
tAATGCGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAata                                                               <  1:550077/69‑1 (MQ=255)
tAATGCGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAata                                                               <  1:554066/69‑1 (MQ=255)
tAATGCGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAata                                                               <  1:2300183/69‑1 (MQ=255)
tAATGCGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAata                                                               <  1:2395222/69‑1 (MQ=255)
     cGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTaaa                                                                   >  1:3416034/1‑60 (MQ=255)
     cGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTaaa                                                                   >  1:1676265/1‑60 (MQ=255)
     cGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTaaa                                                                   >  1:2833254/1‑60 (MQ=255)
     cGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAg                                                          <  1:1387624/69‑1 (MQ=255)
       tATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGc                                                         <  1:2917594/68‑1 (MQ=255)
       tATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGc                                                         <  1:2990787/68‑1 (MQ=255)
       tATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGc                                                         <  1:1280610/68‑1 (MQ=255)
       tATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGc                                                         <  1:2719725/68‑1 (MQ=255)
       tATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGc                                                         <  1:2339370/68‑1 (MQ=255)
       tATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGc                                                         <  1:1144058/68‑1 (MQ=255)
       tATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGc                                                         <  1:2859798/68‑1 (MQ=255)
       tATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCa                                                        >  1:3157429/1‑69 (MQ=255)
       tATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCa                                                        >  1:948426/1‑69 (MQ=255)
           ccgagaagTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTa                                                    <  1:1534611/66‑1 (MQ=255)
           ccAAGAAGTAACGCTGATGGTAATTGTCTGACTGTATTGTTTAATGTCAGTAAAAATATTCAGCATGTa                                                    <  1:2464085/69‑1 (MQ=255)
           ccAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTa                                                    >  1:282750/1‑69 (MQ=255)
           ccAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTa                                                    <  1:789588/69‑1 (MQ=255)
           ccAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTa                                                    <  1:2228063/69‑1 (MQ=255)
           ccAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTa                                                    >  1:1633336/1‑69 (MQ=255)
           ccAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTa                                                    <  1:702499/69‑1 (MQ=255)
           ccAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTa                                                    >  1:282749/1‑69 (MQ=255)
           ccAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTa                                                    <  1:352870/69‑1 (MQ=255)
           ccAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTa                                                    <  1:3049460/69‑1 (MQ=255)
           ccAAGAAGTAACGCTGATGGTAATTGTCTGACTGGAGTGTTTAATGTCAGTAAAAATATTCAGCATGTa                                                    <  1:2253263/69‑1 (MQ=255)
              agaagTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:3027953/68‑1 (MQ=255)
              agaagTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:2276867/68‑1 (MQ=255)
              agaagTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:711384/68‑1 (MQ=255)
              agaagTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:1229782/68‑1 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTc                                                            >  1:349750/1‑57 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCa                                                           >  1:967253/1‑58 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCa                                                           >  1:261639/1‑58 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCa                                                           >  1:2900326/1‑58 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCa                                                           >  1:2920687/1‑58 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCa                                                           >  1:3076549/1‑58 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCa                                                           >  1:3098387/1‑58 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCa                                                           >  1:421650/1‑58 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCa                                                           >  1:518456/1‑58 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCa                                                           >  1:861833/1‑58 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCa                                                           >  1:1920880/1‑58 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCa                                                           >  1:2177321/1‑58 (MQ=255)
               gaagCAATGCGGATGGCGATTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCa                                                           >  1:2324636/1‑58 (MQ=255)
                    aaCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  >  1:2704269/1‑62 (MQ=255)
                    aaCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  >  1:1027670/1‑62 (MQ=255)
                    aaCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  >  1:740964/1‑62 (MQ=255)
                    aaCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  >  1:452090/1‑62 (MQ=255)
                    aaCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  >  1:2184542/1‑62 (MQ=255)
                    aaCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  >  1:2925104/1‑62 (MQ=255)
                    aaCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  >  1:2842049/1‑62 (MQ=255)
                    aaCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  >  1:250704/1‑62 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:3390890/61‑1 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:914121/61‑1 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:1458063/61‑1 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:2712389/61‑1 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:1643219/61‑1 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  >  1:1799801/1‑61 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  >  1:1804418/1‑61 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:1815335/61‑1 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:2032583/61‑1 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:213942/61‑1 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:3083583/61‑1 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTAAATATTACTa                                          <  1:2887916/69‑1 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTAAATATTACTa                                          <  1:1510697/69‑1 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTAAATATTACTa                                          <  1:1795583/69‑1 (MQ=255)
                     aCGCTGATGGTAATTGTCTGACTGGAGTGTTTAATGTCAGTAAAAATATTCAGCATGTAAATATTACTa                                          <  1:1114301/69‑1 (MQ=255)
                      cGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTAAATATTACTa                                          <  1:2795767/68‑1 (MQ=255)
                      cGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTAAATATTACTa                                          <  1:3013356/68‑1 (MQ=255)
                       gCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTa                                                    <  1:1082329/57‑1 (MQ=255)
                       gCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTAAATATTACTaa                                         <  1:2227725/68‑1 (MQ=255)
                       gCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTAAATATTACTaa                                         <  1:2369433/68‑1 (MQ=255)
                       gCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTAAATATTACTaa                                         <  1:465011/68‑1 (MQ=255)
                        cTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGc                                                         <  1:337997/51‑1 (MQ=255)
                                gaTTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCAGCATATAAATATTACTGACCTGCACGGt                               <  1:2355216/68‑1 (MQ=255)
                                gaTTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCAGCATATAAATATTACTGACCTGCACGGt                               <  1:248278/68‑1 (MQ=255)
                                gaTTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCAGCATATAAATATTACTGACCTGCACGGt                               <  1:632094/68‑1 (MQ=255)
                                gaTTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCAGCATATAAATATTACTGACCTGCACGGt                               <  1:621361/68‑1 (MQ=255)
                                gaTTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCAGCATATAAATATTACTGACCTGCACGGt                               <  1:2970557/68‑1 (MQ=255)
                                gaTTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCAGCATATAAATATTACTGACCTGCACGGt                               <  1:2504374/68‑1 (MQ=255)
                                gaTTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCAGCATATAAATATTACTGACCTGCACGGt                               <  1:2767140/68‑1 (MQ=255)
                                gaTTGCCTGACGGGATTGTTTAATGTCAGTGAAAATATTCAGCATATAAATATTACTGACCTGCACGGt                               <  1:3192213/68‑1 (MQ=255)
                                   tGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTaaa                                                  <  1:3025786/47‑1 (MQ=255)
                                                tGTTTAATGTCAGTAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTg               <  1:2238664/69‑1 (MQ=255)
                                                      aTGTCAGTAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTGAAGtt          <  1:2214911/68‑1 (MQ=255)
                                                      aTGTCAGTAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTGAAGtt          <  1:310300/68‑1 (MQ=255)
                                                      aTGTCAGTAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTGAAGtt          <  1:2800305/68‑1 (MQ=255)
                                                        gTCAGTAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTGAAGTTGa        <  1:468743/68‑1 (MQ=255)
                                                        gTCAGTAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTGAAGTTGa        <  1:2107999/68‑1 (MQ=255)
                                                        gTCAGTAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTGAAGTTGa        <  1:1279660/68‑1 (MQ=255)
                                                        gTCAGTAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTGAAGTTGa        <  1:2140810/68‑1 (MQ=255)
                                                         tCAGTAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTGAAGTTGata      <  1:2490092/69‑1 (MQ=255)
                                                             tAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTGAAGTTGATATAtt  >  1:3374020/1‑69 (MQ=255)
                                                             tAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTGAAGTTGATATAtt  >  1:2138323/1‑69 (MQ=255)
                                                              |                                                                   
TAATGCGTATTCCAAGAAGTAACGCTGATGGTAATTGTCTGACTGGATTGTTTAATGTCAGTAAAAATATTCAGCATGTAAATATTACTAATCTGCATGGTCGGGATCTGATTAGTGAAGTTGATATATT  >  minE/901657‑901786

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: