Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,508,937 T→C 52.0% T667T (ACT→ACC yicI → predicted alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,508,9370TC52.0% 0.3 / 43.6 52T667T (ACT→ACCyicIpredicted alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base T (13/12);  new base C (16/11);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

AGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCC  >  minE/2508875‑2509001
                                                              |                                                                
aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:1443447/64‑1 (MQ=255)
aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:2326014/64‑1 (MQ=255)
aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:2385026/64‑1 (MQ=255)
aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:320757/64‑1 (MQ=255)
aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:1201576/64‑1 (MQ=255)
aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:785311/64‑1 (MQ=255)
aGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACActc                                                                 <  1:712987/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:1988900/69‑1 (MQ=255)
                     caCGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCaa                                       <  1:1307918/69‑1 (MQ=255)
                     caCGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCaa                                       <  1:197266/69‑1 (MQ=255)
                     caCGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCaa                                       <  1:2614564/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:2614555/69‑1 (MQ=255)
                                                      gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:732050/1‑69 (MQ=255)
                                                      gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:3198155/1‑69 (MQ=255)
                                                      gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:2566665/1‑69 (MQ=255)
                                                      gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:2821630/1‑69 (MQ=255)
                                                      gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:3069357/1‑69 (MQ=255)
                                                      gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:615236/1‑69 (MQ=255)
                                                      gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      <  1:107203/69‑1 (MQ=255)
                                                      gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      <  1:418241/69‑1 (MQ=255)
                                                      gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      >  1:481158/1‑69 (MQ=255)
                                                      gATAGCACCCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGCGTGGCACGAAGGTACGGCa      <  1:484044/69‑1 (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)
                                                              |                                                                
AGTCGCTGGCATAAACAGCAGCACGGCTTCCTGAGTCTGCCCGTTTATGTGCGTGATAACACTCTACTGGCGCTGGGCAACAACGATCAACGTCCCGATTACGTGTGGCACGAAGGCACGGCATTCC  >  minE/2508875‑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: