Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606260,1810T→G27.0% 10.863P9P (CCA→CCCpepDaminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  new base (17/0):  ref base (9/37):  total (26/37)

GAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTC  >  REL606/260146‑260215
                                   |                                  
gAACAGATTTTGGCAAAAATATCCCACAGCGGCtgt                                    <  1:4769882/36‑1 (MQ=255)
gAACAGATTTCGGCAAAAATATCCCACAGCGGCtgt                                    <  1:1310749/36‑1 (MQ=255)
 aaCAGATTTTGGCAAAAATACCCCACAGCGGCtgtg                                   >  1:5015018/1‑36 (MQ=255)
  aCAGATTTTGGCAAAAATATCCCACAGCGGCTGTgg                                  <  1:4046868/36‑1 (MQ=255)
  aCAGATTTTGGCAAAAATATCCCACAGCGGCTGTgg                                  <  1:2631259/36‑1 (MQ=255)
  aCAGATTTTGGCAAAAATATCCCACAGCGGCTGGgg                                  >  1:7380747/1‑36 (MQ=255)
  aCAGATTTTGGCAAAAATATCCCACAGCGGCTGGgg                                  >  1:4125106/1‑36 (MQ=255)
    aGATTTTGGCAAAAATATCCCACAGCGGCTGGgggg                                >  1:4325686/1‑34 (MQ=25)
    aGATTTTGGCAAAAATATCCCACAGCGGCTGGgggg                                >  1:2457585/1‑34 (MQ=25)
    aGATTTTGGCAAAAATATCCCACAGCGGCTGGgggg                                >  1:2939681/1‑34 (MQ=25)
     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:3883205/36‑1 (MQ=255)
           ggCAAAAATATCCCACAGCGGCTGTGGAGATAATTg                         <  1:6474990/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:4592274/36‑1 (MQ=255)
                          caGCGGCTGTGGAGATAATTGAGACAGTTCAGacac          <  1:1914292/36‑1 (MQ=255)
                          caGCGGCTGTGGAGATAATTGAGACAGTTCAGacac          <  1:2246707/36‑1 (MQ=255)
                          caGCGGCTGTGGAGATAATTGAGACAGTTCAGacac          <  1:5684251/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:6916980/36‑1 (MQ=255)
                             cGGCTGTGGAGATAATTGAGACAGTTCAGACACGtt       <  1:1065273/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)
                                   |                                  
GAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTC  >  REL606/260146‑260215

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