Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 475,288 +G 100% 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.→G100.0% 98.2 28coding (18/1677 nt)ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  new base (9/19):  ref base (0/0):  total (9/19)

GCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGC.GGGGTGGCGTGATGCATTCCGTCTCCTTTT  >  REL606/475253‑475318
                                    |                              
gCGAGCACAAGGCCGCCTACAATGGTGGGGTTAAGc                                 >  1:118373/1‑36 (MQ=37)
gCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGc                                 <  1:5305072/36‑1 (MQ=255)
gCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGc                                 >  1:4478893/1‑36 (MQ=255)
 cGAGCACAAGGCCGCCAACAATGGTGGCGATAAGC.g                               >  1:6717143/1‑36 (MQ=255)
  gAGCACAAGGCCGCCAACAATGGTGGTGATAAGC.gg                              >  1:7695089/1‑36 (MQ=255)
  gAGCACAAGGCCGCCAACAATGGTGGTGATAAGC.gg                              <  1:1379578/36‑1 (MQ=255)
  gAGCACAAGGCCGCCAACAATGGTGGTGATAAGC.gg                              >  1:5050456/1‑36 (MQ=255)
  gAGCACAAGGCCGCCAACAATGGTGGTGATAAGC.gg                              >  1:4654697/1‑36 (MQ=255)
  gAGCACAAGGCCGCCAACAATGGCGGAGATAAGC.gg                              >  1:7602369/1‑36 (MQ=37)
   aGCACAAGGCCGCCAACAATGGTGGTGATAAGC.ggg                             >  1:677087/1‑36 (MQ=255)
   aGCACAAGGCCGCCAACAATGGTGGTGATAAGC.ggg                             >  1:1746207/1‑36 (MQ=255)
   aGCACAAGGCCGCCAACAATGGTGGTGATAAGC.ggg                             >  1:5278529/1‑36 (MQ=255)
   aGCACAAGGCCGCCAACAATGGTGGTGATAAGC.ggg                             <  1:3016363/36‑1 (MQ=255)
    gCACAAGGCCGCCAACAATGGTGGTGATAAGC.gggg                            >  1:5888412/1‑36 (MQ=255)
     cacaAGGCCGCCAACAATGGTGGTGATAAGC.ggggg                           >  1:5321210/1‑35 (MQ=255)
          ggCCGCCAACAATGGTGGTGATAAGCGGGGGTGGCg                       <  1:4641332/36‑1 (MQ=37)
             cgccAACAATGGTGGTGATAAGCGGGGGTGGCGTGa                    <  1:2315555/36‑1 (MQ=37)
             cgccAACAATGGTGGTGATAAGCGGGGGTGGCGTGa                    <  1:568100/36‑1 (MQ=37)
             cgccAACAATGGTGGTGATAAGCGGGGGTGGCGTGa                    <  1:7220341/36‑1 (MQ=37)
             cgccAACAATGGTGGTGATAAGCGGGGGTGGCGTGa                    <  1:4214311/36‑1 (MQ=37)
             cgccAACAATGGTGGTGATAAGCGGGGGTGGCGTGa                    <  1:397318/36‑1 (MQ=37)
              gccAACAATGGTGGTGATAAGCGGGGGTGGCGTGAt                   >  1:1351220/1‑36 (MQ=37)
              gccAACAATGGTGGTGATAAGCGGGGGTGGCGTGAt                   <  1:5274622/36‑1 (MQ=37)
                caacaaTGGTGGTGATAAGCGGGGGTGGCGTGATGc                 <  1:2581644/36‑1 (MQ=37)
                 aacaaTGGCGGTGATAAGCGGGGGTGGCGTGATGCa                >  1:6997256/1‑36 (MQ=21)
                  acaaTGGTGGTGATAAGCGGGGGTGGCGTGATGCAt               <  1:4434932/36‑1 (MQ=37)
                  acaaTGGTGGTGATAAGCGGGGGTGGCGTGATGCAt               <  1:3635593/36‑1 (MQ=37)
                    aaTGGTGGTGATAAGCGGGGGTGGCGTGATGCATTc             >  1:3299874/1‑36 (MQ=37)
                    aaTGGTGGTGATAAGCGGGGGTGGCGTGATGCATTc             <  1:1613895/36‑1 (MQ=37)
                     aTGGTGGTGATAAGCGGGGGTGGCGTGTTGCTTTcc            >  1:5274984/1‑36 (MQ=16)
                     aTGGTGGTGATAAGCGGGGGTGGCGTGATGCATTcc            >  1:2529725/1‑36 (MQ=37)
                      tggcggTGATAAGCGGGGGTGGCGTGATGCATTCCg           <  1:7061471/36‑1 (MQ=21)
                      gggtggtGATAAGCGGGGGTGGCGTGATGCATTCCg           <  1:78337/35‑1 (MQ=37)
                       ggtggtGATAAGCGGGGGTGGCGTGATGCATTCCGt          <  1:3067211/36‑1 (MQ=37)
                         tggtgATAAGCGGGGGTGTCGTGATGCATTCCGtct        <  1:6988261/36‑1 (MQ=25)
                          ggtgATAAGCGGGGGTGGCGTGATGCATTCCGtctc       <  1:5548131/36‑1 (MQ=37)
                          ggtgATAAGCGGGGGTGGCGTGATGCATTCCGtctc       <  1:2904772/36‑1 (MQ=37)
                           gtgATAAGCGGGGGTGGCGTGATGCATTCCGtctca      >  1:4825783/1‑35 (MQ=37)
                           gtgATAAGCGGGGGTGGCGTGATGCATTCCGTCTcc      <  1:2307078/36‑1 (MQ=37)
                           gtgATAAGCGGGGGTGGCGTGATGCATTCCGTCTcc      >  1:29982/1‑36 (MQ=37)
                              aTAAGCGGGGGTGGCGTGATGCATTCCGTCTCCttt   >  1:6387958/1‑36 (MQ=37)
                              aTAAGCGGGGGTGGCGTGATGCATTCCGTCTCCttt   <  1:2347662/36‑1 (MQ=37)
                               tAAGCGGGGGCGGCGTGATGCATTCCGTCTCCtttt  >  1:1639869/1‑36 (MQ=21)
                                    |                              
GCGAGCACAAGGCCGCCAACAATGGTGGTGATAAGC.GGGGTGGCGTGATGCATTCCGTCTCCTTTT  >  REL606/475253‑475318

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG