Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,852,381 A→G 16.1% L142P (CTC→CCC)  yeaM ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,852,3810A→G16.1% 4.2 ‑4.3 0.044425yeaMpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (18/18):  total (26/18)
Fisher's exact test strand distribution p-value = 1.37e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GGGTAGGTGCAGCTGTTGTTGTGGCTGCTGAGGGAGTTCATCAAACAGGACATGAATTAATCTGCGG  >  REL606/1852347‑1852413
                                  |                                
gggTAGGTGCAGCTGTTGTTGTGGCTGGTTAGGGAg                                 >  1:2070650/1‑36 (MQ=37)
 ggTGGGTGCAGCTGTTGTTGTGGCTGCTGAGGGAGt                                <  1:3734471/36‑1 (MQ=255)
 ggTAGGTGCAGCTGTTGTTGTGGCTGCTGAGGGAGt                                <  1:629296/36‑1 (MQ=255)
  gTAGGTGCAGCTGTTGTTGTGGCTGCTGAGGGAGtt                               <  1:2087980/36‑1 (MQ=255)
  gTAGGTGCAGCTGTTGTTGTGGCTGCTGAGGGAGtt                               <  1:3641824/36‑1 (MQ=255)
   tAGGTGCAGCTGTTGTTGTGGCTGGTGGGGGGGttt                              >  1:3961975/1‑35 (MQ=25)
    aGGTGCAGCTGTTGTTGTGGCTTTTGAGGGGGttta                             >  1:1438561/1‑34 (MQ=21)
     ggTGCAGCTGTTGTTGTGGCTGCTGAGGGAGTtcat                            <  1:3440518/36‑1 (MQ=255)
      gTGCAGCTGTTGTTGTGGCTGGTGGGGGGGTttatc                           >  1:5380054/1‑36 (MQ=19)
       tGCAGCTGTTGTTGTGGCTGGTGAGGGGGGttatca                          >  1:5678237/1‑36 (MQ=19)
       tGCAGCTGTTGTTGTGGCTGCTGGGGGGGTtcatca                          >  1:5198400/1‑36 (MQ=37)
        gCAGCTGTTGTTGTGGCTGGTGTGGGGGTTCCTCaa                         >  1:5110497/1‑36 (MQ=19)
        gCAGCTGTTGTTGTGGCTGGTGAGGGAGTTTATCaa                         >  1:4214831/1‑36 (MQ=38)
        gCAGCTGTTGTTGTGGCTGCTGAGGGGGTTTATTaa                         >  1:1573726/1‑36 (MQ=25)
          aGCTGTTGTTGTGGCTGGTGAGGGGGTTTATCaaaa                       >  1:2568029/1‑35 (MQ=21)
          aGCTGTTGTTGTGGCTGCTGAGGGAGTTCATCAAAc                       <  1:1191668/36‑1 (MQ=255)
           gCTGTTGTTGTGGCTGGTGTGGGAGTTCATCATACa                      >  1:409972/1‑36 (MQ=25)
            cTGTTGTTGTGGCTGCTGAGGGAGTTCATCAAACAg                     <  1:4328656/36‑1 (MQ=255)
             tgttgttgtGGCTGCTGAGGGAGTTCATCAAACAgg                    >  1:4136809/1‑36 (MQ=255)
                tgttgtGGCTGCTTAGGGAGTTCATCAAACAGGACa                 >  1:1777374/1‑36 (MQ=255)
                tgttgtGGCTGCTGAGGGAGTTCATCAAACAGGACa                 >  1:2862128/1‑36 (MQ=255)
                 gttgtGGCTGCTGTGGGAGTTCATCAAACAGGGCAt                >  1:1633766/1‑36 (MQ=37)
                  ttgtGGCTGATGAGGGAGTTCATCAAACAGGACATg               >  1:425403/1‑36 (MQ=255)
                   tgtgGCTGCTGAGGGAGTTCATCAAACAGGACATGa              <  1:3124737/36‑1 (MQ=255)
                   tgtgGCTGCTGAGGGAGTTCATCAAACAGGACATGa              >  1:3832039/1‑36 (MQ=255)
                    gtgGCTGCTGAGGGAGTTCATCAAACAGGACATGaa             >  1:2385929/1‑36 (MQ=255)
                    gtgGCTGCTGAGGGAGTTCATCAAACAGGACATGaa             <  1:4008674/36‑1 (MQ=255)
                     tgGCTGCTGAGGGAGTTCATCAAACAGGACATTAAt            <  1:916224/36‑1 (MQ=255)
                     tgGCTGCTGAGGGAGTTCATCAAACAGGACATGAAt            >  1:4963369/1‑36 (MQ=255)
                      ggCTGCTGAGGGAGTTCATCAAACAGGACATGAAtt           >  1:5696417/1‑36 (MQ=255)
                        ctgctgATGGATTTCATCAAACAGGACATGAATTaa         <  1:4668373/36‑1 (MQ=37)
                        ctgctgAGGGAGTTCATCAAACAGGACATGAATTaa         >  1:2355838/1‑36 (MQ=255)
                        ctgctgAGGGAGTTCATCAAACAGGACATGAATTaa         <  1:3397785/36‑1 (MQ=255)
                         tgctgAGGGAGTTCATCAAACAGGACATGAATTAAt        >  1:2290823/1‑36 (MQ=255)
                         tgctgAGGGAGTTAATAAAACAGGACATGCATTAAt        <  1:5784349/36‑1 (MQ=37)
                          gctgAGGGAGTTCATCAAACAGGACATGAATTAATc       <  1:2239578/36‑1 (MQ=255)
                          gctgAGGGAGTTCATCAAACAGAACATGAATTAATc       <  1:4455118/36‑1 (MQ=255)
                          gaagAGGGAGTTCATCAAACAGGACATGAATTAATc       <  1:3686185/33‑1 (MQ=255)
                            tgAGGGAGTTCATCAAACAGGAGATGAATTTATGTg     >  1:4391883/1‑36 (MQ=25)
                             gAGGGAGTTCATCAAACAGGACATGAATTAATCTGc    <  1:1688136/36‑1 (MQ=255)
                             gAGGGAGTTCATCAAACAGGACATGAATTAATCTGc    >  1:1551910/1‑36 (MQ=255)
                              aGGGAGTTCATCAAACAGGATATGAAGTAATCTGCg   >  1:1413747/1‑36 (MQ=37)
                              aGGGAGTTCATCAAACAGGACATGAATTAATCTGCg   <  1:2812573/36‑1 (MQ=255)
                               gggAGTTCATCAAACAGGACATGAATTAATCTGCgg  >  1:433871/1‑36 (MQ=255)
                                  |                                
GGGTAGGTGCAGCTGTTGTTGTGGCTGCTGAGGGAGTTCATCAAACAGGACATGAATTAATCTGCGG  >  REL606/1852347‑1852413

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 14 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG