Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,610,8600A→G17.9% 10.478573yhjJpredicted 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                                    >  2:2066744/1‑36 (MQ=255)
 cGCTTCAACGGGTTGTTTCAGCGGATCGGCAGGATc                                   <  2:903382/36‑1 (MQ=255)
 cGCTTCAACGGGGTGTGTCAGCGGATCGGCAGGATc                                   <  1:313385/36‑1 (MQ=37)
  gCTTCAACGGGTTGTTTCAGCGGATCGGCCCGATCg                                  >  2:892375/1‑36 (MQ=37)
  gCTTCAACGGGTTGTTTCAGCGGATCGGCAGGGTCg                                  >  2:2220778/1‑36 (MQ=255)
  gCTTCAACGGGTTGTTTCAGCGGATCGGCAGGGTCg                                  >  2:2683910/1‑36 (MQ=255)
   cTTCAACGGGTTGTTTCAGCGGATCGGCAGGGTCGt                                 >  2:1867838/1‑36 (MQ=255)
   cTTCAACGGGTTGTTTCAGCGGATCGGCAGGATTGt                                 >  1:2877023/1‑36 (MQ=255)
    ttCAACGTGTTGTTTCAGCGGATCGGCAGGATCGTg                                <  1:183028/36‑1 (MQ=255)
    ttCAACGGGTTGTTTCAGCGGATCGGCAGGGTCGTg                                >  2:1958528/1‑36 (MQ=255)
      cAACGGGTTGTTTCAGCGGATCGGCCGGGTTGTgtc                              >  2:962397/1‑34 (MQ=21)
      cAACGGGTTGTTTCAGCGGATCGGCAGGAACGTGAc                              >  2:1623879/1‑36 (MQ=255)
       aaCGGGTTGTTTCAGCGGATCGGCAGGGTCGTTAcc                             >  1:477089/1‑36 (MQ=25)
       aaCGGGTTGTTTCAGCGGATCGGCAGGGTCCTCAcc                             >  2:1131855/1‑36 (MQ=21)
       aaCGGGTTGTTTCAGCGGATCGGCAGGATCGTgggc                             >  2:3327477/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                             <  2:1088176/36‑1 (MQ=255)
        aCGGGTTGTTTCAGCGGATCGGCAGGATCGTGACCt                            <  2:250506/36‑1 (MQ=255)
         cGGGTTGTTTCAGCGGATCGGCAGGATCGTGAGCTa                           >  2:1049217/1‑36 (MQ=255)
         cGGGTTGTTTCAGCGGATCGGCAGGATCGTGACCTa                           <  1:414185/36‑1 (MQ=255)
         cGGGTTGTTTCAGCGGATCGGCAGGATCGTGACCTa                           <  2:1292292/36‑1 (MQ=255)
         cGGGTGGTTTCAGCGGATCGGCAGGATCGTGACCTa                           <  2:2841527/36‑1 (MQ=255)
          gggTTGTTTCAGCGGATCGGCAGGGTGGTGTGCTaa                          >  2:3678439/1‑36 (MQ=16)
          gggTTGTTTCAGCGGATCGGCAGGGTCGCGACCTaa                          >  2:1419256/1‑36 (MQ=25)
          gggTTGTTTCAGCGGATCGGCAGGGACGTGACCTaa                          >  1:592718/1‑36 (MQ=37)
          gggTTGTTTCAGCGGATCGGCAGGATTGTGACCTaa                          >  2:3488221/1‑36 (MQ=255)
          gggTTGTTTCAGCGGATCGGCAGGATCGTGACCTaa                          <  2:714760/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                         <  2:259066/36‑1 (MQ=255)
            gttgttTCAGCGGATCGGCAGGGTCGGGCTCTAAcc                        >  2:2237166/1‑35 (MQ=12)
            gttgttTCAGCGGATCGGCAGGATCGTGACCTAACa                        <  1:1615568/36‑1 (MQ=255)
            gttgttTCAGCGGATCGGCAGGATCGTGACCTAACa                        <  2:2236822/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                      >  2:3578344/1‑36 (MQ=21)
               gttTCAGCGGATCGGCAGGGTCGTGGCCCAAAAAgg                     >  2:1498473/1‑36 (MQ=16)
                tttCAGCGGATCGGCAGGATCGTGACCTAACAAGGt                    <  1:2910402/36‑1 (MQ=255)
                tttCAGCGGATCGGCAGGATCGTGACCTAACAAGGt                    <  1:1032712/36‑1 (MQ=255)
                   cAGCGGATCGGCAGGATCCGGACCTAAAAACGATGa                 >  2:1075097/1‑36 (MQ=16)
                     gCGGATCGGCAGGGTCGTGGCCTAACAAGGTTGAcc               >  1:2469929/1‑36 (MQ=38)
                     gCGGATCGGCAGGATCGTGACCTAACAAGGTTGAcc               <  1:2151282/36‑1 (MQ=255)
                     gCGGATCGGCAGGATCGTGACCTAACAAGGTTGAcc               <  2:2426710/36‑1 (MQ=255)
                     gCGGATCGGCAGGATCGTGACCTAACAAGGTTGAcc               <  1:845172/36‑1 (MQ=255)
                      cGGATCCGCAGGATCGTGACCGAAAAAGGATGGccc              >  2:2301725/1‑36 (MQ=16)
                       ggATCGGCAGGGTCGTGACCTAACAAGGGTCAcccc             >  1:2514571/1‑35 (MQ=25)
                       ggATCGGCAGGATCGTGACCTAACAAGGTTGACCCt             <  1:2690176/36‑1 (MQ=255)
                        gATCGGCAGGATCGTGACCTAACAAGGTTGATCCtt            >  1:1947611/1‑36 (MQ=255)
                        gATCGGCAGGATCGTGACCTAACAAGGTTGACCCtt            <  1:2134286/36‑1 (MQ=255)
                        gATCGGCAGGATCGTGACCTAACAAGGTTGACCCtt            >  2:418716/1‑36 (MQ=255)
                        gATCGGCAGGATCGTGACCTAACAAGGTTGACCCtt            >  2:120434/1‑36 (MQ=255)
                        gATCGGCAGGATCGTGACCTAAAAAGGTTGAcccct            >  1:1410584/1‑34 (MQ=38)
                        gATAGGCAGGATCGTGACCTAACAAGGTTGACCCtt            <  2:2486007/36‑1 (MQ=255)
                         aTCGGCAGGATCGTGACCTAACAAGGTTGACCCttt           <  2:1596063/36‑1 (MQ=255)
                         aTCGGCAGGATCGTGACCTAACAAGGTTGACCCttt           <  2:2635099/36‑1 (MQ=255)
                          tCGGCAGGATCGTGACCTAACAAGGTTGACCCTTTc          <  1:2364835/36‑1 (MQ=255)
                           cGGCAGGATCGTGACCTAACAAGGTTGACCCTTTCa         <  1:1828177/36‑1 (MQ=255)
                           cGGCAGGATCGTGACCTAACAAGGTTGACCCTTTCa         >  1:3314894/1‑36 (MQ=255)
                            ggCAGGATCGTGACCTAACAAGGTTGACCCTTTCac        >  1:1154886/1‑35 (MQ=255)
                            ggCAGGATCGTGACCTAACAAGGTTGACCCTTTCAg        >  1:1008814/1‑36 (MQ=255)
                            ggCAGGATCGTGACCTAACAAGGTTGACCCTTTCAg        >  2:2618341/1‑36 (MQ=255)
                               aGGATCGTGACCTAACAAGGTTGACCCTTTCAGACg     <  2:3206850/36‑1 (MQ=255)
                               aGGATCGTGACCTAACAAGGTTGACCCTTTCAGACg     <  2:3231838/36‑1 (MQ=255)
                               aGGATCGTGACCTAACAAGGTTGACCCTTTCAGACg     >  2:2381839/1‑36 (MQ=255)
                               aGGATCGTGACCTAACAAGGTTGACCCTTTCAGACg     >  1:2836643/1‑36 (MQ=255)
                               aGGATCGTGACCTAACAAGGTTGACCCTTTCAGACg     >  1:2871115/1‑36 (MQ=255)
                               aGGATCGTGACCTAACAAGGTTGACCCTTTCAGACg     <  1:3030769/36‑1 (MQ=255)
                                ggATCGTGACCTAACAAGGTTGACCCTTTCAGACGa    <  1:1867990/36‑1 (MQ=255)
                                ggATCGTGACCTAACAAGGTTGACCCTTTCAGACGa    >  1:3651842/1‑36 (MQ=255)
                                 gATCGTGACCTAACAAGGTTGACCCTTTCAGACGAt   >  2:25353/1‑36 (MQ=255)
                                 gATCGTGACCTAACAAGGTTGACCCTTTCAGACGAt   <  2:511941/36‑1 (MQ=255)
                                 gATCGTGACCTAACAAGGTTGACCCTTTCAGACGAt   <  2:1615985/36‑1 (MQ=255)
                                  aTCGTGACCTAACAAGGTTGCCCCTTTCAGACGATa  >  2:1096592/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