Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,758,6260A→C32.6% 15.043902yicFNAD‑dependent DNA ligase LigB
Reads supporting (aligned to +/- strand):  new base (0/14):  ref base (25/4):  total (25/18)

CGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCG  >  REL606/3758591‑3758660
                                   |                                  
catcAACCTGAGCTACAAGTTGATATTTCCAGGcca                                    <  1:1977441/33‑1 (MQ=37)
cGGCAACCTGAGCTACAGGTTGATATTTCCAGGcca                                    <  2:427857/36‑1 (MQ=255)
cGGCAACCCGAGCTACAGGATGATATTTACAGGcca                                    >  1:200010/1‑36 (MQ=25)
 ggCAACCTGAGCTACAGGTTGATATTTCCAGGccac                                   >  1:3722375/1‑36 (MQ=255)
 ggCAACCTGAGCTACAGGTTGATATTTCCAGGccac                                   >  1:2524776/1‑36 (MQ=255)
  gCAACCTGAGCTACAGGTTGATATTTCCAGGccacc                                  >  1:3134374/1‑36 (MQ=255)
  gCAACCTGAGCTACAGGTTGATATTTCCAGGccacc                                  >  1:131541/1‑36 (MQ=255)
  gCAACCGGAGCTACAGGTTGATATTTCCAGGccacc                                  <  1:1758311/36‑1 (MQ=255)
   cAACCTGAGCTACAGGTTGATATTTCCAGGccacca                                 <  1:3737375/36‑1 (MQ=255)
   cAACCTGAGCTACAGGTTGATATTTCCAGGccacca                                 <  1:1991353/36‑1 (MQ=255)
   cAACCTGAGCTACAGGTTGATATTTCCAGGccacca                                 <  2:2729625/36‑1 (MQ=255)
    aaCCTGAGCTACAGGTTGATATTTCCAGGCCACCAg                                <  2:74996/36‑1 (MQ=255)
      ccTGAGCTACAGGTTGATATTTCCAGGCCACCAGcc                              <  2:2735931/36‑1 (MQ=255)
      ccTGAGCTACAGGTTGATATTTCCAGGCCACCAGcc                              >  2:2515843/1‑36 (MQ=255)
      ccTGAGCTACAGGTTGATATTTCCAGGCCACCAGcc                              >  2:247329/1‑36 (MQ=255)
      ccTGAGCTACAGGTTGATATTTCCAGGCCACCAGcc                              >  2:2841939/1‑36 (MQ=255)
       cTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCa                             >  1:809703/1‑36 (MQ=255)
       cTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCa                             >  2:3234920/1‑36 (MQ=255)
        tGAGCTACAGGTTGATATTTCCAGGCCACCAGCCAc                            >  2:1643411/1‑36 (MQ=255)
         gAGCTAAGGTTTGATATTTCCGGCCCCCCAGCCAct                           <  1:1797335/36‑1 (MQ=9)
          aGCTACAGGTTGATATTTCCAGGCCACCAGCCActc                          >  1:775849/1‑36 (MQ=255)
          aGCTACAGGTTGATATTCCCAGGCCACCAGCCActc                          <  2:3686505/36‑1 (MQ=255)
            cTACAGGTTTATTTTTCCAGGCCCCCAGCCACTCTg                        <  2:1533315/36‑1 (MQ=18)
            cTACAGGTTGATATTTCCAGGCCACCAGCCActctt                        >  1:1116231/1‑35 (MQ=255)
              ataGGTTGATATTTCCGGCCCCCCAGCCACTCTGcc                      <  1:1658204/34‑1 (MQ=18)
               cAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCt                     >  2:2878975/1‑36 (MQ=255)
                aGGTTGATATTCCCACCCCCCCAGCCACTCTGCCTg                    <  2:1338195/36‑1 (MQ=19)
                 gtttGTTATTCCCAGCCCCCCAGCCACTCTGCCTgg                   <  1:378378/34‑1 (MQ=16)
                 ggTTGATATTTCCAGGCCACCAGCCACTCTGCCTgg                   >  1:2880500/1‑36 (MQ=255)
                    tGATATTTCCAGGCCACCAGCCACTCTGCCTGGccc                >  1:2530906/1‑36 (MQ=255)
                    tGATATTTCCAGGCCACCAGCCACTCTGCCTGGccc                >  2:2031202/1‑36 (MQ=255)
                    tGAGATTGCCAGCCCCCCAGCCACTCTGCCTGGccc                <  2:1177337/36‑1 (MQ=16)
                     gATTTTTCCAGGCCCCCAGCCACTCTGCCTGGCCCg               <  1:352856/36‑1 (MQ=38)
                     gATATTTCGAGGCCACCGGCCACTCTGCCTGGCCCg               <  1:586813/36‑1 (MQ=37)
                     gATATCCCCAGCCCCCCAGCCACTCTGCCTGGCCCg               <  2:1504945/36‑1 (MQ=16)
                     gATATCCCCAGCCCCCCAGCCACTCTGCCTGGCCCg               <  2:2654616/36‑1 (MQ=16)
                      atatTTCCAGGCCACCAGCCACTCTGCCTGGCCCgt              >  2:1136969/1‑35 (MQ=255)
                      atatTTCCAGGCCACCAGCCACTCTGCCTGGCCCgg              >  2:3037713/1‑36 (MQ=255)
                       tatTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGt             >  1:1647451/1‑36 (MQ=255)
                       tatTCCCGGGCCCCCAGCCACTCTGCCTGGCCCGGt             <  1:206875/36‑1 (MQ=37)
                        ttttCGAGGCCCCCAGCCACTCTGCCTGGCCCGGTa            <  2:1482180/35‑1 (MQ=37)
                        atTCCCAGGCCCCCAGCCACTCTGCCTGGCCCGGTa            <  1:3501111/36‑1 (MQ=38)
                         tttCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAg           >  2:2412235/1‑36 (MQ=255)
                         ttcccAGCCCCCCAGCCACTCTGCCTGGCCCGGTAg           <  2:2292952/33‑1 (MQ=25)
                          ttCCAGGCCACCAGCCACTCTGCTTGGCCCGGTAGc          >  2:1721887/1‑36 (MQ=255)
                          ttCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGc          >  1:45917/1‑36 (MQ=255)
                           tCCAGGCCACCAGCCACTCTACCCGGCCCGGTAGcc         >  1:691891/1‑36 (MQ=37)
                            ccAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCa        >  2:2650981/1‑36 (MQ=255)
                            ccAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCa        >  1:2770278/1‑36 (MQ=255)
                            ccAGCCCCCCAGCCACTCTGCCTGGCCCGGTAGCCa        <  1:2670839/36‑1 (MQ=37)
                             cAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCaa       >  1:2201472/1‑36 (MQ=255)
                              aGGCCACCAGCCACTCTGCCTGGCCCGGTAGGCatt      >  1:2527896/1‑34 (MQ=37)
                              aGGCCACCAGCCACTCTGACTGGCCCGGTAGCCAAt      >  1:959029/1‑36 (MQ=255)
                               ggCCACCCGCCACTCTGCCTGGCCCGGTACCCAATg     >  2:654717/1‑36 (MQ=38)
                                  caccaGCCACTCTGCCTGGCCCGGTAGCCAATGGCg  >  1:3584297/1‑36 (MQ=255)
                                  caccaGCCACTCTGCCTGGCCCGGTAGCCAATGGCg  <  2:1032974/36‑1 (MQ=255)
                                   |                                  
CGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCG  >  REL606/3758591‑3758660

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 13 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG