Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606987,4560G→T24.5% 12.249976ycaQhypothetical protein
Reads supporting (aligned to +/- strand):  new base (12/0):  ref base (16/21):  total (28/21)

ATACCCCAGCGCCGAAACGCCAGTATGGCT.ATTTTGTTCTGCCGTTATTACATCGTGGGCAATTAGTTGG  >  REL606/987421‑987490
                                    |                                  
atcccccAGCGCCGAAACGCCAGTATGGCT.ATTTTg                                    <  2:2360567/33‑1 (MQ=255)
ataCCCCAGCGCCGAAACGCCAGTATGGCT.ATTTTg                                    <  1:341919/36‑1 (MQ=255)
  aCCCCAGCGCCGAAACGCCAGTATGGCT.ATTTTGtt                                  <  2:1058048/36‑1 (MQ=255)
  aCCCCAGCGCCGAAACGCCAGTATGGCT.AATTTTtt                                  >  2:1450717/1‑36 (MQ=37)
   ccccAGCGCCGAAACGCCAGTATGGCT.ATTTTGTTc                                 <  2:1858903/36‑1 (MQ=255)
    cccAGCGCCGAAACGCCAGTATGGCT.ATTTTGTTCt                                <  1:2365522/36‑1 (MQ=255)
     ccAGCGCCGAAACGCCAGTATGGCT.ATTTTTTTCTg                               >  1:2644115/1‑36 (MQ=255)
     ccAGCGCCGAAACGCCAGTATGGCT.ATTTTGTTCTg                               <  1:589403/36‑1 (MQ=255)
       aGCGCCGAAGCGCCAGTATGGCT.ATTTTGTTCTGcc                             <  1:2067114/36‑1 (MQ=255)
       aGCGCCGAAACGCCAGTATGGCT.ATTTTTTTCTGcc                             >  2:1399570/1‑36 (MQ=255)
        gcgtCGAAGCGCCAGTATGGCT.ATTTTGTTCTGCCg                            <  1:2991849/36‑1 (MQ=37)
         cgcCGAAACGCCAGTATGGCT.AATTTTTTCTTCCGt                           >  1:3389590/1‑36 (MQ=25)
          gcCGAAACGCCAGTATGGCT.ATTTTGTTCTGCCGtt                          <  1:1171610/36‑1 (MQ=255)
           ccGAAACGCCAGTATGGCT.AATTTTTTTTTCCGttt                         >  1:3254397/1‑35 (MQ=16)
            cGAAACGCCAGTATGGCT.ATTTTGTTTTTTCGttat                        >  2:1770132/1‑36 (MQ=25)
            cGAAACGCCAGTATGGCT.ATTTTGTTCTGCCGttat                        >  2:1764930/1‑36 (MQ=255)
             gAAACGCCAGTATGGCT.ATTTTTTTTTTTCGttatt                       >  2:1665813/1‑36 (MQ=16)
             aaaaCCCCAGTATGGCT.ATTTTGTTCTGCCGttatt                       <  1:2959253/35‑1 (MQ=255)
              aaaCGCCAGTATGGCT.ATTTTTTTCTTCCGttattt                      >  1:3642873/1‑35 (MQ=25)
              aaaCGCCAGTATGGCT.AATTTGTTTTTCCGgtttta                      >  2:378895/1‑36 (MQ=12)
              aaaCGCCAGTATGGCT.AATTTGTTCTTCCCtcatta                      >  1:2592732/1‑36 (MQ=16)
               aaCGTCAATAGGGCT.ATTTTGTTCTGCCGTAATTAc                     <  2:2824792/36‑1 (MQ=16)
                aCGCCAGTATGGCT.ATTTTTTTTTGCCGttttttta                    >  1:2889758/1‑33 (MQ=16)
                aCGCCAGTATGGCT.AATTTTTTTTTCCGTTTTTACa                    >  1:3571605/1‑36 (MQ=16)
                 cGCCAGTATGGCT.ATTTTGTTCTGCCGTTATTACAt                   <  2:1074593/36‑1 (MQ=255)
                 cGCCAGTATGGCT.ATTTTGTTCTGCCGTTATTACAt                   <  1:3693947/36‑1 (MQ=255)
                 cGCCAGTATGCCT.ATTTTGTTCTGCCGTTATTACAt                   <  1:1105615/36‑1 (MQ=255)
                  cccAGTATGGCT.ATTTTGTTCTGCCGTTATTACATc                  <  2:3344130/35‑1 (MQ=255)
                     aGTATGGCTAATTTTTTTCT.CCGTTTTTATATcctt               >  1:3602634/1‑33 (MQ=0)
                      gTATGGCT.ATTTTTTTCTGCCGTTATTACATCGTgg              >  1:603292/1‑36 (MQ=255)
                      gTATGGCT.ATTTTGTTCTGCCGTTATTACATCGTgg              <  1:2666020/36‑1 (MQ=255)
                      gTATGGCT.ATTTTGTTCTGCCGTTATTACAACGTgg              >  2:174399/1‑36 (MQ=255)
                       tATGGCT.ATATTTTTCTTCCGTTATTACAACGttgt             >  2:3434015/1‑33 (MQ=12)
                          ggCT.ATTTTGTTCTGCCGTTATTACATCGTGGGCaa          >  1:254266/1‑36 (MQ=255)
                          ggCG.ATTTTGTTCTGCCGTTATTACATCGTGGGCaa          <  2:1900621/36‑1 (MQ=255)
                           gCT.ATCTTGTTCTGCCGCTATCACATCGTGGGCAAt         >  1:131678/1‑36 (MQ=25)
                            cT.ATTTTGTTCTGCCGTTATTACATCGTGGGCAAtt        <  1:2045970/36‑1 (MQ=255)
                             t.ATTTTGTTCTGCCGTTATTACATCGTGGGCAATTa       <  1:3687967/36‑1 (MQ=255)
                             t.ATTTTGTTCTGCCGTTATTACATCGTGGGCAATTa       <  2:3338461/36‑1 (MQ=255)
                             t.ATTTTGTTCTGCCGTTATTACATCGTGGGAAATTa       >  2:95605/1‑36 (MQ=255)
                               aTTTTGTTCTGCCGTTATTACATCGTGGGCAATTAg      <  1:251732/36‑1 (MQ=255)
                               aTTTTGTTCTGCCGTTATAACATCGTGGGCAATTAg      >  1:1697611/1‑36 (MQ=255)
                                 tttGTTCTGCCGTTATTACATCGTGGGCAATTTGtt    >  1:2322703/1‑36 (MQ=255)
                                 tttGTTCTGCCGTTATTACATCGTGGGCAATTAGtt    >  2:2718392/1‑36 (MQ=255)
                                  ttGTTCTGCCGTTATTACATCGTGGGCAATTAGtcg   >  1:770751/1‑34 (MQ=255)
                                  ttGTTCTGCCGTTATTACATCGTGGGCAATTAGTTg   >  2:2350705/1‑36 (MQ=255)
                                  ttGTTCTGCCGTTATTACATCGTGGGCAATTAGTTg   >  1:1694569/1‑36 (MQ=255)
                                  ttGTTCTGCCGTTATTACATCGTGGGCAATTAGTTg   <  2:2895822/36‑1 (MQ=255)
                                  ttGTTCTGCCGTTATTACATCGTGGGCAATTAGTTg   <  2:1702718/36‑1 (MQ=255)
                                  ttGTTCTGCCGTTATTACATCGTGGGCAATTAATTg   >  2:2096475/1‑36 (MQ=255)
                                   tGTTCTGCCGTTATTACATCGTGGGCAATTAGTTgg  >  1:2914539/1‑36 (MQ=255)
                                    |                                  
ATACCCCAGCGCCGAAACGCCAGTATGGCT.ATTTTGTTCTGCCGTTATTACATCGTGGGCAATTAGTTGG  >  REL606/987421‑987490

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 13 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG