Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606260,1840A→G31.0% 19.858S8S (TCT→TCCpepDaminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  new base (17/1):  ref base (6/34):  total (23/35)

GATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCT  >  REL606/260151‑260219
                                 |                                   
gATTTTGGCAAAAATATCCCACAGCGGCTGTGgagg                                   >  1:5953710/1‑35 (MQ=255)
 aTTTTGGCAAAAATATCCCACAGCGGCTGGGGGGAt                                  >  1:6253391/1‑36 (MQ=21)
  ttttGTCAAAAATATCCCACAGCGGCTGTGGAGATa                                 <  1:5558343/36‑1 (MQ=255)
  ttttGGCAAAAATATCCCACAGCGGCTGTGGAGATa                                 <  1:5963395/36‑1 (MQ=255)
   tttGGCAAAAATATCCCACAGCGGCTGTGGGGGTaa                                >  1:6076472/1‑36 (MQ=21)
   tttGGCAAAAATATCCCACAGCGGCTGTGGAGATaa                                <  1:6464646/36‑1 (MQ=255)
    ttGGAAAAAATATCCCACAGCGGCTGTGGAGATAAt                               <  1:5971333/36‑1 (MQ=255)
     tagCAAAAATATCCCACAGCGGCTGTGGAGATAAtt                              <  1:7132483/34‑1 (MQ=255)
     tGGCAAAAATATCCCACAGCGGCTGTGGGGATAAtt                              <  1:1928793/36‑1 (MQ=255)
     tGGCAAAAATATCCCACAGCGGCTGTGGAGATAAtt                              <  1:3165811/36‑1 (MQ=255)
      ggaaaaaaTATCCCACAGCGGCTGTGGAGATAATTg                             <  1:4551175/33‑1 (MQ=255)
      ggCAAAAATATCCCACAGCGGCTGTGGAGATAATTg                             <  1:6474990/36‑1 (MQ=255)
      ggCAAAAATATCCCACAGCGGCTGTGGAGATAATTg                             <  1:3883205/36‑1 (MQ=255)
      ggCAAAAATATCCCACAGCGGCTGGGGGGGTAATTg                             >  1:3507315/1‑36 (MQ=18)
       gCAAAAATATCCCACAGCGGCTGTGGGGATAAttta                            >  1:3412167/1‑34 (MQ=25)
       gCAAAAATATCCCACAGCGGCTGTGGAGATAATTga                            <  1:4154005/36‑1 (MQ=255)
       gCAAAAATATCCCACAGCGGCTGGGGGGATAATTga                            >  1:4174379/1‑36 (MQ=21)
        ggaaaaTATCGCTCAGCGGCTGTGGAGATAATTgag                           <  1:400788/34‑1 (MQ=21)
        cAAAAATATCCCACAGCGGCTGTGGAGATAATTgag                           <  1:1243762/36‑1 (MQ=255)
        cAAAAATATCCCACAGCGGCTGGGGGGATAATTgag                           >  1:5332539/1‑36 (MQ=21)
        cAAAAACATCCCACAGCGGCTGGGGGGGTACTTgag                           >  1:5040492/1‑36 (MQ=12)
         aaaaaTATCCCACAGCGGCTGGGGGGATAATAgaga                          >  1:903730/1‑36 (MQ=18)
         aaaaaTATCCCACAGCGGCTGGGGGGATAATAgaaa                          >  1:1281853/1‑34 (MQ=18)
         aaaaaTAGCAAACAGCGTCTGTGGAGATAATTgaga                          <  1:2106407/36‑1 (MQ=18)
           taaTATCCCACAGCGGCGGTGGAGATAATTGAGACa                        <  1:4502418/35‑1 (MQ=25)
           aaaTATCCCACAGCGGCTGTGGAGATAATTGAGACa                        <  1:2514378/36‑1 (MQ=255)
           aaaTATCCCACAGCGGCTGGGGGGATAATTtagaaa                        >  1:1215212/1‑34 (MQ=18)
           aaaTATCCCACAGCGGCTGGGGGGATAATTgagaaa                        >  1:7205115/1‑34 (MQ=21)
            aaTATCCCACAGCGGCTGGGGGGATAATTAAGACaa                       >  1:7057241/1‑35 (MQ=21)
            aaTATCCCACAGCGGCTGGGGGGACAATTAAGACAg                       >  1:104454/1‑36 (MQ=18)
             atatCCCACAGCGGCTGTGGAGATAATTGAGACAGt                      <  1:3505996/36‑1 (MQ=255)
             aaatCCCACAGCGGCTGTGGAGATAATTGAGACAGt                      <  1:1222013/34‑1 (MQ=255)
              tatCCCACAGCGGCTGTGGGGGTAATTGAGACAGtt                     >  1:5668707/1‑36 (MQ=25)
              tatCCCACAGCGGCTGTGGAGATAATTGAGACAGtt                     <  1:5615258/36‑1 (MQ=255)
              aatCCCACAGCGGCTGTGGAGATAATTGAGACAGtt                     <  1:3337758/35‑1 (MQ=255)
               atCCCACAGCGGCTGTGGGGGTAATTGAGACAGttt                    >  1:4600470/1‑35 (MQ=25)
               atCCCACAGCGGCTGTGGAGAAAAATTAGACAAtta                    >  1:3578273/1‑35 (MQ=18)
                tCCCACAGCGGCTGTGGGGATAATTGAGACAGtttt                   >  1:965025/1‑34 (MQ=25)
                 cccACAGCGGCTGTGGAGATAATTGAGACAGTTCAg                  <  1:3698297/36‑1 (MQ=255)
                 cccACAGCGGCTGGGGGGGTAATTGAGGCAGttttg                  >  1:5697099/1‑33 (MQ=12)
                   cacaGCGGCTGTGGAGATAATTGAGACAGTTCAGac                <  1:7256308/36‑1 (MQ=255)
                   cacaGCGGCTGTGGAGACAATTGAGACGGTTCAGac                <  1:4635517/36‑1 (MQ=38)
                     caGCGGCTGTGGAGATAATTGAGACAGTTCAGacac              <  1:1914292/36‑1 (MQ=255)
                     caGCGGCTGTGGAGATAATTGAGACAGTTCAGacac              <  1:5684251/36‑1 (MQ=255)
                     caGCGGCTGTGGAGATAATTGAGACAGTTCAGacac              <  1:2246707/36‑1 (MQ=255)
                     caGCGGCTGTGGAGATAATTGAGACAGTTCAGacac              <  1:4592274/36‑1 (MQ=255)
                      aGCGGCTGTGGAGATAATTGAGACAGTTCAGACACg             <  1:342305/36‑1 (MQ=255)
                        cGGGTGTGGAGATAATTGAGACAGTTCAGACACGtt           <  1:3907701/36‑1 (MQ=255)
                        cGGCTGTGGAGATAATTTAGACAGTTCATACACgct           >  1:2030647/1‑34 (MQ=255)
                        cGGCTGTGGAGATAATTGAGACAGTTCAGACACGtt           <  1:1065273/36‑1 (MQ=255)
                        cGGCTGTGGAGATAATTGAGACAGTTCAGACACGtt           <  1:6916980/36‑1 (MQ=255)
                         ggCTGTGGAGATAATTGAGACAGTTCAGACACGTTa          >  1:3802687/1‑36 (MQ=255)
                          gCTGTGGAGATAATTGAGACAGTTCAGACACGTTaa         <  1:6741946/36‑1 (MQ=255)
                           cTGTGGAGATAATTGAGACAGTTCAGACACGTTAAg        <  1:2641446/36‑1 (MQ=255)
                           cTGTGGAGATAATTGAGACAGTTCAGACACGTTAAg        <  1:5817152/36‑1 (MQ=255)
                           cTGTGGAGATAATTGAGACAGTTCAGACACGTTAAg        <  1:5560059/36‑1 (MQ=255)
                            tgtgGAGATAATTGAGACAGTTCAGACACGTTAAGt       >  1:4824255/1‑36 (MQ=255)
                             gtgGAGATAATTGAGACAGTTCAGACACGTTAAGtc      >  1:3237328/1‑36 (MQ=255)
                               ggAGATAATTGAGACAGTTCAGACACGTTAAGtctc    >  1:4756200/1‑36 (MQ=255)
                                gagaTAATTGAGACAGTTCAGACACGTTAAGTCTcc   >  1:5246442/1‑36 (MQ=255)
                                 agaTAATTGAGACAGTTCAGACACGTTAAGTCTCCt  <  1:7063244/36‑1 (MQ=255)
                                 |                                   
GATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCT  >  REL606/260151‑260219

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