Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,893,5501.→G 2.9% 4.0 0.0 ‑0.570intergenic (+5/‑51)kup/insJ‑5potassium transporter/IS150 hypothetical protein
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (1/1):  ref base (33/35):  total (34/36)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.97e-01

TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTACTGCACCCATTTTGTTGGACGATGAAAT  >  REL606/3893515‑3893585
                                    |                                    
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCa                                       <  1:6025886/36‑1 (MQ=255)
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCa                                       <  1:44577/36‑1 (MQ=255)
tATCGAACTGGGCACTCAGGTCGAAATCTAACGCCa                                       >  1:4258227/1‑36 (MQ=255)
 aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..g                                    >  1:903041/1‑36 (MQ=255)
 aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..g                                    <  1:1813694/36‑1 (MQ=255)
   cGAACTGGGTACTCAGGTCGAAATCTAACGCCA..GAc                                  <  1:4619890/36‑1 (MQ=255)
   cGAACTGGGTACTCAGGTCGAAATCTAACGCCA..GAc                                  >  1:4329564/1‑36 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCCA..GACt                                 >  1:3479372/1‑36 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCCA..GACt                                 >  1:6693877/1‑36 (MQ=255)
      aCTGGGTACTCAGGTCTAAATCTAACGCCA..GACTGt                               <  1:2643515/36‑1 (MQ=255)
      aCTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGt                               >  1:5503759/1‑36 (MQ=255)
      aCTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGt                               <  1:2202490/36‑1 (MQ=255)
       cTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTa                              <  1:6799242/36‑1 (MQ=255)
       cTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTa                              <  1:3520982/36‑1 (MQ=255)
       cTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTa                              >  1:6115126/1‑36 (MQ=255)
       cTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTa                              >  1:1339050/1‑36 (MQ=255)
        tGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTAc                             <  1:3359013/36‑1 (MQ=255)
           gTACTCAGGTCGAAATCTAACGCCA..GACTGTACTGc                          >  1:2286584/1‑36 (MQ=255)
             aCTCAGGTCGAAATCTAACGCCA..GACTGTACTGCAc                        <  1:5353453/36‑1 (MQ=255)
             aCTCAGGTCGAAATCTAACGCCA..GACTGTACTGCAc                        <  1:1854745/36‑1 (MQ=255)
              cTCAGGTCGAAATCTAACGCCA..GACTGTACTTCAcc                       <  1:6620667/36‑1 (MQ=255)
              cTCAGGTCGAAATCTAACGCCA..GACTGTACTGCAcc                       >  1:7043161/1‑36 (MQ=255)
              cTCAGGTCGAAATCTAACGCCA..GACTGTACTGCAcc                       >  1:3279528/1‑36 (MQ=255)
              cTCAGGTAGAAATCTAACGCCA..GACTGTACTGCAcc                       <  1:1493422/36‑1 (MQ=255)
               tCAGGTCGAAATCTAACGCCA..GACTGTACTGCAccc                      <  1:4145591/36‑1 (MQ=255)
               tCAGGTCGAAATCTAACGCCA..GACTGTACTGCAccc                      <  1:2701724/36‑1 (MQ=255)
                cAGGTCGAAATCTAACGCCA..GACTGTACTGCACCCa                     <  1:6565609/36‑1 (MQ=255)
                cAGGTCGAAATCTAACGCCA..GACTGTACTGCACCCa                     <  1:5218593/36‑1 (MQ=255)
                 aGGTCGAAATCTAACGCCA..GACTGTACTGCACCCAt                    <  1:4210210/36‑1 (MQ=255)
                 aGGTCGAAATCTAACGCCA..GACTGTACTGCACCCAt                    <  1:5184194/36‑1 (MQ=255)
                 aGGTCGAAATCTAACGCCA..GACTGTACTGCACCCAt                    >  1:7305295/1‑36 (MQ=255)
                   gTCGAAATCTAACGCCA..GACTGTACTGCACCCAttt                  <  1:2571732/36‑1 (MQ=255)
                   gTCGAAATCTAACGCCA..GACTGTACTGCACCCAttt                  >  1:2339688/1‑36 (MQ=255)
                   gTCGAAATCTAACGCCA..GACTGTACTGCACCCAttt                  <  1:423694/36‑1 (MQ=255)
                    tCGAAATCTAACGCCA..GACTGTACTGCACCCAtttt                 >  1:1086055/1‑36 (MQ=255)
                    tCGAAATCTAACGCCA..GACTGTACTGCACCCAtttt                 <  1:2429068/36‑1 (MQ=255)
                     cGAAATCTAACGCCAG.GACTGTACTGCACCCAtttt                 >  1:6198436/1‑36 (MQ=25)
                     cGAAATCTAACGCCA..GACTGTACTGCACCCATTttg                <  1:3455184/36‑1 (MQ=255)
                     cGAAATCTAACGCCA..GACTGTACTGCACCCATTttg                <  1:4975839/36‑1 (MQ=255)
                      gAAATCTAACGCCA..GACTGTACTGCACCCATTttgt               <  1:1123347/36‑1 (MQ=255)
                      gAAATCTAACGCCA..GACTGTACTGCACCCATTttgt               >  1:6075048/1‑36 (MQ=255)
                       aaaTCTAACGCCA..GACTGTACTGCACCCATTttgtt              >  1:1538032/1‑36 (MQ=255)
                        aaTCTAACGCCA..GACTGTACTGCACCCATTttgttg             <  1:51526/36‑1 (MQ=255)
                        aaTCTAACGCCA..GACTGCACTGCACCCATTttgttg             >  1:7300889/1‑36 (MQ=255)
                         aTCTAACGCCAGGGACTGTACTGCACCCATTttgtt              <  1:6984964/36‑1 (MQ=17)
                         aTCTAACGCCA..GACTGTACTGCACCCATTTTGTTgg            >  1:5676672/1‑36 (MQ=255)
                          tCTAACGCCA..GACTGTACTGCACCCATTTTGTTGGa           >  1:5291489/1‑36 (MQ=255)
                          tCTAACGCCA..GACTGTACTGCACCCATTTTGTTGGa           <  1:3253005/36‑1 (MQ=255)
                          tCTAACGCCA..GACTGTACTGCACCCATTTTGTTGGa           >  1:4419196/1‑36 (MQ=255)
                           cTAACGCCA..GACTGTACTGCACCCATTTTGTTGGAc          <  1:4358178/36‑1 (MQ=255)
                             aaCGCCA..GACTGTACTGCACCCATTTTTTTGGACGa        >  1:2136078/1‑36 (MQ=255)
                             aaCGCCA..GACTGTACTGCACCCATTTTGTTGGACGa        <  1:5912004/36‑1 (MQ=255)
                              aCGCCA..GACTGTACTGCACCCATTTTGTTGGACGAt       >  1:32132/1‑36 (MQ=255)
                              aCGCCA..GACTGTACTGCACCCATTTTGTTGGACGAt       >  1:4678343/1‑36 (MQ=255)
                              aCGCCA..GACTGTACTGCACCCATTTTGTTGGACGAt       <  1:2988156/36‑1 (MQ=255)
                               cGCCA..GACTGTACTGCACCCATTTTGTTGGACGATg      >  1:5065330/1‑36 (MQ=255)
                               cGCCA..GACTGTACTGCACCCATTTTGTTGGACGATg      >  1:3168875/1‑36 (MQ=255)
                               cGCCA..GACTGTACTGCACCCATTTTGTTGGACGATg      >  1:6500074/1‑36 (MQ=255)
                               cGCCA..GACTGTACTGCACCCATTTTGTTGGACGATg      <  1:1085652/36‑1 (MQ=255)
                                gCCA..GACTGTACTGCACCCATTTTGTTGGACGATGa     <  1:1079770/36‑1 (MQ=255)
                                 ccA..GACTGTCCTGCACCCATTGTGTTGGACGATGaa    <  1:3930709/36‑1 (MQ=17)
                                 ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa    >  1:5872723/1‑36 (MQ=255)
                                 ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa    <  1:2778370/36‑1 (MQ=255)
                                 ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa    >  1:452964/1‑36 (MQ=255)
                                 ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa    >  1:4008109/1‑36 (MQ=255)
                                 ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa    <  1:315520/36‑1 (MQ=255)
                                 ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa    >  1:681215/1‑36 (MQ=255)
                                  cA..GACTGTACTGCCCCCATTTTGTTGGACGATGaaa   <  1:5518838/36‑1 (MQ=255)
                                  cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa   <  1:1868054/36‑1 (MQ=40)
                                  cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa   >  1:2105405/1‑36 (MQ=40)
                                  cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa   >  1:2823119/1‑36 (MQ=40)
                                  cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa   <  1:6644060/36‑1 (MQ=40)
                                  cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa   >  1:542311/1‑36 (MQ=40)
                                  cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa   >  1:5546571/1‑36 (MQ=40)
                                  cA..GACTGTACTGCACCCATTTTGTTGGACGATGaaa   >  1:5960908/1‑36 (MQ=40)
                                   a..GACTGTACTGCACCCATTTTGTTGGGCGATGAAAt  >  1:7389806/1‑36 (MQ=255)
                                   a..GACTGTACTGCACCCATTTTGTTGGACGATGAAAt  <  1:3422945/36‑1 (MQ=38)
                                   a..GACTGTACTGCACCCATTTTGTTGGACGATGAAAt  >  1:3878815/1‑36 (MQ=38)
                                   a..GACTGTACTGCACCCATTTTGTTGGACGATGAAAt  <  1:6047420/36‑1 (MQ=38)
                                    |                                    
TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTACTGCACCCATTTTGTTGGACGATGAAAT  >  REL606/3893515‑3893585

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG