Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 475,288 +G 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% 101.1 32coding (18/1677 nt)ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  new base (12/18):  ref base (1/1):  total (13/19)

GAGCACAAGGCCGCCAACAATGGTGGTGATAAGC.GG.GGTGGCGTGATGCATTCCGTCTCCTTTTC  >  REL606/475255‑475319
                                  |                                
gAGCACAAGGCCGCCAACTATGGTGGCGCTAAGC.gg                                >  2:2305443/1‑36 (MQ=37)
  gCACAAGGCCGCCAACAATGGTGGTGATAAGC.gg.gg                             <  2:761601/36‑1 (MQ=255)
  gCACAAGGCCGCCAACAATGGTGGTGATAAGC.gg.gg                             >  2:1062864/1‑36 (MQ=255)
   cacaAGGCCGCCAACAATGGTGGTGATAAGC.gg.ggg                            >  2:1749385/1‑35 (MQ=255)
     caAGGCCGCCAACAATGGTGGTGATAAGC.gg.gggtg                          <  2:2467182/36‑4 (MQ=255)
       aGGCCGCCAACAATGGTGGTGATAAGC.GG.GGGgggc                        >  1:2801360/1‑34 (MQ=37)
        ggCCGCCAACAATGGTGGTGATAAGCGGG.GGTGGCg                        >  2:1895108/1‑36 (MQ=37)
        ggCCGCCAACAATGGTGGTGATAAGCGGG.GGTGGCg                        <  1:1881262/36‑1 (MQ=37)
         gccgccAACAATGGTGGTGATAAGCGGG.GGTGGCGt                       <  2:2251437/36‑1 (MQ=37)
            gccAACAATGGTGGTGATAAGCGGG.GGTGGCGTGAt                    <  1:3047435/36‑1 (MQ=37)
             ccAACAATGGTGGTGATAAGCGGG.GGTGGCGTGATg                   <  1:3251921/36‑1 (MQ=37)
             ccAACAATGGTGGTGATAAGCGGG.GGTGGCGTGATg                   <  2:415898/36‑1 (MQ=37)
              caacaaTGGTGGTGATAAGCGGG.GGTGGCGTGATGc                  <  1:2498056/36‑1 (MQ=37)
              caacaaTGGTGGTGATAAGCGGG.GGTGGCGTGATGc                  <  1:2402895/36‑1 (MQ=37)
              caacaaTGGTGGTGATAAGCGGG.GGTGGCGTGATGc                  >  1:1502726/1‑36 (MQ=37)
                acaaTGGTGGTGATAAGCGGG.GGTGGGGTGATGCAt                >  2:2409916/1‑36 (MQ=21)
                acaaTGGTGGTGATAAGCGGG.GGTGGCGTGATGCAt                <  2:3432013/36‑1 (MQ=37)
                acaaTGGTGATGATAAGCGGG.GGTGGCGTGATgact                >  1:3575510/1‑33 (MQ=16)
                 caaTGGTGGTGATAAGCGGG.GGGGGCGTGATGCAtt               >  2:3415632/1‑36 (MQ=21)
                  aaTGGTGGTGATAAGCGGG.GGTGGCGTGATGCATTc              <  1:1509408/36‑1 (MQ=37)
                  aaTGGTGGTGATAAGCGGG.GGTGGCGTGATAGATTc              >  2:2478291/1‑36 (MQ=16)
                   aTGGTGGTGATAAGCGGG.GGTGGCGTGATGCATTcc             <  2:1287664/36‑1 (MQ=37)
                   aTGGTGGTGATAAGCGGG.GGTGGCGTGATGCATTcc             <  1:1708221/36‑1 (MQ=37)
                     ggtggtGATAAGCGGG.GGTGGCGTGATGCATTCCGt           <  2:700197/36‑1 (MQ=37)
                     ggtggtGATAAGCGGG.GGTGGCGTGATGCATTCCGt           <  2:126244/36‑1 (MQ=37)
                      gtggtgATAAGCGGG.GGTGGCGTGATGCATTCCGtc          >  1:2389353/1‑36 (MQ=37)
                       tggagATAAGCGGG.GGTGGCGTGATGCATTCCGtct         <  2:450480/36‑1 (MQ=25)
                        ggtgATAAGCGGG.GGTGGCGTGATGCATTCCGtctc        <  1:2598848/36‑1 (MQ=37)
                         gtgATAAGCGGG.GGTGGCGTGATGCATTTCGTCTcc       >  1:3237855/1‑36 (MQ=25)
                         gtgATAAGCGGG.GGTGGCGTGATGCATTCCGTCTcc       >  2:1366940/1‑36 (MQ=37)
                          tgATAAGCGGG.GGTGGCGTGATGCATTCCGTCTCCt      <  2:1125671/36‑1 (MQ=37)
                          tgATAAGCGGG.GGTGGCGTGATGCATTCCGTCTCCt      >  1:2202101/1‑36 (MQ=37)
                            aTAAGCGGG.GGTGGCGTGATGCATTCCGTCTCCttt    >  1:359037/1‑36 (MQ=37)
                            aTAAGCGGG.GGTGGCGTGATGCATTCCGTCTCCttt    <  1:2099570/36‑1 (MQ=37)
                            aTAAGCGGG.GGTGGCGTGATGCATTCCGTCTCCttt    >  1:1800488/1‑36 (MQ=37)
                             tAAGCGGG.GGTGGCGTGATGCATTCCGTCTCCtttt   <  2:547277/36‑1 (MQ=37)
                              aaGC.GGAGGTGGCGTGATGCATTCCGTCTCCTTTTc  <  1:1319735/36‑1 (MQ=37)
                                  |                                
GAGCACAAGGCCGCCAACAATGGTGGTGATAAGC.GG.GGTGGCGTGATGCATTCCGTCTCCTTTTC  >  REL606/475255‑475319

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 18 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG