Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,735,4291.→G 2.8% 0.3 ‑2.1 ‑0.071coding (1046/1458 nt)ascFfused cellobiose/arbutin/salicin‑specific PTS enzymes: IIB component/IC component
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (23/46):  total (25/46)
Fisher's exact test strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

TAACCTGTCGCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGAAAACGAAAAACCCGGAACTGCGCCAGACG  >  REL606/2735394‑2735463
                                     |                                  
tacccTGTCGCTGGGCGGTTCATCACTGGCGGT.GGc                                     <  1:4016504/33‑1 (MQ=255)
tAACCTGTCGCTGGGCGGTTCATCACTGGCGGT.GGc                                     <  1:7565284/36‑1 (MQ=255)
tAACCTGTCGCTGGGCGGTTCATCACTGGCGGT.GGc                                     <  1:2476536/36‑1 (MQ=255)
  aCCTGTCGCTGGGCGGTTCATCCCTGGCGGT.GGC.Gt                                  <  1:2269319/36‑1 (MQ=255)
  aCCTGTCGCTGGGCGGTTCATCACTGGCGGT.GGG.Gt                                  >  1:5959016/1‑36 (MQ=255)
  aCCTGTCGCTGGGCGGTTCATCACTGGCGGT.GGC.Gt                                  <  1:7078047/36‑1 (MQ=255)
  aCCTGTCGCTGGGCGGTTCATCACTGGCGGT.GGC.Gt                                  <  1:6640695/36‑1 (MQ=255)
  aCCTGTCGCTGGGCGGTTCATCACTGGCGGT.GGC.Gt                                  <  1:6307208/36‑1 (MQ=255)
  aCCTGTCGCTGGGCGGTTCATCACTGGCGGT.GGC.Gt                                  <  1:6224989/36‑1 (MQ=255)
  aCCTGTCGCTGGGCGGTGCATCGCTGGCGTT.GGC.Gt                                  <  1:7745481/36‑1 (MQ=21)
     tGTCGCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGa                               <  1:1448922/36‑1 (MQ=255)
      gTCGCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGaa                              <  1:6194967/36‑1 (MQ=255)
      gTCGCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGaa                              <  1:5665751/36‑1 (MQ=255)
      gTCGCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGaa                              <  1:5226542/36‑1 (MQ=255)
       gcGCTGGTGGGTTCATCACTGGCGGT.GGC.GTGGaaa                             <  1:4638952/35‑1 (MQ=25)
        cGCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGaaaa                            <  1:1040080/36‑1 (MQ=255)
        cGCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGaaaa                            <  1:1205784/36‑1 (MQ=255)
        cGCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGaaaa                            <  1:149341/36‑1 (MQ=255)
         gCTGGGCGGTTCATCACTGGCGGTGGGG.GTGGaaac                            >  1:5199332/1‑35 (MQ=21)
         gCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGAAAAc                           <  1:2195952/36‑1 (MQ=255)
         gCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGAAAAc                           <  1:90594/36‑1 (MQ=255)
         gCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGAAAAc                           <  1:2971022/36‑1 (MQ=255)
         gCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGAAAAc                           <  1:3049433/36‑1 (MQ=255)
         gCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGAAAAc                           <  1:6056211/36‑1 (MQ=255)
         gCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGAAAAc                           <  1:3646695/36‑1 (MQ=255)
         gCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGAAAAc                           <  1:3491215/36‑1 (MQ=255)
         gCTGGGCGGTTCATCACTGGCGGG.GGC.GGGGgaaaa                           >  1:6059195/1‑35 (MQ=21)
          cTGGGCGGTTCATCACTGGCGGT.GGG.GTGGgaaaag                          >  1:5485528/1‑34 (MQ=25)
          cTGGGCGGTTCATCACTGGCGGT.GGG.GGGGGAAACg                          >  1:895604/1‑36 (MQ=25)
           tGGGCGGTTCATCACTGGCGGT.GGG.GGGGGAAAcac                         >  1:404958/1‑34 (MQ=21)
            gggCGGTTCATCACTGGCGGT.GGCGGTGGaaaagaa                         >  1:8442230/1‑33 (MQ=25)
               cGGTTCATCACTGGCGGT.GGG.GTGGGAAAAGaaaaa                     >  1:7134169/1‑36 (MQ=25)
                ggTTCATCACTGGCGGT.GGG.GTGGTGAAAGaaaaaa                    >  1:4611389/1‑35 (MQ=16)
                ggTTCATCACTGGCGGT.GGC.GTGGAAAACGAAAAAc                    <  1:3410048/36‑1 (MQ=255)
                ggTTCATCACTGGCGGT.GGC.GTGGAAAACGAAAAAc                    <  1:444645/36‑1 (MQ=255)
                ggTTCATCACTGGCGGT.GGC.GTGGAAAACGAAAAAc                    <  1:1346977/36‑1 (MQ=255)
                ggTTCATCACTGGCGGT.GGC.GTGGAAAACGAAAAAc                    <  1:3364981/36‑1 (MQ=255)
                ggTTCATCACTGGCGGG.GGCGGGGGAAAACAaaaaa                     >  1:3085585/1‑36 (MQ=12)
                 gTTCATCACTGGCGGT.GGC.GTGGAAAACGAAAAAcc                   <  1:981825/36‑1 (MQ=255)
                 gTTCATCACTGGCGGC.GGG.GTGGGGAAAGAAAAAcc                   >  1:6935813/1‑36 (MQ=12)
                  ttCATCACTGGCGGT.GGC.GTTGAAAACGAAAAAccc                  <  1:1799523/36‑1 (MQ=255)
                  ttCATCACTGGCGGT.GGC.GTGGAAAACGAAAAAccc                  <  1:7109726/36‑1 (MQ=255)
                   tcatcaCTGGCGGT.GGC.GTGGAAAACGAAAAACCCg                 <  1:7131570/36‑1 (MQ=255)
                   tcatcaCTGGCGGT.GGC.GTGGAAAACGAAAAACCCg                 <  1:6539361/36‑1 (MQ=255)
                   tcatcaCTGGCGGT.GGC.GGGGGTAACGAAAAACCCg                 >  1:1225679/1‑36 (MQ=25)
                    catcaCTGGCGGT.GGG.GGGGGAAACGAAAAAaccag                >  1:6029049/1‑34 (MQ=16)
                    catcaCTGGCGGT.GGC.GTGGAAAACGAAAAACCCgg                <  1:7473732/36‑1 (MQ=255)
                     atcaCTGGCGGT.GGG.GGGGGAAACGAAAAAACCggg               >  1:3299536/1‑35 (MQ=16)
                     atcaCTGGCGGT.GGC.GTGGAAAACGAAAAACCCGGa               <  1:2583282/36‑1 (MQ=255)
                     atcaCTGGCGGT.GGC.GTGGAAAACGAAAAACCCGGa               <  1:3094306/36‑1 (MQ=255)
                     atcaCTGGCGGT.GGC.GTGGAAAACGAAAAACCCGGa               <  1:4618046/36‑1 (MQ=255)
                     atcaCTGGCGGT.GGC.GTGGAAAACGAAAAACCCGGa               <  1:77701/36‑1 (MQ=255)
                     atcaCTGGCGGT.GGC.GGGGTAAAAGAAAAAGCCGGa               >  1:1237295/1‑36 (MQ=19)
                      tcaCTGGCGGT.GGC.GTGGAAAACGAAAAACCCGGaa              <  1:2828155/36‑1 (MQ=255)
                       caCTGGCGGT.GGC.GTGGAAAACGAAAAACCCGGAAc             <  1:185852/36‑1 (MQ=255)
                       caCTGGCGGT.GGC.GTGGAAAACGAAAAACCCGGAAc             <  1:5128380/36‑1 (MQ=255)
                       caCTGGCGGT.GGC.GTGGAAAACGAAAAACCCGGAAc             <  1:541153/36‑1 (MQ=255)
                         cTGGCGGT.GGC.GTGGAAAACGAAAAACCCGGAACTg           >  1:7134056/1‑36 (MQ=255)
                         cTGGCGGT.GGC.GTGGAAAACGAAAAACCCGGAACTg           >  1:1212540/1‑36 (MQ=255)
                           ggCGGT.GGC.GTGGAAAACGAAAAACCCGGAACTgcg         >  1:1359553/1‑36 (MQ=255)
                           ggCGGG.GTC.GTGGACAACGAAAAACCCGGAACTgcg         <  1:1391968/36‑1 (MQ=25)
                            gCGGT.GGC.GTGGAAAACGAAAAACCCGGAACTgcgc        >  1:613485/1‑36 (MQ=255)
                            gCGGT.GGC.GTGGAAAACGAAAAACCCGGAACTgcgc        >  1:5388334/1‑36 (MQ=255)
                              ggT.GGG.GTGGAAAACGAAAAACCCGGAACTGCGCCa      <  1:8218082/36‑1 (MQ=255)
                              ggT.GGC.GTGGAAAACGAAAAACCCGGAACTGCGCCa      <  1:6790896/36‑1 (MQ=255)
                              ggT.GGC.GTGGAAAACGAAAAACCCGGAACTGCGCCa      >  1:1263024/1‑36 (MQ=255)
                              ggT.GGC.GTGGAAAACGAAAAACCCGGAACTGCGCCa      >  1:2732406/1‑36 (MQ=255)
                              ggT.GGC.GTGGAAAACGAAAAACCCGGAACTGCGCCa      <  1:1652415/36‑1 (MQ=255)
                              ggT.GGC.GTGGAAAACGAAAAACCCGGAACTGCGCCa      <  1:4265219/36‑1 (MQ=255)
                              ggT.GGC.GTGGAAAACGAAAAACCCGGAACTGCGCCa      >  1:3877753/1‑36 (MQ=255)
                               gT.GGC.GTGGAAAACGAAAAACCCGGAACTGCGCCAg     >  1:5796464/1‑36 (MQ=255)
                               gT.GGC.GTGGAAAACGAAAAACCCGGAACTGCGCCAg     <  1:7726371/36‑1 (MQ=255)
                                  ggC.GTGGAAAACGAAAAACCCGGAACTGCGCCAGAc   >  1:3495864/1‑36 (MQ=255)
                                   gC.GTGGAAAACGAAAAACCCGGAACTGCGCCAGACg  <  1:5452132/36‑1 (MQ=255)
                                     |                                  
TAACCTGTCGCTGGGCGGTTCATCACTGGCGGT.GGC.GTGGAAAACGAAAAACCCGGAACTGCGCCAGACG  >  REL606/2735394‑2735463

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 19 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG