Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,485,8930G→T10.7% 1.4 ‑10.9 0.077787rtcBhypothetical protein
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (11/0):  ref base (21/45):  total (32/45)
Fisher's exact test strand distribution p-value = 1.93e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCA  >  REL606/3485858‑3485928
                                   |                                   
cATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTg                                     <  1:125521/36‑1 (MQ=255)
 aTCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGc                                    <  1:2732447/36‑1 (MQ=255)
  tCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCg                                   <  1:4088704/36‑1 (MQ=255)
   cGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGt                                  <  1:345626/36‑1 (MQ=255)
   cGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGt                                  <  1:5409334/36‑1 (MQ=255)
   cGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGt                                  <  1:2055693/36‑1 (MQ=255)
    gCCAGCGTTTGTGGCTGTCTTACCGTTTTTTTCGtt                                 >  1:364834/1‑35 (MQ=18)
    gCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTa                                 <  1:3926898/36‑1 (MQ=255)
     ccAGCGTTTGTGGCTGTCTGACCGTTTTCTTCTTaa                                >  1:7422704/1‑36 (MQ=21)
     ccAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTaa                                <  1:5367966/36‑1 (MQ=255)
     ccAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTaa                                <  1:6338420/36‑1 (MQ=255)
      cAGCGTTTGTGGCTGTTGTTCCGTTTTCTTCTTAAt                               >  1:7431697/1‑36 (MQ=9)
      cAGCGTTTGTGGCTGTCTTACCGTTTTTTTGGTAAt                               >  1:4345664/1‑36 (MQ=16)
      cAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAAt                               <  1:1834668/36‑1 (MQ=255)
      cAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAAt                               <  1:6203497/36‑1 (MQ=255)
      cAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAAt                               <  1:4707185/36‑1 (MQ=255)
       aGCGTTTGTGGCTGTTTGAGCGTTTTTTTTGTAAtt                              >  1:1114111/1‑35 (MQ=9)
       aGCGTTTGTGGCTGTGTGACCGTTTTTTTGGTAAtt                              >  1:5503342/1‑35 (MQ=16)
       aGCGTTTGTGGCTGTCTTACCGTTTTTTTCGTTAtt                              >  1:3060104/1‑35 (MQ=16)
       aGCGTTTGTGGCTGTCTGGCCCGTTTTTTCGTAAtt                              >  1:2206093/1‑35 (MQ=12)
       aGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATg                              <  1:816757/36‑1 (MQ=255)
       aGCGTTTGTGGCTATCTTACCGTTTTTTTCGTAATg                              >  1:4368633/1‑36 (MQ=19)
        gCGTTTGTGTCTGTCTGACCGTTTTCTGCGTAATGc                             <  1:5208306/36‑1 (MQ=255)
         cGTTTGTGGCTGTCTTACCGTTTTTTTCGTAATTct                            >  1:533454/1‑36 (MQ=19)
         cGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGct                            <  1:5823606/36‑1 (MQ=255)
         cGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGct                            <  1:586546/36‑1 (MQ=255)
          gTTTGTGGCTGTCTGGTCGGTTTTTGCGGCATGctc                           >  1:5832427/1‑36 (MQ=9)
          gTTTGTGGCTGTCTGAACGTTTTCTTCGTTATTctc                           >  1:4142136/1‑36 (MQ=12)
           tttGTGGCTGTCTGACCGTTTTCTGCGTAATGctct                          <  1:5262178/36‑1 (MQ=255)
            ttGTGGCTGTCTGCCCGTTTTCTGCGTAATGCTCTg                         <  1:5228296/36‑1 (MQ=255)
            ttGTGGCTGTCTGACCGGTTTCTGAGTGATGCTCTg                         >  1:4260587/1‑36 (MQ=25)
              gtgGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCa                       <  1:6972858/36‑1 (MQ=255)
               tgGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCaa                      <  1:7270404/36‑1 (MQ=255)
                ggCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAAt                     <  1:3428248/36‑1 (MQ=255)
                ggCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAAt                     >  1:3408636/1‑36 (MQ=255)
                ggCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAAt                     <  1:135492/36‑1 (MQ=255)
                 gCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATg                    <  1:1232522/36‑1 (MQ=255)
                 gCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATg                    <  1:6823307/36‑1 (MQ=255)
                  cTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGc                   >  1:1310862/1‑36 (MQ=255)
                   tGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGcc                  <  1:758055/36‑1 (MQ=255)
                    gTCTGACCGTTTTCTGCTTAATTCTCTGCAATGCCg                 >  1:810214/1‑36 (MQ=37)
                    gTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCg                 >  1:825528/1‑36 (MQ=255)
                     tCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGt                >  1:666303/1‑36 (MQ=255)
                     tCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGt                <  1:5888142/36‑1 (MQ=255)
                     tCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGt                >  1:2886657/1‑36 (MQ=255)
                     tCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGt                <  1:6843479/36‑1 (MQ=255)
                       tGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTa              <  1:5054231/36‑1 (MQ=255)
                       tGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTa              >  1:6023197/1‑36 (MQ=255)
                       tGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTa              >  1:7194093/1‑36 (MQ=255)
                        gACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTAc             <  1:1576010/36‑1 (MQ=255)
                        gACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTAc             <  1:5610439/36‑1 (MQ=255)
                        gACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTAc             >  1:4782824/1‑36 (MQ=255)
                         aCCGTTTTCTGCGTAATGCTCTGCAATGCCGTTAcc            >  1:2276766/1‑36 (MQ=255)
                         aCCGTTTTCTGCGTAATGCTCTGCAATGCCGTTAcc            <  1:4957665/36‑1 (MQ=255)
                         aCCGTTTACTGCGTAATGCTCTGCAATGCCGTTAcc            <  1:4676442/36‑1 (MQ=255)
                           cGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCAc          >  1:2811144/1‑36 (MQ=255)
                           cGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCAc          <  1:711010/36‑1 (MQ=255)
                           cGTTTTCTGAGTAATGCTCTGCAATGCCGTTACCAc          <  1:4147058/36‑1 (MQ=255)
                            gTTTTCTGCGTAATGCTCTGCAATGCCGTTCCCACg         <  1:4510331/36‑1 (MQ=255)
                            gTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACg         <  1:7046471/36‑1 (MQ=255)
                             ttttCTGCGTAATGCTCTGCAATGCCGTTACCACGt        <  1:3802613/36‑1 (MQ=255)
                             ttttCTGCGTAATGCTCTGCAATGCCGTTACCACGt        >  1:2246332/1‑36 (MQ=255)
                             ttttCTGCGTAATGCTCTGCAATGCCGTTACCACGt        <  1:974373/36‑1 (MQ=255)
                              tttCTGCGTAATGCTCTGCAATGCCGTTACCACGtt       <  1:4931649/36‑1 (MQ=255)
                              tttCTGCGTAATGCTCTGCAATGCCGTTACCACGtt       >  1:1166545/1‑36 (MQ=255)
                               ttCTGCGTAATGCTCTGCAATGCCGTTACCACGttt      <  1:6524067/36‑1 (MQ=255)
                               ttCTGCGTAATGCTCTGCAATGCCGTTACCACGttt      <  1:3486634/36‑1 (MQ=255)
                               ttCTGCGTAATGCTCTGCAATGCCGTTACCACGttt      <  1:2344557/36‑1 (MQ=255)
                                tCTGCGTAATGCTCTGCAATGCCGTTACCACGtttt     >  1:5061787/1‑36 (MQ=255)
                                tCTGCGTAATGCTCTGCAATGCCGTTACCACGtttt     <  1:7291614/36‑1 (MQ=255)
                                 cTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTc    <  1:6716202/36‑1 (MQ=255)
                                 cTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTc    >  1:5713514/1‑36 (MQ=255)
                                 cTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTc    <  1:2223498/36‑1 (MQ=255)
                                  tGCGTAATGCTCTGCAATGCCGTTACCACGTTTTcc   >  1:6499069/1‑36 (MQ=255)
                                  tGCGTAATGCTCTGCAATGCCGTTACCACGTTTTcc   <  1:5990803/36‑1 (MQ=255)
                                  tGCGTAATGCTCTGCAATGCCGTTACCACGTTTTcc   >  1:3313662/1‑36 (MQ=255)
                                  tGCGTAATGCTCTGCAATGCCGTTACCACGTTTTcc   >  1:1010425/1‑36 (MQ=255)
                                  tGCGTAATGCTCTGAAATGCCGTTACCACGTTTTcc   >  1:3104702/1‑36 (MQ=255)
                                   gCGTAATGCTCTGCAATGCCGTTACCACGTTTTCca  <  1:2746242/36‑1 (MQ=255)
                                   |                                   
CATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCA  >  REL606/3485858‑3485928

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 22 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG