Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,304,2750C→A13.2% 4.2 ‑9.5 0.086400sspBClpXP protease specificity‑enhancing factor
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (0/15):  ref base (38/33):  total (38/48)
Fisher's exact test strand distribution p-value = 7.22e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATG  >  REL606/3304240‑3304310
                                   |                                   
aTGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                     <  1:7400594/36‑1 (MQ=255)
 tGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                    <  1:6660308/36‑1 (MQ=255)
 tGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                    <  1:4598952/36‑1 (MQ=255)
 tGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                    <  1:4330236/36‑1 (MQ=255)
 tGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                    <  1:4321748/36‑1 (MQ=255)
 tGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                    <  1:1808831/36‑1 (MQ=255)
  gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa                                   >  1:7561913/1‑36 (MQ=255)
  gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa                                   <  1:2591158/36‑1 (MQ=255)
  gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa                                   <  1:3564865/36‑1 (MQ=255)
  gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa                                   <  1:209344/36‑1 (MQ=255)
  gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa                                   <  1:2353235/36‑1 (MQ=255)
   aGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                                  <  1:4216631/36‑1 (MQ=255)
   aGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                                  <  1:3253037/36‑1 (MQ=255)
    gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa                                 <  1:3587971/36‑1 (MQ=255)
    gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa                                 >  1:3808376/1‑36 (MQ=255)
    gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa                                 <  1:1943920/36‑1 (MQ=255)
     tgtCATCATCGTGATCTGGCTTGTCGCCATCAATaa                                >  1:2967477/1‑36 (MQ=255)
     tgtCATCATCGTGATCTGGCTTGTCGCCATCAATaa                                <  1:862914/36‑1 (MQ=255)
      gtCATCATCGTGATCTGGCTTGTCGCCATCAATAAc                               >  1:5637189/1‑36 (MQ=255)
      gtCATCATCGTGATCTGGCTTGTCGCCATCAATAAc                               <  1:4066736/36‑1 (MQ=255)
       tcatcaTCGTGATCTGGCTTGTCGCCATCCATAAcc                              <  1:6643762/36‑1 (MQ=255)
       tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc                              <  1:6809978/36‑1 (MQ=255)
       tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc                              <  1:3415714/36‑1 (MQ=255)
       tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc                              <  1:3463209/36‑1 (MQ=255)
       tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc                              <  1:2484326/36‑1 (MQ=255)
       tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc                              <  1:3651558/36‑1 (MQ=255)
       tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc                              <  1:4890253/36‑1 (MQ=255)
        catcatCGTGATCTGGCTTGTCGCCATCAATAACCg                             >  1:1732722/1‑36 (MQ=255)
        catcatCGTGATCTGGCTTGTCGCCATCAATAACCg                             >  1:7382900/1‑36 (MQ=255)
          tcatcGTGATCTGGCTTGTCGCCATCAATAACCGAc                           <  1:761193/36‑1 (MQ=255)
          tcatcGTGATCTGGCTTGTCGCCATCAATAACCGAc                           <  1:1773492/36‑1 (MQ=255)
          tcatcGTGATCTGGCTTGTCGCCATCAATAACCGAc                           <  1:5054900/36‑1 (MQ=255)
           catcGTGATCTGGCTTGTCGCCATCAATAACCGACa                          >  1:4344950/1‑36 (MQ=255)
            atcGTGATCTGGCTTTTCGCCATCAATAACCGACAt                         <  1:6856034/36‑1 (MQ=255)
            atcGTGATCTGGCTTGTCGCCATCAATAACCGACAt                         >  1:4866153/1‑36 (MQ=255)
            atcGTGATCTGGCTTGTCGCCATCAATAACCGACAt                         >  1:4639302/1‑36 (MQ=255)
            atcGTGATCTGGCTTGTCGCCATCAATAACCGACAt                         <  1:1897014/36‑1 (MQ=255)
            atcGTGATCTGGCTTGTCGCCATCAATAACCGACAt                         >  1:5299033/1‑36 (MQ=255)
            atcGTGATCTGGCTTGTCGCCATCAATAACCGACAt                         >  1:5301662/1‑36 (MQ=255)
            agcGTGATCTGGCTTGTCGCCATCAATAACCGACAt                         <  1:2390317/34‑1 (MQ=255)
             tcGTGATCTGGCTTGTCGCCATCAATAACCGACATa                        >  1:4422839/1‑36 (MQ=255)
              cGTGATCTGGCTTGTCGCCATCAATAACCGACATaa                       >  1:7000089/1‑36 (MQ=255)
               gTGATCTGGCTTGTCGCCATCAATAACCGACATAAc                      >  1:3827992/1‑36 (MQ=255)
               gTGATCTGGCTTGTCGCCAAAAAAAACCGACATAAc                      <  1:7124110/36‑1 (MQ=21)
                tGATCTGGTTTGTCACCATAAATAACCGACATAACg                     <  1:4773768/36‑1 (MQ=25)
                tGATCTGGCTTGTCGCCATCAATAACCGACATAACg                     >  1:3678126/1‑36 (MQ=255)
                 gATCTGTCTTCTCGCAAAAAATAACCGACATAACgg                    <  1:1336613/36‑1 (MQ=12)
                 gATCTGGCTTGTCGCCATCAATAACCGACATAACgg                    >  1:842225/1‑36 (MQ=255)
                 gATCTGGCTTGTCGCCATCAATAACCGACATAACgg                    >  1:7506715/1‑36 (MQ=255)
                 gATCTGGCTTGTCGCCATAAATAACCGACATAACgg                    <  1:59899/36‑1 (MQ=255)
                  aTCTGGCTTGTCGCCATAAAAAACCGACATAACGGt                   <  1:1822065/36‑1 (MQ=37)
                   tCTGGCTTGTCGCCATCAATAACCGACATAACGGtt                  >  1:1957554/1‑36 (MQ=255)
                    cgggCTTGTCGCAATAAATAACCGACATAACGGttt                 <  1:5702652/34‑1 (MQ=37)
                    cTGGCTTGTCGCCATCAATAACCGACATAACGGttt                 >  1:485364/1‑36 (MQ=255)
                    cTGGCTTGTCGCCATCAATAACCGACATAACGGttt                 >  1:2981015/1‑36 (MQ=255)
                    cTGGCTTGTCGCCATCAATAACCGACATAACGGttt                 >  1:6005296/1‑36 (MQ=255)
                    cTGGCTTGTCGCCATCAATAACCGACATAACGGttt                 >  1:7551454/1‑36 (MQ=255)
                    cTGGCTTGTCGCCATAAATAACCGACATAACGGttt                 <  1:6117806/36‑1 (MQ=255)
                    cTGGCTTGTCGCAATAAATAACCGACATAACGGttt                 <  1:587298/36‑1 (MQ=37)
                     tGGCTTGTCGCCATCAATAACCGACATAACGGTTTc                <  1:7587774/36‑1 (MQ=255)
                     tGGCTTGTCGCCATCAATAACCGACATAACGGTTTc                <  1:6628959/36‑1 (MQ=255)
                     tGGCTTGTCGCAATCAATAACCGACATAACGGTTTc                <  1:4883891/36‑1 (MQ=255)
                      ggCTTGTCGCCATCAATAACCGACATAACGGTTTCg               >  1:525499/1‑36 (MQ=255)
                       gCTTGTCGCCATCAAAAACCGACATAACGGTTTCGt              <  1:2551649/36‑1 (MQ=255)
                        cTTGTCGCCATCAATAACCGACATAACGGTTTCGtt             >  1:4439245/1‑36 (MQ=255)
                        cTTGTCGCCATAAATAACCGACATAACGGTTTCGtt             <  1:3313446/36‑1 (MQ=255)
                        cTTGTCGCCATAAATAACCGACATAACGGTTTCGtt             <  1:7812143/36‑1 (MQ=255)
                        cTTGTCACCATAAATAACCGACATAACGGTTTCGtt             <  1:3714796/36‑1 (MQ=37)
                        cTTGCCGCCATAAATAACCGACATAACGGTTTCGtt             <  1:2543404/36‑1 (MQ=37)
                        cTTGCCGCCATAAAAACCCGACATAACGGTTTCGtt             <  1:6368225/36‑1 (MQ=19)
                         ttGTCGCCATCAATAACCGACATAACGGTTTCGTTg            >  1:2651840/1‑36 (MQ=255)
                         ttGTCGCCATCAATAACCGACATAACGGTTTCGTTg            >  1:6995327/1‑36 (MQ=255)
                          tGTCGCCATCAATAACCGACATAACGGTTTCGTTGt           >  1:2896478/1‑36 (MQ=255)
                           gTCGCCATCAATAACCGACATAACGGTTTCGTTGTc          >  1:6637869/1‑36 (MQ=255)
                            tCGCCATCAATAACCGACATAACGGTTTCGTTGTCt         >  1:1549311/1‑36 (MQ=255)
                            tCGCCATCAATAACCGACATAACGGTTTCGTTGTCt         >  1:425605/1‑36 (MQ=255)
                            tCGCCATAAATAACCGACATAACGGTTTCGTTGTCt         <  1:5470358/36‑1 (MQ=255)
                              gCCATCAATAACCGACATAACGGTTTCGTTGTCTGc       >  1:6963054/1‑36 (MQ=255)
                              gCCATCAATAACCGACATAACGGTTTCGTTGTCTGc       >  1:6562222/1‑36 (MQ=255)
                              aaaaTCAATCACCGACATAACGGTTTCGTTGTCTGc       <  1:403527/33‑1 (MQ=37)
                                cATCAATAACCGACATAACGGTTTCGTTGTCTGCCg     >  1:6021335/1‑36 (MQ=255)
                                cATCAATAACCGACATAACGGTTTCGTTGTCTGCCg     >  1:944769/1‑36 (MQ=255)
                                cATAAATAACCGACATAACGGTTTCGTTGTCTGCCg     <  1:5430382/36‑1 (MQ=255)
                                 aacAATAACCGACATAACGGTTTCGTTGTCTGCCGa    <  1:1953447/34‑1 (MQ=255)
                                 aTCAATAACCGACATAACGGTTTCGTTGTCTGCCGa    >  1:739814/1‑36 (MQ=255)
                                 aTCAATAACCGACATAACGGTTTCGTTGTCTGCCGa    >  1:5710992/1‑36 (MQ=255)
                                 aTCAATAACCGACATAACGGTTTCGTTGTCTGCCGa    >  1:1746781/1‑36 (MQ=255)
                                  tCAATAACCGACATAACGGTTTCGTTGTCTGCCGAt   >  1:6816129/1‑36 (MQ=255)
                                   cAATAACCGACATAACGGTTTCGTTGTCTGCCGATg  >  1:5438916/1‑36 (MQ=255)
                                   |                                   
ATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATG  >  REL606/3304240‑3304310

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG