Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,485,8930G→T 9.9% 2.6 ‑9.5 0.078787rtcBhypothetical protein
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (9/0):  ref base (21/48):  total (30/48)
Fisher's exact test strand distribution p-value = 7.85e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCA  >  REL606/3485858‑3485928
                                   |                                   
cATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTg                                     <  1:4720394/36‑1 (MQ=255)
cATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTg                                     <  1:3074282/36‑1 (MQ=255)
 aTCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGc                                    <  1:4517962/36‑1 (MQ=255)
  tCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCg                                   <  1:497388/36‑1 (MQ=255)
  tCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCg                                   <  1:6604081/36‑1 (MQ=255)
   cGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTTCgc                                  >  1:4476701/1‑35 (MQ=255)
   cGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGt                                  <  1:1350235/36‑1 (MQ=255)
   cGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGt                                  <  1:2218448/36‑1 (MQ=255)
   cGCCAGCGTTTGTGGCTGGCTGTCCGTTTTTTGTGt                                  >  1:7222408/1‑36 (MQ=9)
    gCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTa                                 <  1:6671735/36‑1 (MQ=255)
    gCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTa                                 <  1:7762456/36‑1 (MQ=255)
     ccAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTaa                                <  1:3052814/36‑1 (MQ=255)
       aGCGTTTGTGGCTGTCTTGCCGTTTTTTTTGTAAtc                              >  1:4324457/1‑35 (MQ=12)
        gCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGc                             <  1:7277385/36‑1 (MQ=255)
        gCGTTTGTGGCTGGCTGTCCGTTTTCTTCGTAATgt                             >  1:4686068/1‑35 (MQ=18)
         cGTTTGTGGCTGTTTTACCGTTTTCTTCGTTATGct                            >  1:946815/1‑36 (MQ=19)
         cGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGct                            <  1:7175212/36‑1 (MQ=255)
          gTTTGTGGCTGTCTTACCGTTTTTTTCGTAATGctt                           >  1:6073154/1‑35 (MQ=25)
          gTTTGTGGCTGTCTGACCGTTTTCTTCGTAATTctt                           >  1:4013969/1‑35 (MQ=37)
          gTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGctc                           <  1:754178/36‑1 (MQ=255)
           tttGTGGCTGTCTGACCGTTTTTTTCGTAATGctct                          >  1:1630881/1‑36 (MQ=37)
           tttGTGGCTGTCTGACCGTTTTCTGCGTAATGctct                          <  1:3750011/36‑1 (MQ=255)
           tttGTGGCTGTCTGACCGTTTTCTGCGTAATGctct                          <  1:6471695/36‑1 (MQ=255)
            ttGTGGCTGTCTGACCGTTTTTTTCGTTATTctttt                         >  1:5690625/1‑33 (MQ=12)
            ttGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTg                         <  1:1958322/36‑1 (MQ=255)
            ttGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTg                         <  1:6192724/36‑1 (MQ=255)
             tgtgGCTGTCTGACCGTTTTTTTCGTAATTctcttc                        >  1:7602962/1‑34 (MQ=21)
             tgtgGCTGTCTGACCGTTTTCTGCGTAATGCTCTGc                        <  1:5842972/36‑1 (MQ=255)
              gtgTCTGTCTGCCCGTTTTCTGCGTAATGCTCTGCa                       <  1:706937/36‑1 (MQ=37)
              gtgGCTGTCTGACCGTTTTTTGCGTAATGCTCTGCa                       >  1:1848686/1‑36 (MQ=255)
              gtgGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCa                       <  1:5363104/36‑1 (MQ=255)
               ttgCTGTCTGACCGTTTTCTGCGTAATGCTCTGCaa                      <  1:3372904/34‑1 (MQ=255)
               tgGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCaa                      <  1:719587/36‑1 (MQ=255)
               tgGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCaa                      >  1:2665072/1‑36 (MQ=255)
                ggCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAAt                     <  1:6053273/36‑1 (MQ=255)
                ggCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAAt                     <  1:971776/36‑1 (MQ=255)
                ggCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAAt                     <  1:2537304/36‑1 (MQ=255)
                ggCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAAt                     >  1:3452403/1‑36 (MQ=255)
                 gCTGTCTGACCGTTTTCTGCGTAATGCTCTTCAATg                    >  1:6238076/1‑36 (MQ=255)
                 gCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATg                    <  1:5464195/36‑1 (MQ=255)
                 gCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATg                    >  1:907555/1‑36 (MQ=255)
                  cTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGc                   <  1:1765254/36‑1 (MQ=255)
                  cTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGc                   <  1:1999494/36‑1 (MQ=255)
                  cTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGc                   >  1:575327/1‑36 (MQ=255)
                   tGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGcc                  <  1:2159979/36‑1 (MQ=255)
                   tGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGcc                  >  1:629078/1‑36 (MQ=255)
                   tGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGcc                  <  1:3284010/36‑1 (MQ=255)
                   tGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGcc                  <  1:4944220/36‑1 (MQ=255)
                   tGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGcc                  <  1:4192954/36‑1 (MQ=255)
                    gTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCg                 <  1:3233484/36‑1 (MQ=255)
                     tCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGt                >  1:5444072/1‑36 (MQ=255)
                      cTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGtt               <  1:1821618/36‑1 (MQ=255)
                      cTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGtt               <  1:7018022/36‑1 (MQ=255)
                      cTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGtt               >  1:7224733/1‑36 (MQ=255)
                       tGACCGTTTTCCGCGTAATGCTCTGCAATGCCGTTa              >  1:6022008/1‑36 (MQ=255)
                        gACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTAc             >  1:910407/1‑36 (MQ=255)
                        gACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTAc             >  1:1772090/1‑36 (MQ=255)
                        gACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTAc             >  1:5502921/1‑36 (MQ=255)
                         aCCGTTTTCTGCGTAATGCTCTGCAATGCCGTTAcc            >  1:1830794/1‑36 (MQ=255)
                         aCCGTTTTCTGCGTAATGCTCTGCAATGCCGTTAcc            <  1:5432951/36‑1 (MQ=255)
                         aCCGTTTTCTGCGTAATGCTCTGCAATGCCGTTAcc            >  1:7789920/1‑36 (MQ=255)
                          ccGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCa           <  1:7832776/36‑1 (MQ=255)
                           cGTTTTCTGCGTAATGCTCTGCGATGCCGTTACCAc          <  1:918843/36‑1 (MQ=255)
                           cGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCAc          <  1:1757250/36‑1 (MQ=255)
                           cGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCAc          <  1:1446754/36‑1 (MQ=255)
                            gTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACg         <  1:6366744/36‑1 (MQ=255)
                            gTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACg         >  1:6960801/1‑36 (MQ=255)
                             ttttCTGCGTAATGCTCTGCAATGCCGTTACCACGt        >  1:6743663/1‑36 (MQ=255)
                              tttCTGCGTAATGCTCTGCAATGCCGTTACCACGtt       >  1:931578/1‑36 (MQ=255)
                              tttCTGCGTAATGCTCTGCAATGCCGTTACCACGtt       >  1:2200549/1‑36 (MQ=255)
                              tttCTGCGTAATGCTCTGCAATGCCGTTACCACGtt       >  1:112941/1‑36 (MQ=255)
                               ttCTGCGTAATGCTCTGCAATGCCGTTACCACGttt      <  1:3732630/36‑1 (MQ=255)
                               ttCTGCGTAATGCTCTGCAATGCCGTTACCACGttt      <  1:3580853/36‑1 (MQ=255)
                                tCTGCGTAATGCTCTGCAATGCCGTTACCACGtttt     <  1:664852/36‑1 (MQ=255)
                                 cTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTc    <  1:4598820/36‑1 (MQ=255)
                                 cTGCATAATGCTCTGCAATGCCGTTACCACGTTTTc    <  1:5671852/36‑1 (MQ=255)
                                  tGCGTAATGCTCTGCAATGCCGTTACCACGTTTTcc   <  1:5743312/36‑1 (MQ=255)
                                  tGCGTAATGCTCTGCAATGCCGTTACCACGTTTTcc   <  1:2878430/36‑1 (MQ=255)
                                  tGCGTAATGCTCTGCAATGCCGTTACCACGTTTTcc   <  1:1445827/36‑1 (MQ=255)
                                  tGCGTAATGCTCTGCAATGCCGTTACCACGTTTTcc   <  1:982168/36‑1 (MQ=255)
                                   gCGTAATGCTCTGCAATGCCGTTACCACGTTTTCca  <  1:3218567/36‑1 (MQ=255)
                                   gCGTAATGCTCTGCAATGCCGTTACCACGTTTTCca  <  1:2583479/36‑1 (MQ=255)
                                   |                                   
CATCGCCAGCGTTTGTGGCTGTCTGACCGTTTTCTGCGTAATGCTCTGCAATGCCGTTACCACGTTTTCCA  >  REL606/3485858‑3485928

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG