Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,508,991 C→T 22.7% G685G (GGC→GGT yicI → predicted alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,508,9910CT22.7% 162.3 / 36.9 118G685G (GGC→GGTyicIpredicted alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (66/25);  new base T (16/11);  total (82/36)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 118 total aligned reads displayed.

CGTGATAACACTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAGTGCCCGCT  >  minE/2508926‑2509051
                                                                 |                                                            
cGTGATAACACTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCAc                                                            <  1:1483466/68‑1 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:2566665/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:2821630/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:1822817/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        <  1:186877/69‑1 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        <  1:484044/69‑1 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:2173725/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:2204001/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        <  1:2220210/69‑1 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:615236/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:481158/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:1442416/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:1352632/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:1313341/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:732050/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:1235955/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        <  1:1205169/69‑1 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:3069357/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        <  1:1055499/69‑1 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        >  1:1072849/1‑69 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        <  1:107203/69‑1 (MQ=255)
   gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATAACGCGTGGCACGAAGGTACGGCa                                                        >  1:640421/1‑69 (MQ=255)
    aTAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        <  1:2639065/68‑1 (MQ=255)
    aTAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa                                                        <  1:1182204/68‑1 (MQ=255)
    aTAACACTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGc                                                         >  1:1623720/1‑67 (MQ=255)
       acacTCTAGTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc                                                    >  1:2318597/1‑69 (MQ=255)
       acacTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCa                                                             >  1:3137831/1‑60 (MQ=255)
       acacTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc                                                    >  1:1571830/1‑69 (MQ=255)
       acacTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc                                                    >  1:3259043/1‑69 (MQ=255)
       acacTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc                                                    >  1:3072350/1‑69 (MQ=255)
       acacTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc                                                    >  1:2325762/1‑69 (MQ=255)
       acacTCTACTGGCGCTGGGAAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc                                                    >  1:3148172/1‑69 (MQ=255)
            ctACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACctc                                               >  1:2414048/1‑69 (MQ=255)
               cTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTc                                            >  1:729121/1‑69 (MQ=255)
                 ggCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACGTCTTCaa                                          <  1:2152507/69‑1 (MQ=255)
                 ggCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCaa                                          <  1:2906974/69‑1 (MQ=255)
                 ggCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCaa                                          <  1:2608217/69‑1 (MQ=255)
                 ggCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCaa                                          <  1:3253813/69‑1 (MQ=255)
                 ggCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCaa                                          <  1:1541148/69‑1 (MQ=255)
                 ggCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCaa                                          <  1:3197011/69‑1 (MQ=255)
                        ggCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCaa                                   >  1:910802/1‑69 (MQ=255)
                              aacGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:1890866/1‑69 (MQ=255)
                              aacGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:2576321/1‑69 (MQ=255)
                              aacGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:416135/1‑69 (MQ=255)
                              aacGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:2317808/1‑69 (MQ=255)
                              aacGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:3244043/1‑69 (MQ=255)
                              aacGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:537769/1‑69 (MQ=255)
                              aacGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:2771416/1‑69 (MQ=255)
                              aacGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:2026996/1‑69 (MQ=255)
                              aacGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:1683683/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:2686236/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:343572/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:2107909/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:2504865/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:387121/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:421870/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:1110720/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:129250/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:2131566/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:1673159/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:1862198/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:2574630/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:515431/1‑69 (MQ=255)
                                    cAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:1753485/1‑69 (MQ=255)
                                    cAACGTACCGATTACGTGTGGCACGAAGGCACGGCATTCGACCTCTTCAATCTGCAAGACGGGCATGaa                       >  1:3264324/1‑69 (MQ=255)
                                           ccGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGc                     >  1:2280783/1‑64 (MQ=255)
                                           ccGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCt                >  1:1474011/1‑69 (MQ=255)
                                           ccGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCt                >  1:927763/1‑69 (MQ=255)
                                           ccGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCt                >  1:1022797/1‑69 (MQ=255)
                                              aTTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCaa                                          <  1:936500/40‑1 (MQ=255)
                                               ttACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGa            <  1:2470848/69‑1 (MQ=255)
                                               ttACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGa            <  1:353464/69‑1 (MQ=255)
                                                tACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATg                         <  1:3028824/55‑1 (MQ=255)
                                                tACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATg                         <  1:101775/55‑1 (MQ=255)
                                                tACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATg                         <  1:887360/55‑1 (MQ=255)
                                                tACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATg                         <  1:729830/55‑1 (MQ=255)
                                                tACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATg                         <  1:2431797/55‑1 (MQ=255)
                                                tACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGaa           >  1:1111395/1‑69 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCa                                           >  1:718543/1‑36 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:3064284/1‑50 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:1214567/1‑50 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:2987437/1‑50 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:2952370/1‑50 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:2284997/1‑50 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:2275317/1‑50 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:132044/1‑50 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:1502416/1‑50 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACggg                             >  1:3317357/1‑50 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAg          >  1:1639546/1‑69 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAg          >  1:2891393/1‑69 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAg          >  1:397247/1‑69 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAg          >  1:359224/1‑69 (MQ=255)
                                                 aCGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAg          >  1:2738792/1‑69 (MQ=255)
                                                          aCGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAGTGCCCgct  <  1:1766552/68‑1 (MQ=255)
                                                          aCGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAGTGCCCgct  <  1:2009165/68‑1 (MQ=255)
                                                          aCGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAGTGCCCgct  <  1:1315391/68‑1 (MQ=255)
                                                          aCGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAGTGCCCgct  <  1:1271694/68‑1 (MQ=255)
                                                          aCGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAGTGCCCgct  <  1:815608/68‑1 (MQ=255)
                                                          aCGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAGTGCCCgct  <  1:1000995/68‑1 (MQ=255)
                                                          aCGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAGTGCCCgct  <  1:2665627/68‑1 (MQ=255)
                                                                 |                                                            
CGTGATAACACTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCCACCTCTTCAATCTGCAAGACGGGCATGAAGCCGTCTGTGAAGTGCCCGCT  >  minE/2508926‑2509051

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