Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,485,4530A→G 7.9% 0.8 ‑6.0 0.074582ydcVpredicted spermidine/putrescine transporter subunit
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (7/0):  ref base (26/41):  total (33/41)
Fisher's exact test strand distribution p-value = 2.37e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCAT  >  REL606/1485418‑1485488
                                   |                                   
ttGCCGAATCTCAGTTCGGCGTTACTGGCAggagga                                     >  1:4089511/1‑36 (MQ=255)
ttGCCGAATCTCAGTTCGGCGTTACTGGCAggagga                                     >  1:3617990/1‑36 (MQ=255)
ttGCCGAATCTCAGTTCGGCGTTACTGGCAggagga                                     <  1:4536912/36‑1 (MQ=255)
 tGCCGAATCTCAGTTCGGGGTTACTGGCAGGAGGaa                                    <  1:52834/36‑1 (MQ=255)
 tGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGaa                                    <  1:2537028/36‑1 (MQ=255)
 tGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGaa                                    <  1:6337051/36‑1 (MQ=255)
  gCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGGAt                                   >  1:2086991/1‑36 (MQ=255)
  gCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAAt                                   <  1:5267661/36‑1 (MQ=255)
  gCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAAt                                   <  1:6604619/36‑1 (MQ=255)
   ccGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATg                                  <  1:5245128/36‑1 (MQ=255)
   ccGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATg                                  <  1:2166469/36‑1 (MQ=255)
   ccGAATCTCAGTGCGGCGTTACTTGCAGGAGGGATg                                  >  1:4184169/1‑36 (MQ=25)
    cGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGc                                 >  1:4024746/1‑36 (MQ=37)
    cGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGc                                 >  1:2749110/1‑36 (MQ=255)
    cGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGc                                 >  1:6804658/1‑36 (MQ=255)
     gAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCt                                <  1:4365832/36‑1 (MQ=255)
     gAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCt                                <  1:7159254/36‑1 (MQ=255)
     gAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCt                                >  1:294494/1‑36 (MQ=255)
     gAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCt                                >  1:581312/1‑36 (MQ=255)
     gAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCt                                <  1:1981544/36‑1 (MQ=255)
      aaTCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTg                               <  1:4354205/36‑1 (MQ=255)
      aaTCTCAGTTCGGCGTTACGGGCAGGAGGGATGcgg                               >  1:5886424/1‑34 (MQ=37)
       aTCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTgg                              >  1:7457379/1‑36 (MQ=255)
       aTCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTgg                              <  1:7135393/36‑1 (MQ=255)
       aTCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTgg                              <  1:4838068/36‑1 (MQ=255)
       aTCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTgg                              >  1:2003941/1‑36 (MQ=255)
        tctcAGTTCGGCGTTACTGGCAGGAGGAATGCTggg                             >  1:5372467/1‑35 (MQ=255)
        tctcAGTTCGGCGTTACTGGCAGGAGGAATGCTGGc                             >  1:1155337/1‑36 (MQ=255)
        tctcAGTTCGGCGTTACTGGCAGGAGGAATGCTGGc                             <  1:3554517/36‑1 (MQ=255)
        tctcAGTTCGGCGTTACTGGCAGGAGGAATGCTGGc                             <  1:1282649/36‑1 (MQ=255)
        tctcAGTTCGGCGTTACTGGCAGGAGGAATGCTGGc                             <  1:5050718/36‑1 (MQ=255)
        tctcAGTTCGGCGTTACTGGCAGGAGGAATGCTGGc                             <  1:5603867/36‑1 (MQ=255)
        tctcAGTTCGGCGTTACTGGCAGGAGGAATGCTGGc                             <  1:1496219/36‑1 (MQ=255)
         ctcAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCg                            <  1:4177523/36‑1 (MQ=255)
          tcAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGt                           <  1:476620/36‑1 (MQ=255)
          tcAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGt                           <  1:7213818/36‑1 (MQ=255)
          tcAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGt                           >  1:6728063/1‑36 (MQ=255)
           cAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGtt                          <  1:502917/36‑1 (MQ=255)
           cAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGtt                          <  1:3422084/36‑1 (MQ=255)
            aGTTCGGCGTTACTGGCAGGAGGGAAGGTGGCGttt                         >  1:7063368/1‑36 (MQ=25)
            aGTTCGGCGTTACTGGCAGGAGGAATGCTGNGGttt                         >  1:6047384/1‑36 (MQ=255)
            aGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGttt                         <  1:5406395/36‑1 (MQ=255)
            aGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGttt                         <  1:5904321/36‑1 (MQ=255)
             gTTCGGCGTTACTGGCAGGAGGAATGCTGGCGtttt                        >  1:5031701/1‑35 (MQ=255)
              ttCGGCGTTACTGGCAGGAGGAATGCTGGCGttttc                       >  1:5617397/1‑34 (MQ=255)
              ttCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGc                       <  1:1427370/36‑1 (MQ=255)
              ttCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGc                       >  1:4332662/1‑36 (MQ=255)
                cGGCGTTACTGGCAGGGGGGATGGCGGGGTTTTCCt                     >  1:3112235/1‑36 (MQ=9)
                cGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCt                     <  1:6997117/36‑1 (MQ=255)
                 ggCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCtt                    >  1:1500894/1‑36 (MQ=255)
                 ggCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCtt                    <  1:5643769/36‑1 (MQ=255)
                 ggCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCtt                    <  1:5809696/36‑1 (MQ=255)
                  gCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTg                   >  1:6942044/1‑36 (MQ=255)
                  gCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTg                   <  1:5613100/36‑1 (MQ=255)
                   cGTTACTGGCAGGAGGGATGCTGGCGTTTGGCTTGt                  >  1:4920544/1‑36 (MQ=37)
                   cGTTACTGGCAGGAGGAATGCTNGCGTTTGCCTTGt                  <  1:4369789/36‑1 (MQ=255)
                   cGTTACTGGCAGGAGGAATGCTGGCGTTTTCCTTGt                  >  1:6273208/1‑36 (MQ=255)
                    gTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTTTc                 >  1:1269236/1‑36 (MQ=255)
                    gTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTc                 <  1:5208655/36‑1 (MQ=255)
                    gTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTc                 <  1:7232259/36‑1 (MQ=255)
                    gTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTc                 <  1:1368337/36‑1 (MQ=255)
                     ttACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCg                <  1:2894199/36‑1 (MQ=255)
                     ttACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCg                <  1:2323906/36‑1 (MQ=255)
                       aCTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGtt              >  1:2776462/1‑36 (MQ=255)
                       aCTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGtt              <  1:5898360/36‑1 (MQ=255)
                        cTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTc             <  1:789238/36‑1 (MQ=255)
                         tGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCg            >  1:3177711/1‑36 (MQ=255)
                          ggCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGa           <  1:60356/36‑1 (MQ=255)
                          ggCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGa           <  1:3782152/36‑1 (MQ=255)
                          ggCAGGAGGAATGCTGGCGTTTGCCTTGTCGCTCGa           <  1:1657675/36‑1 (MQ=255)
                          ggCAGGAGGAATGCTGGCGTTTGCCTTGTAGTTCGa           >  1:5043109/1‑36 (MQ=255)
                          ggCAGGAGGAATGCTGGCGNTTGCCTTGTCGTTCGa           >  1:2151540/1‑36 (MQ=255)
                           gCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGAt          >  1:4383805/1‑36 (MQ=255)
                           gCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGAt          >  1:3283846/1‑36 (MQ=255)
                             aggaggAATGCTGGCGTTTGCCTTGTCGTTCGATgc        >  1:4943216/1‑35 (MQ=255)
                             aggaggAATGCTGGCGTTTGCCTTGTCGTTCGATGa        <  1:5854287/36‑1 (MQ=255)
                             aggaggAATGCTGGCGTTTGCCTTGTCGTTCGATGa        <  1:5230549/36‑1 (MQ=255)
                             aggaggAATGCTGGCGTTTGCCTTGTCGTTCGATGa        >  1:629471/1‑36 (MQ=255)
                             aggaggAATGCTGGCGTTTGCCTTGTCGTTCGATGa        <  1:1527272/36‑1 (MQ=255)
                             aggaggAATGCTGGCGTTTGCCTTGTAGTTCGGTGa        >  1:1183961/1‑36 (MQ=37)
                              ggaggaATGCTGGCGTTTGCCTTGTGGTGCGATgga       >  1:6314857/1‑34 (MQ=25)
                              ggaggaATGCTGGCGTTTGCCTTGTCGTTCGATGaa       >  1:4881872/1‑36 (MQ=255)
                              ggaggaATGCTGGCGTTTGCCTTGTCGTTCGATGaa       <  1:4168350/36‑1 (MQ=255)
                              ggaggaATGCTGGCGTTTGCCTTGTCGTTCGATGaa       >  1:574943/1‑36 (MQ=255)
                              ggaggaATGCTGGCGTTTGCCTTGTCGTTCGATGaa       >  1:5246421/1‑36 (MQ=255)
                              ggaggaATGCTGGCGTTTGCCTTGTCGTTCGATGaa       >  1:712410/1‑36 (MQ=255)
                               gaggaATGCTGGCGTTTGCCTTGTCGTTCGATGaaa      >  1:228495/1‑36 (MQ=255)
                                aggaATGCTGGCGTTTGCCTTGTCGTTCGATGAAat     >  1:7465793/1‑36 (MQ=255)
                                aggaATGCTGGCGTTTGCCTTGTCGTTCGATGAAat     >  1:2133624/1‑36 (MQ=255)
                                 ggaATGCTGGCGTTTGCCTTGTCGTTCGATGAAatc    >  1:7028409/1‑36 (MQ=255)
                                 ggaATGCTGGCGTTTGCCTTGTCGTTCGATGAAatc    >  1:3070716/1‑36 (MQ=255)
                                 ggaATGCTGGCGTTTGCCTTGTCGTTCGATGAAatc    >  1:1783935/1‑36 (MQ=255)
                                  gaATGCTGGCGTTTGCCTTGTCGTTCGATGAAatca   >  1:3136885/1‑36 (MQ=255)
                                  gaATGCTGGCGTTTGCCTTGTCGTTCGATGAAatca   >  1:4703194/1‑36 (MQ=255)
                                  gaATGCTGGCGTTTGCCTTGTCGTTCGATGAAatca   >  1:1633562/1‑36 (MQ=255)
                                  gaATGCTGGCGTTTGCCTTGTCGTTCGATGAAatca   >  1:355852/1‑36 (MQ=255)
                                  gaATGCTGGCGTTTGCCTTGTCGTGCGATGAAatca   <  1:381947/36‑1 (MQ=255)
                                   aaTGCTGGCGTTTGCCTTGTCGTTCGATGAAatcat  >  1:6917730/1‑36 (MQ=255)
                                   aaTGCTGGCGTTTGCCTTGTCGTTCGATGAAatcat  >  1:1770803/1‑36 (MQ=255)
                                   |                                   
TTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCAT  >  REL606/1485418‑1485488

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG