Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606260,1810T→G28.1% 12.06427pepDaminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  new base (18/0):  ref base (9/37):  total (27/37)

GAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTC  >  REL606/260146‑260215
                                   |                                  
gAACAGATTTTGGCAAAAATATCCCACAGCGGCtgt                                    <  2:1030672/36‑1 (MQ=255)
gAACAGATTTCGGCAAAAATATCCCACAGCGGCtgt                                    <  1:1310749/36‑1 (MQ=255)
 aaCAGATTTTGGCAAAAATACCCCACAGCGGCtgtg                                   >  2:1275808/1‑36 (MQ=255)
  aCAGATTTTGGCAAAAATATCCCACAGCGGCTGTgg                                  <  1:2631259/36‑1 (MQ=255)
  aCAGATTTTGGCAAAAATATCCCACAGCGGCTGTgg                                  <  2:307658/36‑1 (MQ=255)
  aCAGATTTTGGCAAAAATATCCCACAGCGGCTGGgg                                  >  2:3641537/1‑36 (MQ=255)
  aCAGATTTTGGCAAAAATATCCCACAGCGGCTGGgg                                  >  2:385896/1‑36 (MQ=255)
    aGATTTTGGCAAAAATATCCCACAGCGGCTGGgggg                                >  1:2939681/1‑34 (MQ=37)
    aGATTTTGGCAAAAATATCCCACAGCGGCTGGgggg                                >  2:586476/1‑34 (MQ=37)
    aGATTTTGGCAAAAATATCCCACAGCGGCTGGgggg                                >  1:2457585/1‑34 (MQ=37)
     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          <  2:853064/36‑1 (MQ=255)
                          caGCGGCTGTGGAGATAATTGAGACAGTTCAGacac          <  1:1914292/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)
                                   |                                  
GAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTC  >  REL606/260146‑260215

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