Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,508,933 A→G 55.1% N666S (AAC→AGC)  yicI → predicted alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,508,9330AG55.1% 7.9 / 36.4 49N666S (AAC→AGC) yicIpredicted alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (3/19);  new base G (16/11);  total (19/30)
Fisher's exact test for biased strand distribution p-value = 1.34e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GACGGTAGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCC  >  minE/2508869‑2509001
                                                                |                                                                    
gACGGTAGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAacac                                                                   >  1:2634940/1‑68 (MQ=255)
      aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:1443447/64‑1 (MQ=255)
      aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:2385026/64‑1 (MQ=255)
      aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:2326014/64‑1 (MQ=255)
      aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:320757/64‑1 (MQ=255)
      aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:785311/64‑1 (MQ=255)
      aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:712987/64‑1 (MQ=255)
      aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:1201576/64‑1 (MQ=255)
      aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTg                                                            >  1:3104251/1‑69 (MQ=255)
         cGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGgcg                                                         <  1:1433163/69‑1 (MQ=255)
          gCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGgcg                                                         <  1:3197401/68‑1 (MQ=255)
                  aaaCAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGcaac                                                <  1:2038430/69‑1 (MQ=255)
                   aaCAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGcatca                                               >  1:1377159/1‑69 (MQ=255)
                           caCGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCaa                                       <  1:1307918/69‑1 (MQ=255)
                           caCGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCaa                                       <  1:2614564/69‑1 (MQ=255)
                           caCGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCaa                                       <  1:197266/69‑1 (MQ=255)
                           caCGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCaa                                       <  1:1988900/69‑1 (MQ=255)
                           caCGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCaa                                       <  1:2799362/69‑1 (MQ=255)
                                 ttCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCAACGTccc                                 <  1:3216161/69‑1 (MQ=255)
                                                 tttATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACg                 <  1:3223602/69‑1 (MQ=255)
                                                   tATGTGCGTGATAACACTCTACTGGCGCTGGGcaac                                                <  1:131984/36‑1 (MQ=255)
                                                         cGTGATAACACTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCAc          <  1:1483466/68‑1 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      <  1:484044/69‑1 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:732050/1‑69 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:2566665/1‑69 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      <  1:2614555/69‑1 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:481158/1‑69 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:2821630/1‑69 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:3069357/1‑69 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      <  1:418241/69‑1 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:3198155/1‑69 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      <  1:107203/69‑1 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:615236/1‑69 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:1822817/1‑69 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      <  1:1245171/69‑1 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:1313341/1‑69 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:1352632/1‑69 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:1442416/1‑69 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:1235955/1‑69 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      <  1:1205169/69‑1 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      <  1:186877/69‑1 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:1072849/1‑69 (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:1055499/69‑1 (MQ=255)
                                                            gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATAACGCGTGGCACGAAGGTACGGCa      >  1:640421/1‑69 (MQ=255)
                                                             aTAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      <  1:1182204/68‑1 (MQ=255)
                                                             aTAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      <  1:2639065/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:3259043/1‑69 (MQ=255)
                                                                acacTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc  >  1:1571830/1‑69 (MQ=255)
                                                                acacTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc  >  1:3072350/1‑69 (MQ=255)
                                                                acacTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc  >  1:1101852/1‑69 (MQ=255)
                                                                acacTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc  >  1:1944094/1‑69 (MQ=255)
                                                                acacTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc  >  1:740447/1‑69 (MQ=255)
                                                                acacTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc  >  1:2325762/1‑69 (MQ=255)
                                                                acacTCTACTGGCGCTGGGAAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTcc  >  1:3148172/1‑69 (MQ=255)
                                                                |                                                                    
GACGGTAGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCC  >  minE/2508869‑2509001

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: