Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,369,8661.→G 4.8% 3.5 ‑2.8 0.063coding (138/390 nt)rpsK30S ribosomal protein S11
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (22/38):  total (25/38)
Fisher's exact test strand distribution p-value = 5.79e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GATTTGCGAGAACCACGGAAACCGGAACCACCGGCT.GTTGCCCAACCCAACGCGTTACCCTGACGATCAGT  >  REL606/3369831‑3369901
                                    |                                   
gATTTGCGAGAACCACGGAAACCGGAACCACCGGCt                                      >  1:1062321/1‑36 (MQ=255)
gATTTGCGAGAACCACGGAAACCGGAACCACCGGCt                                      >  1:3040418/1‑36 (MQ=255)
gATTTGCGAGAACCACGGAAACCGGAACCACCGGCt                                      <  1:601451/36‑1 (MQ=255)
 aTTTGCGAGAACCACGGAAACCGGAACCACCGGCT.g                                    <  1:3236394/36‑1 (MQ=255)
 aTTTGCGAGAACCACGGAAACCGGAACCACCGGCT.g                                    <  1:2835015/36‑1 (MQ=255)
 aTTTGCGAGAACCACGGAAACCGGAACCACCGGCT.g                                    <  1:649459/36‑1 (MQ=255)
 aTTTGCGAGAACCACGGAAACCGGAACCACCGGCT.g                                    <  1:2503717/36‑1 (MQ=255)
  tttGCGAGACCCACGGAAACCGGAACCACCGGCT.Gt                                   <  1:2864632/36‑1 (MQ=255)
  tttGCGAGAACCACGGAAACCGGAACCACCGGCT.gg                                   >  1:5190697/1‑35 (MQ=255)
   ttGCGAGAACCACGGAAACCGGAACCACCGGCT.ggt                                  >  1:1089697/1‑34 (MQ=255)
   ttGCGAGAACCACGGAAACCGGAACCACCGGCT.ggt                                  >  1:492362/1‑34 (MQ=255)
   ttGCGAGAACCACGGAAACCGGAACCACCGGCT.ggt                                  >  1:3708655/1‑34 (MQ=255)
   ttGCGAGAACCACGGAAACCGGAACCACCGGCT.ggt                                  >  1:1711575/1‑34 (MQ=255)
   ttGCGAGAACCACGGAAACCGGAACCACCGGCT.Gtt                                  <  1:2545415/36‑1 (MQ=255)
   ttGCGAGAACCACGGAAACCGGAACCACCGGCT.Gtt                                  >  1:7536861/1‑36 (MQ=255)
   ttGCGACAACCACTGAAACCAGAACCACCGGCT.Gtt                                  <  1:2226042/36‑1 (MQ=25)
    tGCGAGAACCACGGAAACCGGAACCACCGGCT.GTTg                                 <  1:2110040/36‑1 (MQ=255)
    tGCGAGAACCACGGAAACCGGAACCACCGGCT.GGTg                                 >  1:7220869/1‑36 (MQ=255)
     gCGAGAACCACGGAAACCGGAACCACCGGCT.GTTGc                                <  1:4847609/36‑1 (MQ=255)
      cGAGAACCACGGAAACCGGAACCACCGGCT.GTTGcc                               <  1:3233149/36‑1 (MQ=255)
      cGAGAACCACGGAAACCGGAACCACCGGCT.GGTggc                               >  1:6281746/1‑34 (MQ=38)
       gagaACCACGGAAACCGGAACCACCGGCTGGTTGcc                               >  1:3869945/1‑36 (MQ=37)
       gagaACCACGGAAACCGGAACCACCGGCT.GTTGccc                              <  1:7281480/36‑1 (MQ=255)
       gagaACCACGGAAACCGGAACCACCGGCT.GTTGccc                              <  1:3374064/36‑1 (MQ=255)
       gagaACCACGGAAACCGGAACCACCGGCT.GTTGccc                              <  1:5755020/36‑1 (MQ=255)
       gagaACCACGGAAACCGGAACCACCGGCT.GGTgggc                              >  1:253822/1‑33 (MQ=37)
       gagaACCACGGAAACCGGAACCACCGGCT.GGTGgcc                              >  1:2798554/1‑36 (MQ=37)
       gagaACCACGGAAACCGGAACCACCGGCT.GGGGccc                              >  1:7040780/1‑36 (MQ=37)
        agaACCACGGAAACCGGAACCACCGGCT.GGTGgccc                             >  1:374119/1‑35 (MQ=37)
        agaACCACGGAAACCGGAACCACCGGCT.GGTGgccc                             >  1:7640119/1‑35 (MQ=37)
         gaACCACGGAAACCGGAACCACCGGCT.GTTGcccaa                            <  1:142914/36‑1 (MQ=255)
         gaACCACGGAAACCGGAACCACCGGCT.GTTGcccaa                            <  1:4264597/36‑1 (MQ=255)
          aaCCACGGAAACCGGAACCACCGGCTGGTGGcccaa                            >  1:6701119/1‑36 (MQ=25)
          aaCCACGGAAACCGGAACCACCGGCT.GTTGcccaac                           <  1:3435631/36‑1 (MQ=255)
          aaCCACGGAAACCGGAACCACCGGCT.GGTTccccac                           >  1:1431866/1‑36 (MQ=25)
           aCCACTGAAACCGGAACCACCGGCT.TTTGcccaacc                          <  1:2795069/36‑1 (MQ=37)
           aCCACGGAAACCGGAACCACCGGCT.GTTGcccaacc                          <  1:6159247/36‑1 (MQ=255)
             cACGGAAACCGGAACCACCGGCTGGTTGcccaccca                         >  1:6607893/1‑35 (MQ=25)
             cACGGAAACCGGAACCACCGGCT.GTTGcccaaccca                        <  1:6263321/36‑1 (MQ=255)
               cGGAAACCGGAACCACCGGCT.GTTTctccacccccc                      >  1:3396065/1‑33 (MQ=19)
               cGGAAACCGGAACCACCGGCT.GGTGGccaaccccac                      >  1:629487/1‑36 (MQ=25)
                ggAAACCGGAACCACCGGCT.GTTGCCCAACCCAAcg                     <  1:3567423/36‑1 (MQ=255)
                ggAAACCGGAACCACCGGCT.GTTGCCCAACCCAAcg                     <  1:6366466/36‑1 (MQ=255)
                ggAAACCGGAACCACCGGCT.GTTGCCCAACCCAAcg                     <  1:7321224/36‑1 (MQ=255)
                ggAAACCGGAACCAACGGCT.GGTGGccaacccaaca                     >  1:4452172/1‑35 (MQ=25)
                 gAAACCGGAACCACCGGCT.GTTGCCCAACCCAAcgc                    <  1:4398495/36‑1 (MQ=255)
                 aaaaCCGGAACCACCGGCT.GTTGCCCAACCCAAcgc                    <  1:2647556/35‑1 (MQ=255)
                  aaaCCGGAACCACCGGCT.GTTGCCCAACCCAAcgcg                   <  1:5137226/36‑1 (MQ=255)
                   aaCCGGAACCACCGGCT.GTTGCCCAACCCAACGCGt                  <  1:1446389/36‑1 (MQ=255)
                    aCCGGAACCACCGGCT.GTTGCCCAACCCAACGCGtt                 <  1:5730726/36‑1 (MQ=255)
                      cGGAACCACCGGCT.GTTGCCCAACCCCACCCGTCAc               >  1:1737102/1‑36 (MQ=37)
                       ggAACCACCGGCT.GTTGCCCCACCCCACGCCGTAcc              >  1:2970484/1‑36 (MQ=19)
                       ggAACCACCGGCT.GTTGCCCAACCCAACGCGTTAcc              <  1:7823464/36‑1 (MQ=255)
                       ggAACCACCGGCT.GTTGCCCAACCCAACGCGTTAcc              <  1:4604145/36‑1 (MQ=255)
                       ggAACCACCGGCT.GGTGCCCCAACCCACCCGTTAcc              >  1:2254318/1‑36 (MQ=12)
                        gAACCACCGGCT.GTTGCCCAACCCAACGCGTTAccc             <  1:5513277/36‑1 (MQ=255)
                        gAACCACCGGCT.GTTGCCCAACCCAACGCGTTAccc             <  1:185215/36‑1 (MQ=255)
                         aaCCACCGGCT.GTTGCCCAACCCAACGCGTTACCCt            <  1:304607/36‑1 (MQ=255)
                         aaCCACCGGCT.GTTGCCCAACCCAACGCGTTACCCt            <  1:4672883/36‑1 (MQ=255)
                          accaccGGCT.GTTGCCCAACCCCACGCGTTACCCTg           >  1:4819230/1‑36 (MQ=255)
                          accaccGGCT.GTTGCCCAACCCAACGCGTTACCCTg           <  1:4805006/36‑1 (MQ=255)
                           ccaccGGCT.GTTGCCCAACCCAACGCGTTACCCTGa          <  1:3546594/36‑1 (MQ=255)
                            caccGGCT.GTTGCCCAACCCAACGCGTTACCCTGAc         <  1:2319098/36‑1 (MQ=255)
                            caccGGCT.GTTGCCCAACCCAACGCGTTACCCTGAc         <  1:4532716/36‑1 (MQ=255)
                             accGGCT.GTTGCCCAACCCAACGCGTTACCCTGACg        <  1:7243759/36‑1 (MQ=255)
                             accGGCT.GTTGCCCAACCCAACGCGTTACCCTGACg        <  1:7316621/36‑1 (MQ=255)
                               cGGCT.GTTGCCCAACCCAACGCGTTACCCTGACGAt      >  1:6594164/1‑36 (MQ=255)
                               cGGCT.GTTGCCCAACCCAACGCGTTACCCTGACGAt      <  1:4334675/36‑1 (MQ=255)
                               cGGCT.GTTGCCCAACCCAACGCGTTACCCTGACGAt      <  1:257300/36‑1 (MQ=255)
                                ggCT.GTTGCCCAACCCAACGCGTTACCCTGACGATc     >  1:5262285/1‑36 (MQ=255)
                                ggCT.GTTGCCCAACCCAACGCGTTACCCTGACGATc     >  1:4263619/1‑36 (MQ=255)
                                ggCT.GTTGCCCAACCCAACGCGTTACCCTGACGATc     >  1:2805122/1‑36 (MQ=255)
                                 gCT.GTTGCCCAACCCAACGCGTTACCCTGACGATCa    >  1:7361149/1‑36 (MQ=255)
                                  cT.GTTGCCCAACCCAACGCGTTACCCTGACGATCAg   <  1:7781840/36‑1 (MQ=255)
                                  cT.GTTGCCCAACCCAACGCGTTACCCTGACGATCAg   >  1:7929490/1‑36 (MQ=255)
                                   t.GTTGCCCAACCCAACGCTTTACTCTGTCGATCAGt  <  1:2779217/36‑1 (MQ=18)
                                   t.GTTGCCCAACCCAACGCGTTACCCTGACGATCAGt  <  1:5004128/36‑1 (MQ=255)
                                   t.GTTGCCCAACCCAACGCGTTACCCTGACGATCAGt  <  1:4665354/36‑1 (MQ=255)
                                   t.GTTGCCCAACCCAACGCGTTACCCTGACGATCAGt  >  1:7224747/1‑36 (MQ=255)
                                   t.GTTGCCCAACCCAACGCGTTACCCTGACGATCAGt  >  1:1942597/1‑36 (MQ=255)
                                   t.GTTGCCCAACCCAACGCGTTACCCTGACGATCAGt  <  1:1286848/36‑1 (MQ=255)
                                    |                                   
GATTTGCGAGAACCACGGAAACCGGAACCACCGGCT.GTTGCCCAACCCAACGCGTTACCCTGACGATCAGT  >  REL606/3369831‑3369901

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG