Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,807,0330T→. 2.9% 0.5 ‑1.6 ‑0.070coding (67/1392 nt)uhpTsugar phosphate antiporter
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (29/38):  total (31/39)
Fisher's exact test strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01

CACCAGGTAGGATTGCATGAACGGTTTGAACCACATTTTGCGCCGCACTTCGAGCGGAAGGTCCAGGGTCG  >  REL606/3806998‑3807068
                                   |                                   
caccaGGTAGGATTGCATTAACGGTTTGAACCACAt                                     <  1:5705773/36‑1 (MQ=255)
caccaGGTAGGATTGCATGAACGGTTTGAACCACAt                                     <  1:3168478/36‑1 (MQ=255)
caccaGGTAGGATTGCATGAACGGTTTGAACCACAt                                     <  1:7100337/36‑1 (MQ=255)
 accgGGTATGATTGCATGAACGGTTTGAACCACAtt                                    <  1:2690895/36‑1 (MQ=37)
 accaGGTAGGATTGCATGAACGGTTTGAACCACAtt                                    >  1:436172/1‑36 (MQ=255)
 accaGGTAGGATTGCATGAACGGTTTGAACCACAtt                                    >  1:7797827/1‑36 (MQ=255)
 accaGGTAGGATTGCATGAAAGGTTTGAACCACAtt                                    >  1:2232509/1‑36 (MQ=255)
 accaGGTAGGATTGAATGAACGGTTTGAACCACAtt                                    <  1:4703452/36‑1 (MQ=255)
  ccgggTAGGATTGCCTGAACGGTTTGAACCACAttt                                   <  1:2615426/33‑1 (MQ=37)
  ccaGGTAGGATTGCATGAACGGTTTGAACCACAttt                                   <  1:2539282/36‑1 (MQ=255)
  ccaGGTAGGATTGCATGAACGGTTTGAACCACAttt                                   <  1:1984431/36‑1 (MQ=255)
  ccaGGTAGGATTGCATGAACGGTTTGAACCACAttt                                   <  1:1122400/36‑1 (MQ=255)
   caGGTAGGATTGCATGAACGGTTTGAACCACAtttt                                  <  1:508284/36‑1 (MQ=255)
    aGGTAGGATTGCATGAACGGTTTGAACCACATTTTg                                 <  1:4095845/36‑1 (MQ=255)
    aGGTAGGATTGCATGAACGGTTTGAACCACATTTTg                                 <  1:7233105/36‑1 (MQ=255)
     ggTAGGATTGCATGAACGGTTTGAACCACATTTTgc                                >  1:3861896/1‑36 (MQ=255)
     ggTAGGATTGCATGAACGGTTTGAACCACATTTTgc                                <  1:2167004/36‑1 (MQ=255)
     ggTAGGATTGCATGAACGGTTTGAACCACATTTTgc                                <  1:7749647/36‑1 (MQ=255)
     ggTAGGATTGCATGAACGGTTTGAACCACATTTTgc                                <  1:1695569/36‑1 (MQ=255)
     ggTAGGATTGCATGAACGGTTTGAACCACATTTTga                                >  1:5220931/1‑35 (MQ=255)
      gTAGGATTGCATGAACGGTTTGAACCACATTTTgcg                               >  1:2487391/1‑36 (MQ=255)
      gTAGGATTGCATGAACGGTTTGAACCACATTTTgcg                               >  1:6142669/1‑36 (MQ=255)
      gTAGGATTGCATGAACGGTTTGAAACACATTTTgcg                               <  1:575631/36‑1 (MQ=255)
       tAGGATTGCATGAACGGTTTGAACCACATTTTgcgc                              <  1:7118631/36‑1 (MQ=255)
        aGGATTGCATGAACGGTTTGAACCGCATTTTGcgcc                             <  1:2333023/36‑1 (MQ=255)
        aGGATTGCATGAACGGTTTGAACCACATTTTGcgcc                             >  1:7422370/1‑36 (MQ=255)
        aGGATTGCATGAACGGTTTGAACCACATTTTGcgcc                             <  1:3404742/36‑1 (MQ=255)
        aGGATTGCATGAACGGTTTGAACCACATTTTGcgcc                             <  1:4441613/36‑1 (MQ=255)
         ggATTGCATGAACGGTTTGAACCACATTTTGcgccg                            <  1:7001180/36‑1 (MQ=255)
          gATTGCATGAACGGTTTGCACCACATTTTGcgccgc                           >  1:3895812/1‑36 (MQ=255)
          gATTGCATGAACGGTTTGAACCACATTTTGcgccgc                           >  1:5397827/1‑36 (MQ=255)
          gATTGCATGAACGGTTTGAACCACATTTTGcgccgc                           <  1:2416600/36‑1 (MQ=255)
          gATTGCATGAACGGTTTGAACCACATTTTGcgccgc                           <  1:3414307/36‑1 (MQ=255)
          gATTGCATGAACGGTTTGAACCACATTTTGcgccgc                           <  1:3741823/36‑1 (MQ=255)
            ttGCATGAACGGTTTGAACCACATTTTGCGCCGCAc                         <  1:4384401/36‑1 (MQ=255)
            ttGCATGAACGGTTTGAACCACATTTTGCGCCGCAc                         <  1:5902032/36‑1 (MQ=255)
            ttGCATGAACGGTTTGAACCACATTTTGCGCCGCAc                         >  1:2435594/1‑36 (MQ=255)
             tGCATGAACGGTTTGAACCACATTTTGCGCCGCACt                        >  1:4336026/1‑36 (MQ=255)
             tGCATGAACGGTTTGAACCACATTTTGCGCCGCACt                        >  1:5175719/1‑36 (MQ=255)
             tGCATGAACGGTTTGAACCACACTTTGCGCCGCACt                        <  1:3174386/36‑1 (MQ=255)
              gCATGAACGGTTTGAACCACATTTTGCGCCGCACtt                       >  1:5276136/1‑36 (MQ=255)
              gCATGAACGGTTTGAACCACATTTTGCGCCGCACtt                       <  1:1201814/36‑1 (MQ=255)
              gCATGAACGGTTTGAACCACATTTTGCGCCGCACtt                       <  1:2538013/36‑1 (MQ=255)
               cATGAACGGTTTGAACCACATTTTGCGCCGCACTTc                      <  1:1915930/36‑1 (MQ=255)
               cATGAACGGTTTGAACCACATTTTGCGCCGCACTTc                      <  1:5226090/36‑1 (MQ=255)
                aTGAACGGTTTGAACCACATTTTGCGCCGCACTTCg                     >  1:741243/1‑36 (MQ=255)
                aTGAACGGTTTGAACCACATTTTGCGCCGCACTTCg                     >  1:5032543/1‑36 (MQ=255)
                 tGAACGGTTTGAACCACATTTTGCGCCGCACTTCGa                    <  1:7705502/36‑1 (MQ=255)
                  gAACGGTTTGAACCACATTTTGCGCCGCACTTCGAg                   <  1:5243250/36‑1 (MQ=255)
                  gAACGGTTTGAACCACATTTTGCGCCGCACTTCGAg                   >  1:3419456/1‑36 (MQ=255)
                   aaCGGTTTGAACCACATTTTGCGCCGCACTTCGAGc                  <  1:5225402/36‑1 (MQ=255)
                   aaCGGTTTGAACCACATTTTGCGCCGCACTTCGAGc                  >  1:4780312/1‑36 (MQ=255)
                   aaCGGTTTGAACCACATTTTGCGCCGCACTTCGAGc                  >  1:7337788/1‑36 (MQ=255)
                   aaCGGTTTGAACCACATTTTGCGCCGCACTTCGAGc                  <  1:6061914/36‑1 (MQ=255)
                   aaCGGTTTGAACCACATTTTGCGCCGCACTTCGAGc                  <  1:1696576/36‑1 (MQ=255)
                   aaCGGTTTGAACCACATTTTGCGCCGCACTTCGAGc                  >  1:3948641/1‑36 (MQ=255)
                   aaCAGTTTGAACCACATTTTGCGCCGCACTTCGAGc                  <  1:3019839/36‑1 (MQ=255)
                    aCGGTTTGAACCACATTTTGCGCCGCACTTCGAGCg                 >  1:1516241/1‑36 (MQ=255)
                    aCGGTTTGAACCACATTTTGCGCCGCACTTCGAGCg                 >  1:6369877/1‑36 (MQ=255)
                    aCGGTTTGAACCACATTTTGCGCCGCACTTCGAGCg                 >  1:388357/1‑36 (MQ=255)
                    aCGGTTTGAACCACATTTTACGGCGCACTTCGAGCg                 >  1:1755929/1‑36 (MQ=37)
                     cGGTTTGAACCACATTTTGCTCCGCACTTCGAGCgg                <  1:4673624/36‑1 (MQ=255)
                     cGGTTTGAACCACATTTTGCGCCGCACTTCGAGCgg                >  1:4681682/1‑36 (MQ=255)
                      ggTTTGAACCACATTTTGCGCCGCACTTCGAGCGGa               >  1:4775998/1‑36 (MQ=255)
                      ggTTTGAACCACATTTTGCGCCGCACTTCGAGCGGa               <  1:1403572/36‑1 (MQ=255)
                        tttGAACCACATTTTGCGCCGCACTTCGAGCGGAAg             <  1:7633460/36‑1 (MQ=255)
                         ttGAACCACATTTTGCGCCGCACTTCGAGCGGAAgg            >  1:3284320/1‑36 (MQ=255)
                         ttGAACCACATTTTGCGCCGCACTTCGAGCGGAAgg            <  1:7248115/36‑1 (MQ=255)
                          tGAACCACATTTTGCGCCGCACTTCGAGCGGAAGGt           <  1:2084884/36‑1 (MQ=255)
                            aaCCACATTTTGCGCCGCACTTCGAGCGGAAGGTcc         <  1:7744361/36‑1 (MQ=255)
                            aaCCACATTTTGCGCCGCACTTCGAGCGGAAGGTcc         <  1:911309/36‑1 (MQ=255)
                            aaCCACA.TTTGCGCCGCACTTCGAGCGGTAGGTCCa        >  1:7396355/1‑36 (MQ=37)
                            aaCCACA.TTTGCGCCGCACTTCGAGCGGAAGGTCCa        >  1:518407/1‑36 (MQ=38)
                             aCCACATTTTGCGCCGCTCTTCGAGCGGAAGGTCCa        <  1:7765835/36‑1 (MQ=255)
                             aCCACATTTTGCGCCGCACTTCGAGCGGAAGGTCCa        <  1:4329626/36‑1 (MQ=255)
                              ccACATTTTGCGCCGCACTTCGAGCGGAAGGTCCAg       <  1:2112467/36‑1 (MQ=255)
                               cacaTTTTGCGCCGCACTTCGAGCGGAAGGTCCAgg      >  1:4175387/1‑36 (MQ=255)
                                acaTTTTGCGCCGCACTTCGAGCGGAAGGTCCAggg     >  1:1899038/1‑36 (MQ=255)
                                acaTTTTGCGCCGCACTTCGAGCGGAAGGTCCAggg     >  1:6454252/1‑36 (MQ=255)
                                acaTTTTGCGCCGCACTTCGAGCGGAAGGTCCAggg     >  1:5121033/1‑36 (MQ=255)
                                 caTTTTGCGCCGCACTTCGAGCGGAAGGTCCGGGGt    <  1:5201996/36‑1 (MQ=255)
                                  aTTTTGCGCCGCACTTCGAGCGGAAGGTCCAGGGTc   >  1:3474457/1‑36 (MQ=255)
                                  aTTTTGCGCCGCACTTCGAGCGGAAGGTCCAGGGTc   <  1:3735183/36‑1 (MQ=255)
                                   ttttGCGCCGCACTTCGAGGGGAAGGTCGAGGGTca  >  1:6257149/1‑35 (MQ=37)
                                   ttttGCGCCGCACTTCGAGCGGAAGGTCCAGGGTCg  >  1:2747810/1‑36 (MQ=255)
                                   |                                   
CACCAGGTAGGATTGCATGAACGGTTTGAACCACATTTTGCGCCGCACTTCGAGCGGAAGGTCCAGGGTCG  >  REL606/3806998‑3807068

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG