Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 475,288 +G 92.1% coding (18/1677 nt) ybaL ← predicted transporter with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606475,2881.→G92.1% 5.4 ‑2.5 ‑0.738coding (18/1677 nt)ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  new base (21/14):  ref base (0/3):  total (21/17)
Fisher's exact test strand distribution p-value = 8.06e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.85e-01

GCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGC.GGG.G.TGGCGTGATGCATTCCGTCTCCTTTT  >  REL606/475253‑475318
                                    |                                
gagAGCACAAGGCCGAAAACAATGGTGGTGAGAAGc                                   >  1:5142946/3‑36 (MQ=21)
gCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGc                                   <  1:6565165/36‑1 (MQ=255)
gCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGc                                   <  1:5550914/36‑1 (MQ=255)
  gAGCACAAGGCCGCCAACAATGGTGGTGATAAGC.gg                                <  1:1805138/36‑1 (MQ=255)
   aacacaAGGCCGCCAACAATGGGGGTGATAAGC.ggg                               >  1:5476873/3‑36 (MQ=38)
    gCACAAGGCCGCCAACAATGGTGGTGATAAGC.ggg.g                             >  1:4973954/1‑36 (MQ=255)
    gCACAAGGCCGCCAACAATGGTGGTGATAAGC.ggg.g                             >  1:691971/1‑36 (MQ=255)
    gCACAAGGCCGCCAACAATGGTGGTGATAAGC.ggg.g                             >  1:5272184/1‑36 (MQ=255)
    gCAAAAGGCCGCCAACAATGGTGGTGATAAGC.ggg.g                             <  1:2362228/36‑1 (MQ=255)
     cacaAGGCCGCCAACAATGGTGGTGATAAGC.ggg.gg                            >  1:7076796/1‑35 (MQ=255)
      acaAGGCCGCCAACAATGGTGGTGATAAGC.ggg.ggt                           <  1:1905281/36‑3 (MQ=255)
      acaAGGCCGCCAACAATGGTGGTGATAAGC.ggg.ggt                           <  1:2375662/36‑3 (MQ=255)
       caAGGCCGCCAACAATGGTGGTGATAAGC.ggg.ggtg                          <  1:3365565/36‑4 (MQ=255)
       caAGGCCGCCAACAATGGTGGTGATAAGC.ggg.ggtg                          <  1:650518/36‑4 (MQ=255)
          ggCCGCCAACAATGTTGGTGATAAGC.GGGTG.TGGCg                       <  1:1495238/36‑1 (MQ=21)
           gccgccAACAATGGTGGTGATAAGC.GGG.G.TGGCGTg                     <  1:1343609/36‑1 (MQ=255)
             cgccAACAATGGTGGTGATAAGCGGGG.G.TGGCGTGa                    <  1:6515538/36‑1 (MQ=37)
             cgccAACAATGGTGGTGATAAACGGGG.G.TGGCGTGa                    >  1:4354824/1‑36 (MQ=25)
              gccAACAATGGTGGTGATAAGCGGGG.G.TGGCGTGAt                   <  1:2770024/36‑1 (MQ=37)
              gccAACAATGGTGGTGATAAGCGGGG.G.TGGCGTGAt                   >  1:5736750/1‑36 (MQ=37)
              gccAACAATGGTGGTGATAAGCGGGG.G.TGGCGTGAt                   >  1:5966941/1‑36 (MQ=37)
                caacaaTGGTGGTGATAAGCGGGG.G.TGGCGTGATGc                 >  1:70535/1‑36 (MQ=37)
                caacaaTGGTGGTGATAAGCGGGG.G.TGGCGTGATGc                 <  1:5264209/36‑1 (MQ=37)
                caacaaTGGTGGTGATAAGCGGGG.G.CGGCGTGATGc                 >  1:1183998/1‑36 (MQ=21)
                 aacaaTGGTGGTGATAAGCGGGG.G.TGGCGTGATGCa                >  1:1994616/1‑36 (MQ=37)
                  acaaTGGTGGTGATAAGCGGGG.G.TGGCGTGATGCAt               <  1:6460633/36‑1 (MQ=37)
                   gaaTGGTGGTGATAAGCGGGG.G.TGGCGTGATGCAtt              <  1:4461230/35‑1 (MQ=25)
                    aaTGGTGGTGATAAGCGGGG.G.GGGCGTGATGCATTc             >  1:3093820/1‑36 (MQ=21)
                     aTGGTGGTGATAAGCGGGG.G.TGGCGTGATGCATTcc            >  1:7034463/1‑36 (MQ=37)
                     aTGGTGGGGATAAGCGGGG.G.CGGCGTGATGCATTcc            >  1:4959795/1‑36 (MQ=16)
                      tggtggTGATAAGCGGGG.G.TGGCGTGATGCATTCCg           <  1:6758436/36‑1 (MQ=37)
                       ggtggtGATAAGCGGGG.G.TGGCGTGATGCATTCCGt          >  1:4055059/1‑36 (MQ=37)
                       ggtggtGATAAGCGGGG.G.TGGCGTGATGCATTCCGt          >  1:1856147/1‑36 (MQ=37)
                        gtggtgATAAGCGGGG.G.TGGCGTGATGCATTCCGtc         >  1:3646302/1‑36 (MQ=37)
                        gtggtgATAAGCGGGG.G.TGGCGTGATGCATTCCGtc         >  1:2941234/1‑36 (MQ=37)
                        gtggagATAAGCGGGG.G.TGGCGTGATGCATTCCGtc         <  1:537596/36‑1 (MQ=25)
                         tggtgATAAGCGGGG.G.TGGCGTGATGCATTCCGtct        <  1:7579302/36‑1 (MQ=37)
                          ggtgATAAGCGGGG.G.TGGCGTGATGCATTCCGtctc       <  1:5159910/36‑1 (MQ=37)
                          ggtgATAAGCGGGG.G.TGGCGTGATGCATTCCGtctc       >  1:3483192/1‑36 (MQ=37)
                          ggtgATAAGCGGGG.G.TGGCGTGATGCATTCCGtctc       <  1:3341633/36‑1 (MQ=37)
                           gtgATAAGCGGGG.G.TGGCGTGATGCATTCCGTCTcc      <  1:377483/36‑1 (MQ=37)
                           gtgATAAGCGGGG.G.TGGCGTGATGCATTCCGTCTcc      >  1:7826292/1‑36 (MQ=37)
                            tgATAAGCGGGG.G.TGGCGTGATGCATTCCGTCTCCt     >  1:6170120/1‑36 (MQ=37)
                            tgATAAGCGGGG.G.TGGCGTGATGCATTCCGTCTCCt     <  1:6793506/36‑1 (MQ=37)
                            tgATAAGAGGGG.G.TGGAGTGATGCATTCCGTCTCCt     >  1:7381885/1‑36 (MQ=19)
                             gATAAGCGGGG.G.TGGCGTGATGCATTCCGTCTCCtt    <  1:7876116/36‑1 (MQ=37)
                             gATAAGCGGGG.G.GGGCGTGATGCATTCAGGCTCCtt    >  1:697967/1‑36 (MQ=16)
                              aTAAGCGGGG.G.TGGCGTGATGCATTCCGTCTCCttt   >  1:1782391/1‑36 (MQ=37)
                              aTAAGCGGGG.G.TGGCGTGATGCATTCCGTCTCCttt   >  1:1287061/1‑36 (MQ=37)
                               tAAGCGGGG.G.TGGCGTGATGCATTCCGTCTCCtttt  <  1:7902692/36‑1 (MQ=37)
                               tAAGCGGGG.G.TGGCGTGATGCATTCCGTCTCCtttt  >  1:5249600/1‑36 (MQ=37)
                               tAAGC.GGG.GTTGGGGTGATGCATTCCGTCTCCtttt  <  1:6400857/36‑1 (MQ=21)
                                    |                                
GCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGC.GGG.G.TGGCGTGATGCATTCCGTCTCCTTTT  >  REL606/475253‑475318

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG