Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,735,5000A→G13.9% 8.5 ‑15.6 0.0971117ascFfused cellobiose/arbutin/salicin‑specific PTS enzymes: IIB component/IC component
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (15/0):  ref base (24/58):  total (39/58)
Fisher's exact test strand distribution p-value = 1.62e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATACGGCGTGGCGATCCGC  >  REL606/2735465‑2735535
                                   |                                   
cgcTGGCTGCGGCGGCATCAGCCATTATGGCGGGGa                                     <  1:3854254/36‑1 (MQ=255)
 gctggctgCGGCGGCATCAGCCATTATGGCGGGGAt                                    <  1:3994968/36‑1 (MQ=255)
  ctggctgCGGCGGCATCAGCCATTATGGCGGGGGtt                                   >  1:4645913/1‑36 (MQ=255)
  ctggctgCGGCGGCATCAGCCATTATGGCGGGGAtt                                   <  1:1422639/36‑1 (MQ=255)
  ctggctgCGGCGGCATCAGCCATTATGGCGGGGAtt                                   <  1:6561902/36‑1 (MQ=255)
   tggctgCGGCGGCATCAGCCATTATGGCGGGGGttt                                  >  1:6972481/1‑36 (MQ=255)
   tggctgCGGCGGAATCAGCCGTTATGGCGGGGAttt                                  <  1:26112/36‑1 (MQ=37)
    ggctgCGGCGGCATCAGCCATTATGGCGGGGATTTc                                 <  1:280013/36‑1 (MQ=255)
    ggcggCGGCGGCATCAGCCATTATGGCGGGGATTTc                                 <  1:3841456/36‑1 (MQ=255)
     gctgCGGCGGCATCAGCCATTATGGCGGGGGttttc                                >  1:6932088/1‑34 (MQ=37)
     gctgCGGCGGCATCAGCCATTATGGCGGGGATTTcc                                <  1:4097395/36‑1 (MQ=255)
     gctgCGGCGGCATCAGCCATTATGGCGGGGATTTcc                                <  1:4953612/36‑1 (MQ=255)
     gctgCGGCGGCATCAGCCATTATGGCGGGGATTTcc                                <  1:4436665/36‑1 (MQ=255)
     gctgCGGCGGCATCAGCCATTATGGCGGGGATTTcc                                <  1:1617342/36‑1 (MQ=255)
     gctgCGGCGGCATCAGCCATTATGGCGGGGATTTcc                                <  1:3492313/36‑1 (MQ=255)
     gcggcggcgGCATCAGCTGTTATGGCGGGGATTTcc                                <  1:4744159/33‑1 (MQ=21)
      ctgCGGCGGCATCCGCCATTATGGCGGGGGttttcc                               >  1:161239/1‑33 (MQ=25)
        gcggcgGCATCAGCCATTATGGCGGGGGTTTTCGaa                             >  1:7937747/1‑36 (MQ=37)
        gcggcgGCATCAGCCATTATGGCGGGGATTTCCGaa                             <  1:1825656/36‑1 (MQ=255)
         gggcggcATCAGCCATTATGGCGGGGATTTCCGAAc                            <  1:7608883/35‑1 (MQ=255)
         cggcggCATCAGCCATTATGGCGGGGATTTCCGAAc                            <  1:454811/36‑1 (MQ=255)
         cggcggCATCAGCCATTATGGCGGGGATTTCCGAAc                            <  1:6756327/36‑1 (MQ=255)
         cggcggCACCAGCCATTATGGCGGGGATTTCCGAAc                            <  1:8269614/36‑1 (MQ=255)
          ggcggcATCAGCCATTATGGCGGGGGTNTCCGaaac                           >  1:8360727/1‑34 (MQ=37)
          ggcggcATCAGCCATTATGGCGGGGATTTCCGAAcc                           <  1:1620356/36‑1 (MQ=255)
          ggcggcATCAGCCATTATGGCGGGGATTTCCGAAcc                           <  1:2769911/36‑1 (MQ=255)
           gcgtcATCAGCCATTATGGCGGGGATTTCCGAACCg                          <  1:1389456/36‑1 (MQ=255)
           gcggcATCAGCCATTATGGCGGGGGTTTTCGAACCg                          >  1:2728874/1‑36 (MQ=37)
           gcggcATCAGCCATTATGGCGGGGATTTCCGAACCg                          <  1:240679/36‑1 (MQ=255)
             ggcATCAGCCATTATGGCGGGGGTTTTCTGACCGGc                        >  1:2426755/1‑36 (MQ=16)
             ggcATCAGCCATTATGGCGGGGGTTTTCGAAcccgg                        >  1:6924192/1‑33 (MQ=25)
             ggcATCAGCCATTATGGCGGGGGTTTTCGAACCGGc                        >  1:6071331/1‑36 (MQ=37)
             ggcATCAGCCATTATGGCGGGGATTTTCGAACCggg                        >  1:2053744/1‑35 (MQ=255)
             ggcATCAGCCATTATGGCGGGGATTTCCGCACCGGc                        >  1:7519439/1‑36 (MQ=255)
             ggcATCAGCCATTATGGCGGGGATTTCCGAACCGGc                        <  1:1190886/36‑1 (MQ=255)
             ggcATCAGCCATTATGGCGGGGATTTCCGAACCGGc                        <  1:1524597/36‑1 (MQ=255)
             ggcATCAGCCATTATGGCGGGGATTTCCGAACCGGc                        >  1:8427725/1‑36 (MQ=255)
               cATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt                      <  1:3639249/36‑1 (MQ=255)
               cATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt                      <  1:1412332/36‑1 (MQ=255)
               cATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt                      <  1:8075231/36‑1 (MQ=255)
                aTCAGCCATTATGGCGGGGGTTTTCGAACCCGCGtt                     >  1:6304714/1‑36 (MQ=25)
                aTCAGCCATTATGGCGGGGGTTTCCCAACCCGCggt                     >  1:2693246/1‑34 (MQ=21)
                aTCAGCCATTATGGCGGGGATTTCCGAACCGGCGtt                     <  1:1932971/36‑1 (MQ=255)
                 tCAGCCATTATGGCGGGGATTTCCGAAACGGCGTta                    >  1:8308553/1‑36 (MQ=255)
                  cAGCCATTATGGCGGGGGTTTTCGCACCGGCGTtat                   >  1:1972302/1‑36 (MQ=25)
                  cAGCCATTATGGCGGGGATTTCCGAACCGGCGTtat                   <  1:1337459/36‑1 (MQ=255)
                  cAGCCATTATGGCGGGGATTTCCGAACCGGCGTtat                   <  1:6636527/36‑1 (MQ=255)
                  cAGCCATTATGGCGGGGATTTCCGAACCGGCGTtat                   <  1:4756686/36‑1 (MQ=255)
                   aGCCATTATGGCGGGGATTTTCCAAACGTCGGtata                  >  1:6130782/1‑36 (MQ=9)
                   aGCCATTATGGCGGGGATTTCCGAACCGGCGTttta                  >  1:7087091/1‑36 (MQ=255)
                    gCCATTATGGCGGGGATTTCCGAACCGGCGTTATAc                 <  1:5892332/36‑1 (MQ=255)
                    gCCATTATGGCGGGGATTTCCGAACCGGCGTTATAc                 <  1:1353971/36‑1 (MQ=255)
                    gCCATTATGGCGGGGATTTCCGAACCGGCGTTATAc                 <  1:7607930/36‑1 (MQ=255)
                     ccATTATGGCGGGGGTTTCCCAACCGGCGTTATACg                >  1:4386133/1‑36 (MQ=37)
                     ccATTATGGCGGGGATTTTCGAACCGGCGTTATACg                >  1:8045653/1‑36 (MQ=255)
                     ccATTATGGCGGGGATTTCCGCACCGGCGTTATACg                >  1:4207349/1‑36 (MQ=255)
                     ccATTATGGCGGGGATTTCCGAACCGGCGTTATACg                <  1:7691772/36‑1 (MQ=255)
                     ccATTATGGCGGGGATTTCCGAACCGGCGTTATACg                <  1:4126886/36‑1 (MQ=255)
                      cATTATGGCGGGGGTTTTCGAAACGGCGTTATAAgg               >  1:2307975/1‑36 (MQ=19)
                      cATTATGGCGGGGATTTCCGAACCGGCGTTATACgg               <  1:6199744/36‑1 (MQ=255)
                      cATTATGGCGGGGATTTCCGAACCGGCGTTATACgg               >  1:1912904/1‑36 (MQ=255)
                       aTTATGGCGGGGATTTCCGAACCGGCGTTATACGGc              <  1:4091898/36‑1 (MQ=255)
                       aTTATGGCGGGGATTTCCGAACCGGCGTTATACGGc              <  1:7015184/36‑1 (MQ=255)
                         tATTGCGGGGATTTCCGAACCGGCGTTATACGGCGt            <  1:7574630/36‑1 (MQ=255)
                         tATGGCGGGGATTTCCGACCCGGCGTTATACGGCGt            <  1:1701060/36‑1 (MQ=255)
                         tATGGCGGGGATTTCCGAACCGGCGTTATACGGCGt            <  1:8164406/36‑1 (MQ=255)
                         tATGGCGGGGATTTCCGAACCGGCGTTATACGGCGt            <  1:841012/36‑1 (MQ=255)
                         tATGGCGGGGATTTCCGAACCGGCGTTATACGGCGt            <  1:2971915/36‑1 (MQ=255)
                         tATGGCGGGGATTTCCGAACCGGCGTTATACGGCGt            <  1:2439200/36‑1 (MQ=255)
                         tATGGCGGGGATTTCCGAACCGGCGTTATACGGCGt            <  1:2442582/36‑1 (MQ=255)
                          aTGGCGGGGATTTCCGAACCGGCGTTATACGGCGTg           <  1:7315468/36‑1 (MQ=255)
                           tGGCGGGGATTTCCGAACCGGCGTTATACGGCGTgg          >  1:6106698/1‑36 (MQ=255)
                           tGGCGGGGATTTCCGAACCGGCGTTATACGGCGTgg          >  1:631537/1‑36 (MQ=255)
                            ggCGGGGATTTCCGAACCGGCGTTATACGGCGTGGc         >  1:6410010/1‑36 (MQ=255)
                            ggCGGGGATTTCCGAACCGGCGTTATACGGCGTGGc         >  1:3882467/1‑36 (MQ=255)
                             gCGGGGATTTCCGAACCGGCGTTATACGGCGTGGCg        <  1:6802005/36‑1 (MQ=255)
                             gCGGGGATTTCCGAACCGGCGTTATACGGCGTGGCg        >  1:6664886/1‑36 (MQ=255)
                             gCGGGGATTTCCGAACCGGCGTTATACGGCGTGGCg        >  1:785758/1‑36 (MQ=255)
                             gCGGGGATTTCCGAACCGGCGTTATACGGCGTGGCg        <  1:5328226/36‑1 (MQ=255)
                             gCGGGGATTGCCGAACCGGCGGTATACGGCGTggag        >  1:614477/1‑34 (MQ=25)
                              cGGGGATTTCCGAACCGGCGTTATACGGCGTGGCGa       <  1:4366864/36‑1 (MQ=255)
                              cGGGGATTTCCGAACCGGCGTTATACGGCGTGGCGa       <  1:7140515/36‑1 (MQ=255)
                              cGGGGATTTCCGAACCGGCGTTATACGGCGTGGCGa       <  1:5226966/36‑1 (MQ=255)
                               ggggATTTCCGAACCGGCGTTATACGGCGTGGCGAt      >  1:7935644/1‑36 (MQ=255)
                               ggggATTTCCGAACCGGCGTTATACGGCGTGGCGAt      <  1:923429/36‑1 (MQ=255)
                               ggggATTTCCGAACCGGCGTTATACGGCGTGGCGAt      <  1:7075741/36‑1 (MQ=255)
                               ggggATTTCCGAACCGGCGTTATACGGCGTGGCGAt      >  1:3418686/1‑36 (MQ=255)
                               ggggATTTCCGAACCGGCGTTATACGGCGTGGCGAt      >  1:266566/1‑36 (MQ=255)
                               ggggATTTCCGAACCGGCGTTATACGGCGTGGCGAt      >  1:1602894/1‑36 (MQ=255)
                               ggggATTTCCGAACCGGCGTTATACGGCGTGGCGAt      >  1:1023529/1‑36 (MQ=255)
                                gggATTTCCGAACCGGCGTTATACGGCGTGGCGATc     <  1:645356/36‑1 (MQ=255)
                                gggATTTCCGAACCGGCGTTATACGGCGTGGCGATc     <  1:6088132/36‑1 (MQ=255)
                                gggATTTCCGAACCGGCGTTATACGGCGTGGCGATc     >  1:5730525/1‑36 (MQ=255)
                                gggATTTCCGAACCGGCGTTATACGGCGTGGCGATc     <  1:5568215/36‑1 (MQ=255)
                                gggATTTCCGAACCGGCGTTATACGGCGTGGCGATc     <  1:49175/36‑1 (MQ=255)
                                gggATTTCCGAACCGGCGTTATACGGCGTGGCGATc     >  1:2040756/1‑36 (MQ=255)
                                gggATTTCCGAACCGGCGTTATACGGCGTGGCGATc     <  1:1458988/36‑1 (MQ=255)
                                gggATTTCCGAACAGGCGTTATACGGCGTTGCGCTc     >  1:487404/1‑36 (MQ=25)
                                   aTTTCCGAACCGGCGTTATACGGCGTGGCGATCCGc  <  1:5478923/36‑1 (MQ=255)
                                   |                                   
CGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATACGGCGTGGCGATCCGC  >  REL606/2735465‑2735535

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 19 ≤  ATCG  < 28 ≤  ATCG  < 37 ≤  ATCG