Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,454,742 +G 5.3% coding (841/927 nt) yfeR ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,454,7421.→G 5.3% 3.7 ‑2.2 ‑0.238coding (841/927 nt)yfeRpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (11/25):  total (13/25)
Fisher's exact test strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.02e-01

CTGGTCTCGTACCACGTCCCATAACGCTTCAGCGGC.GGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCA  >  REL606/2454707‑2454777
                                    |                                   
cTGGTCTCGTACCTCGTCCCATAACGCTTCAgcggc                                      <  1:5235941/36‑1 (MQ=255)
cTGGTCTCGTACCACGTCCCATAACGCTTCAgcggc                                      <  1:6967664/36‑1 (MQ=255)
cTGGTCTCGTACCACGTCCCATAACGCTTCAgcggc                                      <  1:4986504/36‑1 (MQ=255)
  ggTCTCGTACCACGTCCCATAACGCTTCAGcggc.gg                                   >  1:2336986/1‑36 (MQ=255)
   gTCTCGTACCACGTCCCATAACGCTTCAGcggc.ggc                                  >  1:5467560/1‑35 (MQ=255)
   gTCTCGTACCACGTCCCATAACGCTTCAGCGGC.GGt                                  <  1:898336/36‑1 (MQ=255)
   gTCTCGTACCACGTCCCATAACGCTTCAGCGGC.GGt                                  <  1:7737486/36‑1 (MQ=255)
     ctcGTACCACGTCCCATAACGCTTCAGCGGC.GGTgg                                <  1:6177638/36‑1 (MQ=255)
      tcGTACCACGTCCCATAGCGCTTCAGCGGC.GGTggc                               <  1:323310/36‑2 (MQ=39)
      tcGTACCACGTCCCATAACGCTTCAGCGGC.GGTGGa                               <  1:4174034/36‑1 (MQ=255)
      tcGTACCACGTCCCATAACGCTTCAGCGGC.GGTGGa                               <  1:6345545/36‑1 (MQ=255)
      tcGTACCACGTCCCATAACGCTTCAGCGGC.GGGggg                               >  1:6263372/1‑35 (MQ=255)
        gTACCACGTCCCATAACGCTTCAGCGGC.GGTGGaaa                             <  1:3377898/36‑1 (MQ=255)
        gTACCACGTCCCATAACGCTTCAGCGGC.GGGGGgaa                             >  1:7037617/1‑36 (MQ=37)
         tACCACGTCCCATAACGCTTCAGCGGC.GGTGGGAAg                            >  1:7387036/1‑36 (MQ=255)
         tACCACGTCCCATAACGCTTCAGCGGC.GGTGGAAAg                            <  1:7477612/36‑1 (MQ=255)
          aCCACGTCCCATAACGCTTCAGCGGC.GGTGGgaaag                           >  1:3088092/1‑34 (MQ=38)
          aCCACGTCCCATAACACTTCAGCGGC.GGGGGGAAgg                           >  1:6125073/1‑35 (MQ=25)
           ccACGTCCCATAACGCTTCAGCGGC.GGTGGGAAggg                          >  1:1500118/1‑34 (MQ=38)
            cACGTCCCATAACGCTTCAGCGGC.GGGGGGAAGTgg                         >  1:3852294/1‑35 (MQ=37)
             aCGTCCCATAACGCTTCAGCGGC.GGTGGAAAGTGaa                        <  1:2509913/36‑1 (MQ=255)
             aCGTCCCATAACGCTTCAGCGGC.GGTGGAAAGTGaa                        <  1:576241/36‑1 (MQ=255)
               gTCCCATAACGCTTCAGCGGC.GGGGGGAAATGaaag                      >  1:5896601/1‑34 (MQ=21)
                 cccATAACGCTTCAGCGGC.GGTGGAAAGTGAACGGt                    <  1:3628667/36‑1 (MQ=255)
                 cccATAACGCTTCAGCGGC.GGTGGAAAGTGAACGGt                    <  1:1835897/36‑1 (MQ=255)
                  ccATAACGCTTGAGCGGC.GGTGGAAAGTGAACGGtt                   <  1:7421022/36‑1 (MQ=255)
                  ccATAACGCTTCAGCGGC.GGTGGAAAGTGAACGGtt                   <  1:7799941/36‑1 (MQ=255)
                  ccATAACGCTTCAGCGGC.GGTGGAAAGTGAACGGtt                   <  1:582740/36‑1 (MQ=255)
                   cATAACGCTTCAGCGGC.GGTGGAAAGTGAACGGttc                  <  1:6151636/36‑1 (MQ=255)
                    aTAACGCTTCAGCGGC.GGTGGGAAGGGAAAGGtttt                 >  1:2142059/1‑34 (MQ=21)
                     tAACGCTTCAGCGGC.GGTGGAAAGTGAACGGttctt                <  1:3012683/36‑1 (MQ=255)
                     tAACGCTTCAGCGGC.GGTGGAAAGTGAACGGttctt                <  1:7271007/36‑1 (MQ=255)
                      aaCGCTTCAGCGGC.GGGGGGAAGTGAACGGgttttc               >  1:5873597/1‑36 (MQ=19)
                       aCGCTTCAGCGGC.GGTGGGAAGAGAAAGGttttttc              >  1:6363755/1‑34 (MQ=16)
                        cGCTTCAGCGGCGGGGGGAAAGGGACGGGttcttct              >  1:4105876/1‑35 (MQ=9)
                        cGCTTCAGCGGC.GGTGGAAAGTGAACGGTTCTTCCt             <  1:1899618/36‑1 (MQ=255)
                         gCTTCAGCGGCGGGGGGAGAGGGAAGGGTTTTTccg             >  1:3116103/1‑35 (MQ=2)
                         gCTTCAGCGGC.GGTGGAAAGTGAACGGTTCTTCCTg            <  1:3365655/36‑1 (MQ=255)
                          ctgcAGCGGC.GGTGGAATGTGAACGGTTCTTCCTgc           <  1:3985035/33‑1 (MQ=37)
                           ttCAGCGGC.GGTGGAAAGTGAACGGTTCTTCCTgcg          <  1:3144461/36‑1 (MQ=255)
                           ttCAGCGGC.GGTGGAAAGTGAACGGTTCTTCCTgcg          <  1:7309916/36‑1 (MQ=255)
                             cAGCGGA.GGTGGAACGTGAACGGTTCTTCCTGCGCa        <  1:2813121/36‑1 (MQ=37)
                              aGCGGC.GGTGGAAAGTGAACGGTTCTTCCTGCGCAc       <  1:2865195/36‑1 (MQ=255)
                               gcggc.gGTGGGAAGTGAAAGGGTTTTCCTGCGCGcc      >  1:7744588/1‑36 (MQ=16)
                                 ggc.ggTGGAAAGTGAACGGTTCTTCCTGCGCACCAg    >  1:4661234/1‑36 (MQ=255)
                                  gc.ggTGGAAAGTGAACGGTTCTTCCTGCGCACCAGc   <  1:2479701/36‑1 (MQ=255)
                                  gc.ggTGGAAAGTGAACGGTTCTTCCTGCGCACCAGc   >  1:1537336/1‑36 (MQ=255)
                                  gc.ggTGGAAAGTGAACGGTTCTTCCTGCGCACCAGc   >  1:5976149/1‑36 (MQ=255)
                                   c.ggTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCa  >  1:1561282/1‑36 (MQ=255)
                                    |                                   
CTGGTCTCGTACCACGTCCCATAACGCTTCAGCGGC.GGTGGAAAGTGAACGGTTCTTCCTGCGCACCAGCA  >  REL606/2454707‑2454777

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG