Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,610,8600A→G17.9% 11.378D191D (GAT→GACyhjJpredicted zinc‑dependent peptidase
Reads supporting (aligned to +/- strand):  new base (14/0):  ref base (26/38):  total (40/38)

CCGCTTCAACGGGTTGTTTCAGCGGATCGGCAGGATCGTGACCTAACAAGGTTGACCCTTTCAGACGATA  >  REL606/3610826‑3610895
                                  |                                   
ccGCTTCAACGGGTTGTTTCAGCGGATCGGCAGGAt                                    <  1:3278244/36‑1 (MQ=255)
ccGCTTCAACGGGTTGTTTCAGCGGATCAGCAGGAt                                    >  1:5805954/1‑36 (MQ=255)
 cGCTTCAACGGGTTGTTTCAGCGGATCGGCAGGATc                                   <  1:4642592/36‑1 (MQ=255)
 cGCTTCAACGGGGTGTGTCAGCGGATCGGCAGGATc                                   <  1:313385/36‑1 (MQ=25)
  gCTTCAACGGGTTGTTTCAGCGGATCGGCCCGATCg                                  >  1:4631585/1‑36 (MQ=25)
  gCTTCAACGGGTTGTTTCAGCGGATCGGCAGGGTCg                                  >  1:5959988/1‑36 (MQ=255)
  gCTTCAACGGGTTGTTTCAGCGGATCGGCAGGGTCg                                  >  1:6423120/1‑36 (MQ=255)
   cTTCAACGGGTTGTTTCAGCGGATCGGCAGGGTCGt                                 >  1:5607048/1‑36 (MQ=255)
   cTTCAACGGGTTGTTTCAGCGGATCGGCAGGATTGt                                 >  1:2877023/1‑36 (MQ=255)
    ttCAACGTGTTGTTTCAGCGGATCGGCAGGATCGTg                                <  1:183028/36‑1 (MQ=255)
    ttCAACGGGTTGTTTCAGCGGATCGGCAGGGTCGTg                                >  1:5697738/1‑36 (MQ=255)
      cAACGGGTTGTTTCAGCGGATCGGCCGGGTTGTgtc                              >  1:4701607/1‑34 (MQ=18)
      cAACGGGTTGTTTCAGCGGATCGGCAGGAACGTGAc                              >  1:5363089/1‑36 (MQ=255)
       aaCGGGTTGTTTCAGCGGATCGGCAGGGTCGTTAcc                             >  1:477089/1‑36 (MQ=21)
       aaCGGGTTGTTTCAGCGGATCGGCAGGGTCCTCAcc                             >  1:4871065/1‑36 (MQ=18)
       aaCGGGTTGTTTCAGCGGATCGGCAGGATCGTgggc                             >  1:7066687/1‑33 (MQ=255)
       aaCGGGTTGTTTCAGCGGATCGGCAGGATCGTgggc                             >  1:2772137/1‑33 (MQ=255)
       aaCGGGTTGTTTCAGCGGATCGGCAGGATCGTGAcc                             <  1:2220982/36‑1 (MQ=255)
       aaCGGGTTGTTTCAGCGGATCGGCAGGATCGTGAcc                             <  1:1750333/36‑1 (MQ=255)
       aaCGGGTTGTTTCAGCGGATCGGCAGGATCGTGAcc                             <  1:4827386/36‑1 (MQ=255)
        aCGGGTTGTTTCAGCGGATCGGCAGGATCGTGACCt                            <  1:3989716/36‑1 (MQ=255)
         cGGGTTGTTTCAGCGGATCGGCAGGATCGTGAGCTa                           >  1:4788427/1‑36 (MQ=255)
         cGGGTTGTTTCAGCGGATCGGCAGGATCGTGACCTa                           <  1:5031502/36‑1 (MQ=255)
         cGGGTTGTTTCAGCGGATCGGCAGGATCGTGACCTa                           <  1:414185/36‑1 (MQ=255)
         cGGGTGGTTTCAGCGGATCGGCAGGATCGTGACCTa                           <  1:6580737/36‑1 (MQ=255)
          gggTTGTTTCAGCGGATCGGCAGGGTGGTGTGCTaa                          >  1:7417649/1‑36 (MQ=17)
          gggTTGTTTCAGCGGATCGGCAGGGTCGCGACCTaa                          >  1:5158466/1‑36 (MQ=21)
          gggTTGTTTCAGCGGATCGGCAGGGACGTGACCTaa                          >  1:592718/1‑36 (MQ=25)
          gggTTGTTTCAGCGGATCGGCAGGATTGTGACCTaa                          >  1:7227431/1‑36 (MQ=255)
          gggTTGTTTCAGCGGATCGGCAGGATCGTGACCTaa                          <  1:4453970/36‑1 (MQ=255)
          gggTTGTTTCAGCGGATCGGCAGGATCGTGACCTaa                          <  1:1283292/36‑1 (MQ=255)
           tgttgttTCAGCGGATCGGCAGGATCGTGACCTAAc                         <  1:69454/35‑1 (MQ=255)
           ggTTGTTTCAGCGGATCGGCAGGATCGTGACCTAAc                         <  1:622473/36‑1 (MQ=255)
           ggTTGTTTCAGCGGATCGGCAGGATCGTGACCTAAc                         <  1:3998276/36‑1 (MQ=255)
            gttgttTCAGCGGATCGGCAGGGTCGGGCTCTAAcc                        >  1:5976376/1‑35 (MQ=17)
            gttgttTCAGCGGATCGGCAGGATCGTGACCTAACa                        <  1:5976032/36‑1 (MQ=255)
            gttgttTCAGCGGATCGGCAGGATCGTGACCTAACa                        <  1:1615568/36‑1 (MQ=255)
             ttgttTCAGCGGATCGGCAGGATCGTGACCTAACaa                       <  1:941153/36‑1 (MQ=255)
             ttgttTCAGCGGATCGGCAGGATCGTGACCCAACaa                       >  1:419704/1‑36 (MQ=255)
             ttgttTCAGCGGATCGGCAGGATCGTCACCTAACaa                       >  1:2174206/1‑36 (MQ=255)
              tgttTCAGCGGATCGGCAGGATCGTGACCTAACAAg                      <  1:3054936/36‑1 (MQ=255)
              tgttTCAGCGGATCGGCACGATCGTGACCCAAAAAg                      >  1:7317554/1‑36 (MQ=18)
               gttTCAGCGGATCGGCAGGGTCGTGGCCCAAAAAgg                     >  1:5237683/1‑36 (MQ=17)
                tttCAGCGGATCGGCAGGATCGTGACCTAACAAGGt                    <  1:2910402/36‑1 (MQ=255)
                tttCAGCGGATCGGCAGGATCGTGACCTAACAAGGt                    <  1:1032712/36‑1 (MQ=255)
                   cAGCGGATCGGCAGGATCCGGACCTAAAAACGATGa                 >  1:4814307/1‑36 (MQ=12)
                     gCGGATCGGCAGGGTCGTGGCCTAACAAGGTTGAcc               >  1:2469929/1‑36 (MQ=255)
                     gCGGATCGGCAGGATCGTGACCTAACAAGGTTGAcc               <  1:2151282/36‑1 (MQ=255)
                     gCGGATCGGCAGGATCGTGACCTAACAAGGTTGAcc               <  1:6165920/36‑1 (MQ=255)
                     gCGGATCGGCAGGATCGTGACCTAACAAGGTTGAcc               <  1:845172/36‑1 (MQ=255)
                      cGGATCCGCAGGATCGTGACCGAAAAAGGATGGccc              >  1:6040935/1‑36 (MQ=12)
                       ggATCGGCAGGGTCGTGACCTAACAAGGGTCAcccc             >  1:2514571/1‑35 (MQ=21)
                       ggATCGGCAGGATCGTGACCTAACAAGGTTGACCCt             <  1:2690176/36‑1 (MQ=255)
                        gATCGGCAGGATCGTGACCTAACAAGGTTGATCCtt            >  1:1947611/1‑36 (MQ=255)
                        gATCGGCAGGATCGTGACCTAACAAGGTTGACCCtt            <  1:2134286/36‑1 (MQ=255)
                        gATCGGCAGGATCGTGACCTAACAAGGTTGACCCtt            >  1:3859644/1‑36 (MQ=255)
                        gATCGGCAGGATCGTGACCTAACAAGGTTGACCCtt            >  1:4157926/1‑36 (MQ=255)
                        gATCGGCAGGATCGTGACCTAAAAAGGTTGAcccct            >  1:1410584/1‑34 (MQ=25)
                        gATAGGCAGGATCGTGACCTAACAAGGTTGACCCtt            <  1:6225217/36‑1 (MQ=255)
                         aTCGGCAGGATCGTGACCTAACAAGGTTGACCCttt           <  1:6374309/36‑1 (MQ=255)
                         aTCGGCAGGATCGTGACCTAACAAGGTTGACCCttt           <  1:5335273/36‑1 (MQ=255)
                          tCGGCAGGATCGTGACCTAACAAGGTTGACCCTTTc          <  1:2364835/36‑1 (MQ=255)
                           cGGCAGGATCGTGACCTAACAAGGTTGACCCTTTCa         >  1:3314894/1‑36 (MQ=255)
                           cGGCAGGATCGTGACCTAACAAGGTTGACCCTTTCa         <  1:1828177/36‑1 (MQ=255)
                            ggCAGGATCGTGACCTAACAAGGTTGACCCTTTCac        >  1:1154886/1‑35 (MQ=255)
                            ggCAGGATCGTGACCTAACAAGGTTGACCCTTTCAg        >  1:6357551/1‑36 (MQ=255)
                            ggCAGGATCGTGACCTAACAAGGTTGACCCTTTCAg        >  1:1008814/1‑36 (MQ=255)
                               aGGATCGTGACCTAACAAGGTTGACCCTTTCAGACg     >  1:2836643/1‑36 (MQ=255)
                               aGGATCGTGACCTAACAAGGTTGACCCTTTCAGACg     >  1:2871115/1‑36 (MQ=255)
                               aGGATCGTGACCTAACAAGGTTGACCCTTTCAGACg     >  1:6121049/1‑36 (MQ=255)
                               aGGATCGTGACCTAACAAGGTTGACCCTTTCAGACg     <  1:3030769/36‑1 (MQ=255)
                               aGGATCGTGACCTAACAAGGTTGACCCTTTCAGACg     <  1:6946060/36‑1 (MQ=255)
                               aGGATCGTGACCTAACAAGGTTGACCCTTTCAGACg     <  1:6971048/36‑1 (MQ=255)
                                ggATCGTGACCTAACAAGGTTGACCCTTTCAGACGa    >  1:3651842/1‑36 (MQ=255)
                                ggATCGTGACCTAACAAGGTTGACCCTTTCAGACGa    <  1:1867990/36‑1 (MQ=255)
                                 gATCGTGACCTAACAAGGTTGACCCTTTCAGACGAt   <  1:5355195/36‑1 (MQ=255)
                                 gATCGTGACCTAACAAGGTTGACCCTTTCAGACGAt   >  1:3764563/1‑36 (MQ=255)
                                 gATCGTGACCTAACAAGGTTGACCCTTTCAGACGAt   <  1:4251151/36‑1 (MQ=255)
                                  aTCGTGACCTAACAAGGTTGCCCCTTTCAGACGATa  >  1:4835802/1‑36 (MQ=255)
                                  aTCGTGACCTAACAAGGTTGACCCTTTCAGACGATa  <  1:2310191/36‑1 (MQ=255)
                                  |                                   
CCGCTTCAACGGGTTGTTTCAGCGGATCGGCAGGATCGTGACCTAACAAGGTTGACCCTTTCAGACGATA  >  REL606/3610826‑3610895

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 14 ≤  ATCG  < 28 ≤  ATCG  < 32 ≤  ATCG