Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,304,2750C→A16.4% 5.9 ‑17.0 0.092400sspBClpXP protease specificity‑enhancing factor
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (0/18):  ref base (50/23):  total (50/42)
Fisher's exact test strand distribution p-value = 4.28e-08
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

ATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATG  >  REL606/3304240‑3304310
                                   |                                   
aTGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                     <  1:395794/36‑1 (MQ=255)
aTGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                     <  1:4064381/36‑1 (MQ=255)
aTGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                     <  1:2721306/36‑1 (MQ=255)
aTGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                     >  1:7414375/1‑36 (MQ=255)
aTGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                     <  1:7435873/36‑1 (MQ=255)
aTGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                     >  1:7478239/1‑36 (MQ=255)
aTGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                     <  1:5342932/36‑1 (MQ=255)
 tGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                    <  1:6534376/36‑1 (MQ=255)
 tGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCa                                    <  1:350849/36‑1 (MQ=255)
  gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa                                   >  1:7885704/1‑36 (MQ=255)
  gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa                                   <  1:135077/36‑1 (MQ=255)
  gAGTGTCATCATCGTGATCTGGCTTGTCGCCATCaa                                   >  1:2636077/1‑36 (MQ=255)
   aGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                                  >  1:746019/1‑36 (MQ=255)
   aGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                                  >  1:2149307/1‑36 (MQ=255)
   aGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                                  >  1:4336986/1‑36 (MQ=255)
   aGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                                  >  1:2547961/1‑36 (MQ=255)
   aGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                                  <  1:6027099/36‑1 (MQ=255)
   aGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                                  <  1:6389772/36‑1 (MQ=255)
    gtgtCCTCATCGTGATCTGGCTTGTCGCCATCAATa                                 <  1:7540031/36‑1 (MQ=255)
    gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa                                 >  1:3030663/1‑36 (MQ=255)
    gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa                                 <  1:407909/36‑1 (MQ=255)
    gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa                                 <  1:2895319/36‑1 (MQ=255)
    gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa                                 >  1:6415661/1‑36 (MQ=255)
    gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa                                 <  1:5856313/36‑1 (MQ=255)
    gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa                                 >  1:2252107/1‑36 (MQ=255)
    gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa                                 <  1:7699081/36‑1 (MQ=255)
    gtgtCATCATCGTGATCTGGCTTGTCGCCATCAATa                                 <  1:5711536/36‑1 (MQ=255)
     tgtCATCATCGTGATCTGGCTTGTCGCCATCAATaa                                <  1:3135766/36‑1 (MQ=255)
     tgtCATCATCGTGATCTGGCTTGTCGCCATCAATaa                                >  1:3936755/1‑36 (MQ=255)
     tgtCATCATCGTGATCTGGCTTGTCGCCATCAATaa                                <  1:1707258/36‑1 (MQ=255)
      gtCATCATCGTGATCTGGCTTGTCGCCATCAATAAc                               <  1:5913987/36‑1 (MQ=255)
      gtCATCATCGTGATCTGGCTTGTCGCCATCAATAAc                               >  1:216557/1‑36 (MQ=255)
       tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc                              >  1:3059715/1‑36 (MQ=255)
       tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc                              >  1:5672010/1‑36 (MQ=255)
       tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc                              <  1:5593206/36‑1 (MQ=255)
       tcatcaTCGTGATCTGGCTTGTCGCCATCAATAAcc                              >  1:7291229/1‑36 (MQ=255)
        catcatCGTGATCTGGCTTGTCGCCATCAATAACCg                             <  1:7006990/36‑1 (MQ=255)
        catcatCGTGATCTGGCTTGTCGCCATCAATAACCg                             <  1:6430463/36‑1 (MQ=255)
        catcatCGTGATCTGGCTTGTCGCCATCAATAACCg                             <  1:1338320/36‑1 (MQ=255)
          tcatcGTGATCTGGCTTGTCGCCATCAATAACCGAc                           >  1:4954796/1‑36 (MQ=255)
          tcatcGTGATCTGGCTTGTCGCCATCAATAACCGAc                           <  1:3692636/36‑1 (MQ=255)
           catcGTGATCTGGCTTGTCGCCATCAATAACCGACa                          >  1:7235860/1‑36 (MQ=255)
           catcGTGATCTGGCTTGTCGCCATCAATAACCGACa                          >  1:5306040/1‑36 (MQ=255)
            atcGTGATCTGGCTTGTCGCCATCAATAACCGACAt                         >  1:3494842/1‑36 (MQ=255)
            atcGTGATCTGGCTTGTCGCCATCAATAACCGACAt                         >  1:3015127/1‑36 (MQ=255)
            atcGTGATCTGGCTTGTCGCCATAAATAACCGACAt                         <  1:2549628/36‑1 (MQ=255)
             tcGTGATCTGGCTTGTCGCCATCAATAACCGACATa                        >  1:4122086/1‑36 (MQ=255)
             tcGTGATCTGGCTTGTCGCCATCAATAACCGACATa                        >  1:2823853/1‑36 (MQ=255)
             tcGTGATCTGGCTTGTCGCCATCAATAACCGACATa                        >  1:4796643/1‑36 (MQ=255)
              cGTGATCTGGCTTGTCGCCATCAATAACCGACATaa                       >  1:7228444/1‑36 (MQ=255)
              cGTGATCTGGCTTGTCGCCATCAATAACCGACATaa                       >  1:4478399/1‑36 (MQ=255)
               gTGATCTGGCTTGTCGCCATCAATAACCGACATAAc                      >  1:4945108/1‑36 (MQ=255)
               gTGATCTGGCTTGTCGCCATCAATAACCGACATAAc                      >  1:911624/1‑36 (MQ=255)
                tGATCTGGCTTGTCGCCATCAATAACCGACATAACg                     >  1:730915/1‑36 (MQ=255)
                tGAGCTGGCTTGTCGCAATAAATAACCGACATAACg                     <  1:199519/36‑1 (MQ=25)
                 gATCTGGCTTGTCGCCATCAATAACCGACATAACgg                    >  1:91334/1‑36 (MQ=255)
                 gATCTGGCTTGTCGCCATCAATAACCGACATAACgg                    >  1:30515/1‑36 (MQ=255)
                 gATCTGGCTTGTCGCCATAAATAACCGACATAACgg                    <  1:1424257/36‑1 (MQ=255)
                  aTCTGGCTTGTCGCCATCAATAACCGACATAACGGt                   >  1:1271938/1‑36 (MQ=255)
                  aTCTGGCTTGTCGCCATCAATAACCGACATAACGGt                   >  1:1284563/1‑36 (MQ=255)
                  aTCTGGCTTGTCGCCATCAATAACCGACATAACGGt                   >  1:1741161/1‑36 (MQ=255)
                  aTCTGGCTTGTCGCCATCAATAACCGACATAACGGt                   >  1:217338/1‑36 (MQ=255)
                  aTCTGGCTTGTCGCAATAAATAACCGACATAACGGt                   <  1:3090/36‑1 (MQ=37)
                   tCTGGCTTTTCACCATAAATAACCGACATAACGGtt                  <  1:3009931/36‑1 (MQ=25)
                   tCTGGCTTGTCGCCATCAATAACCGACATAACGGtt                  >  1:3222644/1‑36 (MQ=255)
                    cTGGCTTGTCGCCATAAATAACCGACATAACGGttt                 <  1:7255935/36‑1 (MQ=255)
                    cTGGCTTGTCCCAATAAATAACCGACATAACGGttt                 <  1:2079243/36‑1 (MQ=25)
                    cTGGCTTGGCTCCATAAATACCCGACATAACGGttt                 <  1:6227448/36‑1 (MQ=19)
                     ttgCTTGCCTCAATCAATAACCGACATAACGGTTTc                <  1:4946301/34‑1 (MQ=25)
                     tGGCTTGTCGCCATAAATAACCGACATAACGGTTTc                <  1:4184535/36‑1 (MQ=255)
                     tGGCTGGCCGCCAACAATAACCGACATAACGGTTTc                <  1:3272834/36‑1 (MQ=21)
                      ggCTTGTCGCCATCAATAACCGACATAACGGTTTCg               >  1:1695675/1‑36 (MQ=255)
                       gCTTGTCGCCATCAATAACCGACATAACGGTTTCGt              >  1:5104467/1‑36 (MQ=255)
                       gCTTGTCGCCATCAATAACCGACATAACGGTTTCGt              >  1:2018078/1‑36 (MQ=255)
                        ctcgTCCCCATCAAAAACCGACATAACGGTTTCGtt             <  1:531102/33‑1 (MQ=25)
                        cTTGTCGCCATCAATAACCGACATAACGGTTTCGtt             >  1:3325210/1‑36 (MQ=255)
                        cTTGTCGCCATCAATAACCGACATAACGGTTTCGtt             >  1:7542759/1‑36 (MQ=255)
                        cTTGTCGCCATCAATAACCGACATAACGGTTTCGtt             >  1:7311384/1‑36 (MQ=255)
                        cTTGTCGCAATAAATAACCGACATAACGGTTTCGtt             <  1:4244552/36‑1 (MQ=37)
                        cTTGTCCCAATAAATAACCGACATAACGGTTTCGtt             <  1:7569599/36‑1 (MQ=25)
                         ttttCGCCAAAAAAAACCGACATAACGGTTTCGTTg            <  1:4640114/33‑1 (MQ=16)
                         ttGTCGCCATCAATAACCGACATAACGGTTTCGTTg            >  1:6095399/1‑36 (MQ=255)
                         ttGTCGCAATAAATAACCGACATAACGGTTTCGTTg            <  1:3465149/36‑1 (MQ=37)
                         ctGTCGCCATCAATAACCGACATAACGGTTTCGttt            >  1:2769352/2‑35 (MQ=255)
                          tgccGCAATGAATAACCGACATAACGGTTTCGTTGt           <  1:2539344/33‑1 (MQ=37)
                          tGTCGCCATACAAAACCGACATAACGGTTTCGTTGt           <  1:7020817/36‑1 (MQ=25)
                           gTCGCCATCAATAACCGACATAACGGTTTCGTTGTc          >  1:2018713/1‑36 (MQ=255)
                           gTCGCAATAAATAACCGACATAACGGTTTCGTTGTc          <  1:3277653/36‑1 (MQ=37)
                            tCGCCATCAATAACCGACATAACGGTTTCGTTGTCt         >  1:7833325/1‑36 (MQ=255)
                            tCGCCATCAATAACCGACATAACGGTTTCGTTGTCt         >  1:7908070/1‑36 (MQ=255)
                             cGCCATAAATAACCGACATAACGGTTTCGTTGTCTg        <  1:6353298/36‑1 (MQ=255)
                             cGCCATAAATAACCGACATAACGGTTTCGTTGTCTg        <  1:6108165/36‑1 (MQ=255)
                              gcaaTCAATAACCGACATAACGGTTTCGTTGTCTGc       <  1:3364523/33‑1 (MQ=255)
                              gCCATCAATAACCGACATAACGGTTTCGTTGTCTGc       >  1:5188149/1‑36 (MQ=255)
                              gCCATCAATAACCGACATAACGGTTTCGTTGTCTGc       >  1:6237627/1‑36 (MQ=255)
                              gCCATCAATAACCGACATAACGGTTTCGTTGTCTGc       >  1:4853726/1‑36 (MQ=255)
                               ccATCAATAACCGACATAACGGTTTCGTTGTCTGcc      >  1:3551894/1‑36 (MQ=255)
                               ccATAAAAAACCGACATAACGGTTTCGTTGTCTGcc      <  1:1939806/36‑1 (MQ=38)
                                  acAATACCCGACATAACGGTTTCGTTGTCTGCCGAt   <  1:5804564/35‑1 (MQ=255)
                                   cAATAACCGACATAACGGTTTCGTTGTCTGCCGATg  >  1:4808604/1‑36 (MQ=255)
                                   |                                   
ATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGCCGATG  >  REL606/3304240‑3304310

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG