Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606260,1840A→G32.2% 22.15924pepDaminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  new base (18/1):  ref base (6/34):  total (24/35)

GATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCT  >  REL606/260151‑260219
                                 |                                   
gATTTTGGCAAAAATATCCCACAGCGGCTGTGgagg                                   >  2:2214500/1‑35 (MQ=255)
 aTTTTGGCAAAAATATCCCACAGCGGCTGGGGGGAt                                  >  2:2514181/1‑36 (MQ=37)
  ttttGTCAAAAATATCCCACAGCGGCTGTGGAGATa                                 <  2:1819133/36‑1 (MQ=255)
  ttttGGCAAAAATATCCCACAGCGGCTGTGGAGATa                                 <  2:2224185/36‑1 (MQ=255)
   tttGGCAAAAATATCCCACAGCGGCTGTGGGGGTaa                                >  2:2337262/1‑36 (MQ=37)
   tttGGCAAAAATATCCCACAGCGGCTGTGGAGATaa                                <  2:2725436/36‑1 (MQ=255)
    ttGGAAAAAATATCCCACAGCGGCTGTGGAGATAAt                               <  2:2232123/36‑1 (MQ=255)
     tagCAAAAATATCCCACAGCGGCTGTGGAGATAAtt                              <  2:3393273/34‑1 (MQ=255)
     tGGCAAAAATATCCCACAGCGGCTGTGGGGATAAtt                              <  1:1928793/36‑1 (MQ=255)
     tGGCAAAAATATCCCACAGCGGCTGTGGAGATAAtt                              <  1:3165811/36‑1 (MQ=255)
      ggaaaaaaTATCCCACAGCGGCTGTGGAGATAATTg                             <  2:811965/33‑1 (MQ=255)
      ggCAAAAATATCCCACAGCGGCTGTGGAGATAATTg                             <  2:143995/36‑1 (MQ=255)
      ggCAAAAATATCCCACAGCGGCTGTGGAGATAATTg                             <  2:2735780/36‑1 (MQ=255)
      ggCAAAAATATCCCACAGCGGCTGGGGGGGTAATTg                             >  1:3507315/1‑36 (MQ=21)
       gCAAAAATATCCCACAGCGGCTGTGGGGATAAttta                            >  1:3412167/1‑34 (MQ=37)
       gCAAAAATATCCCACAGCGGCTGTGGAGATAATTga                            <  2:414795/36‑1 (MQ=255)
       gCAAAAATATCCCACAGCGGCTGGGGGGATAATTga                            >  2:435169/1‑36 (MQ=37)
        ggaaaaTATCGCTCAGCGGCTGTGGAGATAATTgag                           <  1:400788/34‑1 (MQ=37)
        cAAAAATATCCCACAGCGGCTGTGGAGATAATTgag                           <  1:1243762/36‑1 (MQ=255)
        cAAAAATATCCCACAGCGGCTGGGGGGATAATTgag                           >  2:1593329/1‑36 (MQ=37)
        cAAAAACATCCCACAGCGGCTGGGGGGGTACTTgag                           >  2:1301282/1‑36 (MQ=16)
         aaaaaTATCCCACAGCGGCTGGGGGGATAATAgaga                          >  1:903730/1‑36 (MQ=21)
         aaaaaTATCCCACAGCGGCTGGGGGGATAATAgaaa                          >  1:1281853/1‑34 (MQ=18)
         aaaaaTAGCAAACAGCGTCTGTGGAGATAATTgaga                          <  1:2106407/36‑1 (MQ=19)
           taaTATCCCACAGCGGCGGTGGAGATAATTGAGACa                        <  2:763208/35‑1 (MQ=255)
           aaaTATCCCACAGCGGCTGTGGAGATAATTGAGACa                        <  1:2514378/36‑1 (MQ=255)
           aaaTATCCCACAGCGGCTGGGGGGATAATTtagaaa                        >  1:1215212/1‑34 (MQ=21)
           aaaTATCCCACAGCGGCTGGGGGGATAATTgagaaa                        >  2:3465905/1‑34 (MQ=25)
            aaTATCCCACAGCGGCTGGGGGGATAATTAAGACaa                       >  2:3318031/1‑35 (MQ=25)
            aaTATCCCACAGCGGCTGGGGGGACAATTAAGACAg                       >  1:104454/1‑36 (MQ=19)
             atatCCCACAGCGGCTGTGGAGATAATTGAGACAGt                      <  1:3505996/36‑1 (MQ=255)
             aaatCCCACAGCGGCTGTGGAGATAATTGAGACAGt                      <  1:1222013/34‑1 (MQ=255)
              tatCCCACAGCGGCTGTGGGGGTAATTGAGACAGtt                     >  2:1929497/1‑36 (MQ=38)
              tatCCCACAGCGGCTGTGGAGATAATTGAGACAGtt                     <  2:1876048/36‑1 (MQ=255)
              aatCCCACAGCGGCTGTGGAGATAATTGAGACAGtt                     <  1:3337758/35‑1 (MQ=255)
               atCCCACAGCGGCTGTGGGGGTAATTGAGACAGttt                    >  2:861260/1‑35 (MQ=37)
               atCCCACAGCGGCTGTGGAGAAAAATTAGACAAtta                    >  1:3578273/1‑35 (MQ=16)
                tCCCACAGCGGCTGTGGGGATAATTGAGACAGtttt                   >  1:965025/1‑34 (MQ=38)
                 cccACAGCGGCTGTGGAGATAATTGAGACAGTTCAg                  <  1:3698297/36‑1 (MQ=255)
                 cccACAGCGGCTGGGGGGGTAATTGAGGCAGttttg                  >  2:1957889/1‑33 (MQ=16)
                   cacaGCGGCTGTGGAGATAATTGAGACAGTTCAGac                <  2:3517098/36‑1 (MQ=255)
                   cacaGCGGCTGTGGAGACAATTGAGACGGTTCAGac                <  2:896307/36‑1 (MQ=255)
                    acaGCGGCTGGGGGGATAATTGAGAAAGGTCACacc               >  2:2159737/1‑35 (MQ=16)
                     caGCGGCTGTGGAGATAATTGAGACAGTTCAGacac              <  2:1945041/36‑1 (MQ=255)
                     caGCGGCTGTGGAGATAATTGAGACAGTTCAGacac              <  1:1914292/36‑1 (MQ=255)
                     caGCGGCTGTGGAGATAATTGAGACAGTTCAGacac              <  2:853064/36‑1 (MQ=255)
                     caGCGGCTGTGGAGATAATTGAGACAGTTCAGacac              <  1:2246707/36‑1 (MQ=255)
                      aGCGGCTGTGGAGATAATTGAGACAGTTCAGACACg             <  1:342305/36‑1 (MQ=255)
                        cGGGTGTGGAGATAATTGAGACAGTTCAGACACGtt           <  2:168491/36‑1 (MQ=255)
                        cGGCTGTGGAGATAATTTAGACAGTTCATACACgct           >  1:2030647/1‑34 (MQ=37)
                        cGGCTGTGGAGATAATTGAGACAGTTCAGACACGtt           <  1:1065273/36‑1 (MQ=255)
                        cGGCTGTGGAGATAATTGAGACAGTTCAGACACGtt           <  2:3177770/36‑1 (MQ=255)
                         ggCTGTGGAGATAATTGAGACAGTTCAGACACGTTa          >  2:63477/1‑36 (MQ=255)
                          gCTGTGGAGATAATTGAGACAGTTCAGACACGTTaa         <  2:3002736/36‑1 (MQ=255)
                           cTGTGGAGATAATTGAGACAGTTCAGACACGTTAAg        <  1:2641446/36‑1 (MQ=255)
                           cTGTGGAGATAATTGAGACAGTTCAGACACGTTAAg        <  2:1820849/36‑1 (MQ=255)
                           cTGTGGAGATAATTGAGACAGTTCAGACACGTTAAg        <  2:2077942/36‑1 (MQ=255)
                            tgtgGAGATAATTGAGACAGTTCAGACACGTTAAGt       >  2:1085045/1‑36 (MQ=255)
                             gtgGAGATAATTGAGACAGTTCAGACACGTTAAGtc      >  1:3237328/1‑36 (MQ=255)
                               ggAGATAATTGAGACAGTTCAGACACGTTAAGtctc    >  2:1016990/1‑36 (MQ=255)
                                gagaTAATTGAGACAGTTCAGACACGTTAAGTCTcc   >  2:1507232/1‑36 (MQ=255)
                                 agaTAATTGAGACAGTTCAGACACGTTAAGTCTCCt  <  2:3324034/36‑1 (MQ=255)
                                 |                                   
GATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCT  >  REL606/260151‑260219

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG