Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,724,2660C→A14.1% 12.3 ‑15.1 0.096360srlRDNA‑bindng transcriptional repressor
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (1/14):  ref base (62/18):  total (63/33)
Fisher's exact test strand distribution p-value = 2.74e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

CTGCTCTCGCGCTTTAATAACATCACGGTGATGACCAACAGCCTGCATATCGTCAATGCGCTATCCGAACT  >  REL606/2724231‑2724301
                                   |                                   
cTGCTCTCGCGCTTTAATAACATCACGGTGATGAcc                                     >  1:936640/1‑36 (MQ=255)
 tGCTCTCGCGCTTTAATAACATCGCGGTGATGTCCa                                    <  1:3777010/36‑1 (MQ=37)
  gCTCTCGCGCTTTAATAACATCACGGTGATGACCaa                                   >  1:4965462/1‑36 (MQ=255)
  gCTCTCGCGCTTTAATAACATCACGGTGATGACCaa                                   <  1:2666672/36‑1 (MQ=255)
    tctcGCGCTTTAATAACATCACGGTGATGACCAACa                                 <  1:1956044/36‑1 (MQ=255)
      tcGCGCTTTAATAACATCACGGTGATGACCAACAGc                               <  1:4651736/36‑1 (MQ=255)
      tcGCGCTTTAATAACATCACGGTGATGACCAACAGc                               >  1:472750/1‑36 (MQ=255)
      tcGCGCTTTAATAACATCACGGTGATGACCAACAGc                               <  1:4344605/36‑1 (MQ=255)
      tcGCGCTTTAATAACATCACGGTGATGACCAACAGc                               <  1:4295937/36‑1 (MQ=255)
       cgcgCTTTAATAACATCACGGTGATGACCAACAGcc                              >  1:1176235/1‑36 (MQ=255)
       cgcgCTTTAATAACATCACGGTGATGACCAACAGcc                              >  1:1484763/1‑36 (MQ=255)
        gcgcTTTAATAACATCACGGTGATGACCAACAGCCt                             >  1:3105044/1‑36 (MQ=255)
        gcgcTTTAATAACATCACGGTGATGACCAACAGCCt                             >  1:5129803/1‑36 (MQ=255)
        gcgcTTTAATAACATCACGGTGATGACCAACAGCCt                             <  1:2640188/36‑1 (MQ=255)
        gcgcTTTAATAACATCACGGTGATGACCAACAGCCt                             <  1:5708093/36‑1 (MQ=255)
         cgcTTTAATAACATCACGGGGATGACCAACAGCCTg                            >  1:204797/1‑36 (MQ=255)
          gcTTTAATAACATCACGGTGATGACCAACAGCCTGc                           >  1:1312581/1‑36 (MQ=255)
          gcTTTAATAACATCACGGTGATGACCAACAGCCTGc                           <  1:4768676/36‑1 (MQ=255)
          gcTTTAATAACATCACGGTGATGACCAACAGCCTGc                           >  1:5865424/1‑36 (MQ=255)
          gcTTTAATAACATCACGGTGATGACCAACAGCCTGc                           >  1:4932494/1‑36 (MQ=255)
          gcTTTAATAACATCACGGTGATGACCAACAGCCTGc                           >  1:1441797/1‑36 (MQ=255)
          gcTTTAATAACATCACGGTGATGACCAACAGCCTGc                           >  1:1650816/1‑36 (MQ=255)
          gcTTTAATAACATCACGGTGATGACCAACAGCCTGc                           >  1:5208717/1‑36 (MQ=255)
          gcTTTAATAACATCACGGCGATGACCAACAGCCTGc                           >  1:4472223/1‑36 (MQ=255)
            tttGATAACATCACGGTGATGACCAACAGCCTGCat                         <  1:4644102/36‑1 (MQ=255)
             ttAATAACATCAGGGTGATGACCAACAGCCTGCata                        >  1:4717168/1‑36 (MQ=255)
             ttAATAACATCACGTTCATGACCAACAGCCTGCata                        <  1:534745/36‑1 (MQ=37)
             ttAATAACATCACGGTGATGACCAACAGCCTGCata                        >  1:3443581/1‑36 (MQ=255)
             ttAATAACATCACGGTGATGACCAACAGCCTGCata                        >  1:2137523/1‑36 (MQ=255)
             ttAATAACATCACGGTGATGACAAACAGCCTGCata                        <  1:3898872/36‑1 (MQ=255)
              taataaCATCACGGTGATGACCAACAGCCTGCatat                       <  1:4011003/36‑1 (MQ=255)
              taataaCACCACGTTAATAACCAACAGCCTGCatat                       <  1:2053948/36‑1 (MQ=16)
               aataaCATCACGGTGATGACAAACAGCCTGCATATc                      <  1:3541304/36‑1 (MQ=255)
               aataaCATCACGGCGACGACAAACAGCCTGCATATc                      >  1:2299705/1‑36 (MQ=21)
               aataaCAGCACGGTGATGACCAACAGCCTGCATATc                      <  1:2373828/36‑1 (MQ=255)
                ataaCATCACGGTGATGACCAACAGCCTGCATATCg                     <  1:3662949/36‑1 (MQ=255)
                 taaCATCGCGGTGAGCAAAAACAGCCTGCATATCGt                    <  1:4186779/36‑1 (MQ=12)
                 taaCATCACGGTGATGACCAACAGCCTGCATATCGt                    >  1:868722/1‑36 (MQ=255)
                 taaCATCACGGTGATGACCAACAGCCTGCATATCGt                    >  1:1466052/1‑36 (MQ=255)
                 taaCATCACGGTGATGACCAACAGCCTGCATATCGt                    >  1:1533055/1‑36 (MQ=255)
                 taaCATCACGGTGATGACCAACAGCCTGCATATCGt                    >  1:1349032/1‑36 (MQ=255)
                 taaCATCACGGTGATGACCAACAGCCTGCATATCGt                    >  1:1156064/1‑36 (MQ=255)
                 taaCATCACGGTGATCACAAACAGCCTGCATATCGt                    <  1:3249211/36‑1 (MQ=37)
                 aaaCATCAGGGTGATGACAAACAGCCTGCATATCGt                    <  1:3820315/35‑1 (MQ=25)
                  aaaaTCCCGGTGATAACAAAAAGCCTGCATATCGTc                   <  1:4514130/33‑1 (MQ=16)
                  aaCATCACGGTGATGACCAACAGCCTGCATATCGTc                   >  1:563511/1‑36 (MQ=255)
                  aaCATCACGGTGATGACCAACAGCCTGCATATCGTc                   >  1:5604075/1‑36 (MQ=255)
                    cATCACGGTGATGACCAACAGCCTGCATATCGTCaa                 >  1:711149/1‑36 (MQ=255)
                     aacACGGTGATGACAAACAGCCTGCATATCGTCAAt                <  1:1951937/34‑1 (MQ=37)
                     aTCACGGTGATGACCAACAGCCTGCATATCGTCAAt                >  1:4897389/1‑36 (MQ=255)
                     aTCACGGTGATGACCAACAGCCTGCATATCGTCAAt                >  1:179772/1‑36 (MQ=255)
                     aTCACGGTGATGACCAACAGCCTGAATATCGTCAAt                >  1:4123235/1‑36 (MQ=255)
                     aTCACGGTGATGACCAACAGCCTCCATATCGTCAAt                >  1:4986825/1‑36 (MQ=255)
                     aTCACGGTGATGACCAAAAGCCTGCATAGCGACAAt                >  1:4813205/1‑36 (MQ=25)
                     aTCACGGTGATGAACAACAGCCTGCATATCGTCAAt                >  1:5179176/1‑36 (MQ=255)
                      tCACGGTGATGACCAACAGCCTGCATATCGTCAATg               >  1:514557/1‑36 (MQ=255)
                      ccACGGGAATGACAAACAGCCTGCATATCGTCAATg               <  1:4977988/35‑1 (MQ=25)
                       cACGGTGATGACCAACAGCCTGCATATCGTCAATgc              >  1:1844508/1‑36 (MQ=255)
                       cACGGTGATGACCAACAGCCTGCATATAGTCAATgc              >  1:4852796/1‑36 (MQ=255)
                       cACGGTGATCCAAAACAGCCTGCATATCGTCAATgc              <  1:801991/36‑1 (MQ=19)
                         cGGTGATGACCAACAGCCTGCATATCGTCAATgcgc            >  1:5311954/1‑36 (MQ=255)
                         cGGAGATGACAAACAGCCTGCATATCGTCAATgcgc            <  1:5100363/36‑1 (MQ=37)
                          ggTGATGACCAACAGCCTGCATATCGTCAATgcgca           >  1:565577/1‑35 (MQ=255)
                          ggTGATGACCAACAGCCTGCATATCGTCAATGCGCt           >  1:5428421/1‑36 (MQ=255)
                          ggTGATGACCAACAGCCTGCATATCGTCAATGCGCt           >  1:3923398/1‑36 (MQ=255)
                          ggTGATGACCAACAGCCTGCATATCGTCAATGCGCt           >  1:2423106/1‑36 (MQ=255)
                          ggTGATGACCAACAGCCTGCATATCGTCAATGCGCt           >  1:1134863/1‑36 (MQ=255)
                          ggTGATGACCAACAGCCTGCATATCGTAAATGCGCt           >  1:1393792/1‑36 (MQ=255)
                          ggTGATGACAAACAGCCTGCATATCGTCAATGCGCt           <  1:5178154/36‑1 (MQ=255)
                          ggTGATGAACAACAGCCTGCATATCATCAAATCGCt           >  1:4344749/1‑36 (MQ=16)
                          ggTGATAACAAACAGCCTGCATATCGTCAATGCGCt           <  1:1713224/36‑1 (MQ=37)
                          ggTGAGGACAAACAGCCTGCATATCGTCAATGCGCt           <  1:4773413/36‑1 (MQ=37)
                          ggTAATGCCGAACAGCCTGCATATCGTCAATGCGCt           <  1:449125/36‑1 (MQ=21)
                           gTGATGACCAACAGCCTGCATATCGTCAATGCGCTa          >  1:2047574/1‑36 (MQ=255)
                           gTGATGACCAACAGCCTGCATATCGTCAATGCGCTa          >  1:845520/1‑36 (MQ=255)
                            tgatgaCCAACAGCCTGCATATCGTCAATGCGCTAt         >  1:520901/1‑36 (MQ=255)
                            tgatgaCCAACAGCCTGCATATCGTCAATGCGCTAt         >  1:2728510/1‑36 (MQ=255)
                            tgatgaCCAACAGCCTGCATATCGTCAATGCGCTAt         >  1:2793383/1‑36 (MQ=255)
                             taggaCCAACAGCATTCATATCGTCAATGCGCTATc        <  1:140112/33‑1 (MQ=21)
                             gttgaCAAACAGCCTGCATATCGTCAATGCGCTATc        <  1:5314227/34‑1 (MQ=37)
                             gatgaCCAACAGCCTGCATATCGTCAATGCGCTATc        >  1:893747/1‑36 (MQ=255)
                             gatgaCCAACAGCCTGCATATCGTCAATGCGCTATc        >  1:2363396/1‑36 (MQ=255)
                             aatgaCCAACAGCCTGCATATCGTCAATGCGCTATc        <  1:3169503/35‑1 (MQ=255)
                              atgaCCAACAGCCTGCATATCGTCAATGCGCTATcc       >  1:4180735/1‑36 (MQ=255)
                              atgaCCAACAGCCTGCATATCGTCAATGCGCTATcc       >  1:1886688/1‑36 (MQ=255)
                              atgaCCAACAGCCTGCATATCGTCAATGCGCTATcc       >  1:2851753/1‑36 (MQ=255)
                              atgaCCAACAGCCTGCATATCGTCAATGCGCTATcc       >  1:2806899/1‑36 (MQ=255)
                                gaCCAACAGCCTGCATATCGTCAATGCGCTATCCGa     >  1:2174372/1‑36 (MQ=255)
                                gaCCAACAGCCTGCATATCGTCAATGCGCTATCCGa     >  1:1877970/1‑36 (MQ=255)
                                gaCCAACAGCCTGCATATCGTAAATGCGCTATCCGa     >  1:4770110/1‑36 (MQ=255)
                                caCCAACAGCCTGACTATCGTCAATGCGCTATCCGa     <  1:4497045/35‑1 (MQ=37)
                                 aCCAACAGCCTGCATATCGTCAATGCGCTATCCGaa    >  1:4893786/1‑36 (MQ=255)
                                 aCCAACAGCCTGCATATCGTCAATGCGCTATCCGaa    >  1:3617567/1‑36 (MQ=255)
                                 aCCAACAGCCTGCATATCGTCAATGCGCTATCCGaa    >  1:2443108/1‑36 (MQ=255)
                                 aCCAACAGCCTGCATATCGTCAATGCGCTATCCGaa    >  1:3994416/1‑36 (MQ=255)
                                  ccACCAGCCTGCATATCGTCAATGCGCTATCCGAAc   <  1:1040662/36‑1 (MQ=255)
                                  ccAACAGCCTGCATATCGTCAATGCGCTATCCGAAc   >  1:2881011/1‑36 (MQ=255)
                                  ccAACAGCCTGCATATCGTCAATGCGCTATCCGAAc   <  1:2973665/36‑1 (MQ=255)
                                   cAACAGCCTGCATATCGTCAATGCGCTATCCGAACt  >  1:2294007/1‑36 (MQ=255)
                                   cAACAGCCTGCATATCGTCAATGCGCTATCCGAACt  <  1:252336/36‑1 (MQ=255)
                                   |                                   
CTGCTCTCGCGCTTTAATAACATCACGGTGATGACCAACAGCCTGCATATCGTCAATGCGCTATCCGAACT  >  REL606/2724231‑2724301

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG