Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,231,6190T→C 8.4% 0.7 ‑5.3 0.076182ubiCchorismate pyruvate lyase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (37/30):  total (37/39)
Fisher's exact test strand distribution p-value = 5.19e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGTCTC  >  REL606/4231584‑4231654
                                   |                                    
gatgatCCGCGAAGGGTTTGTCGAGCAGAATGAAAt                                      <  1:1296134/36‑1 (MQ=255)
 atgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATc                                     >  1:1820223/1‑36 (MQ=255)
 atgatCCGCGAAGGGTTTATCGAGCAGAATGAAATc                                     <  1:1923218/36‑1 (MQ=255)
  tgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATcc                                    <  1:5184081/36‑1 (MQ=255)
  tgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATcc                                    <  1:1617024/36‑1 (MQ=255)
  tgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATcc                                    >  1:7075270/1‑36 (MQ=255)
  tgatCCGCGAAGGGTTTGTCGAGCAGAATGAAATcc                                    >  1:6011475/1‑36 (MQ=255)
  tgatCCGCGAAGGGCTTGTCGAGCAGAATGAAATcc                                    >  1:3070492/1‑36 (MQ=255)
    atCCGCGAAGGGTTTGTCGAGCAGAATGAAATcccc                                  >  1:438127/1‑36 (MQ=255)
     tCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCC.g                                <  1:3146230/36‑1 (MQ=255)
     tCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCC.g                                <  1:62223/36‑1 (MQ=255)
      ccGCGAAGGGTTTGTCGAGCAGAATGAAATCCCC.ga                               <  1:3300668/36‑1 (MQ=255)
       cgcgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaa                              >  1:2207599/1‑36 (MQ=255)
        gcgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaag                             <  1:5701730/36‑1 (MQ=255)
        gcgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaag                             <  1:3003671/36‑1 (MQ=255)
        gcgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaag                             >  1:7588897/1‑36 (MQ=255)
         cgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaaga                            >  1:4050961/1‑36 (MQ=255)
         cgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaaga                            <  1:5372313/36‑1 (MQ=255)
         cgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaaga                            <  1:3791239/36‑1 (MQ=255)
         cgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaaga                            >  1:7086415/1‑36 (MQ=255)
         cgAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaaga                            <  1:6614411/36‑1 (MQ=255)
         cgAAGGGTTTGTAGAGCAGAATGAAATCCCC.gaaga                            <  1:1729607/36‑1 (MQ=255)
          gAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaagaa                           <  1:4448258/36‑1 (MQ=255)
          gAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaagaa                           >  1:5500895/1‑36 (MQ=255)
          gAAGGGTTTGTCGAGCAGAATGAAATCCCC.gaagaa                           <  1:3687686/36‑1 (MQ=255)
          gAAGGGTATGTCGAGCAGAATGAAATCCCC.gaagaa                           <  1:5766309/36‑1 (MQ=255)
           aaGGGTTTGTCGAGCAGAATGAAATCCCC.GAAGAAc                          <  1:5821870/36‑1 (MQ=255)
            aGGGTTTGTCGAGCAGAATGAAATCCCC.GAAGAACt                         >  1:2235948/1‑36 (MQ=255)
              ggTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGc                       >  1:2553364/1‑36 (MQ=255)
              ggTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGc                       <  1:2958566/36‑1 (MQ=255)
               gTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGcc                      <  1:387084/36‑1 (MQ=255)
               gTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGcc                      >  1:1756389/1‑36 (MQ=255)
               gTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGcc                      <  1:4765371/36‑1 (MQ=255)
                tttGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCg                     >  1:487098/1‑36 (MQ=255)
                tttGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCg                     >  1:1241407/1‑36 (MQ=255)
                tttGTCGAGAAGAATGAAATCCCC.GAAGAACTGCCg                     <  1:1956424/36‑1 (MQ=255)
                 ttGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCgc                    >  1:1627714/1‑36 (MQ=255)
                 ttGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCgc                    >  1:6345918/1‑36 (MQ=255)
                 ttGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCgc                    <  1:3042340/36‑1 (MQ=255)
                  tGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCgct                   >  1:3308157/1‑36 (MQ=255)
                   gTCGAGCAGAATGAAATCCCC.GAAGAACTGCCgctg                  >  1:4774873/1‑36 (MQ=255)
                   gTCGAGCAGAATGAAACCCCC.GAAGAACTGCCgctg                  <  1:7000058/36‑1 (MQ=255)
                    tcaaGCACAAAGAAACCCCC.GAAGATCTGCCgctgc                 <  1:5432075/33‑1 (MQ=12)
                     cGAGCAGAATGAAATCCCC.GAAGAACTGCCgctgct                >  1:3643663/1‑36 (MQ=255)
                     cGAGAAGAATGAAATCCCC.GAAGAACTGCCgctgct                <  1:1354217/36‑1 (MQ=255)
                      taGCAGAATGAAATCCCC.TAAGAACTGCCGctgctg               <  1:1028800/35‑1 (MQ=38)
                      gAGCAGAATGAAATCCCC.GAAGAACTGCCGctgctg               >  1:928545/1‑36 (MQ=255)
                      gAGAAGAATGAAATCCCC.GAAGAACTGCCGctgctg               <  1:2289630/36‑1 (MQ=255)
                       aGCAGAATGAAATCCCC.GAAGAACTGCCGCtgctgc              >  1:6947472/1‑36 (MQ=255)
                       aGCAGAATGAAATCCCC.GAAGAACTGCCGCtgctgc              <  1:39330/36‑1 (MQ=255)
                       aGCAGAATGAAATCCCC.GAAGAACTGCCGCtgctgc              >  1:7107741/1‑36 (MQ=255)
                       aGCAGAAGGAAACCCCC.GAGGATCTGCCGCtgctgc              <  1:4732461/36‑1 (MQ=16)
                        gCAGAATGAAATCCCC.GAAGAACTGCCGCTGCTGcc             <  1:3850201/36‑1 (MQ=255)
                        gCAGAATGAAATCCCC.GAAGAACTGCCGCTGCTGcc             >  1:5517815/1‑36 (MQ=255)
                        gCAGAATGAAATCCCC.GAAGAACTGCCGCTGCTGcc             >  1:7817380/1‑36 (MQ=255)
                        gCAGAATGAAATCCCC.GAAGAACTGCCGCTGCTGcc             >  1:5175558/1‑36 (MQ=255)
                        gCAAAATAAAACCCCC.GAAGAACTGCCGCTGCTGcc             <  1:928737/36‑1 (MQ=25)
                          aaaaTGAAACCCCC.GAAAAGCTGCCGCTGCTGCCGa           <  1:4969901/34‑1 (MQ=21)
                          aGAATGAAATCCCC.GAAGAACTGCCGCTGCTGCCGa           >  1:3490657/1‑36 (MQ=255)
                          aGAATGAAATCCCC.GAAGAACTGCCGCTGCTGCCGa           >  1:2057780/1‑36 (MQ=255)
                          aGAATGAAATCCCC.GAAGAACTGCCGCTGCTGCCGa           >  1:1826831/1‑36 (MQ=255)
                           gAATGAAATCCCC.GAAGAACTGCCGCTGCTGCCGaa          >  1:6252613/1‑36 (MQ=255)
                           aaaTGAAACCCCC.GAAGAACTGCCGCTGCTGCCGaa          <  1:7624821/35‑1 (MQ=255)
                            aaTGAAATCCCC.GAAGAACTGCCGCTGCTGCCGaaa         >  1:3660818/1‑36 (MQ=255)
                            aaTAAAATCCCC.AAAGAACTGCCGCTGCTGCCGaaa         <  1:7293996/36‑1 (MQ=37)
                             cataaaTCCCCAGAATAACTGCCGCTGCTGCCGaaa         <  1:426905/33‑1 (MQ=19)
                             aTGAAATCCCC.GAAGAACTGCCGCTGCTGCCGAAag        >  1:148389/1‑36 (MQ=255)
                              taaaaCCCCC.AAAGAACTGCCGCTGCTGCCGAAaga       <  1:6289103/34‑1 (MQ=25)
                              tGAAATCCCC.GAAGAACTGCCGCTGCTGCTGAAaga       >  1:7914348/1‑36 (MQ=255)
                              cgAAAACCCC.GAAAAACTGCCGCTGCTGCCGAAaga       <  1:1651026/35‑1 (MQ=37)
                               aaaaTCCCC.GAAGAACTGCCGCTGCTGCCGAAagag      <  1:714315/35‑1 (MQ=255)
                                aaaCCCCC.AAAGAACTGCCGCTGCTGCCGAAAGAGt     <  1:6624950/36‑1 (MQ=37)
                                 aaTCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGtc    >  1:2215552/1‑36 (MQ=255)
                                 aaTCCCC.AAAAAACTGCCGCTGCTGCCGAAAGAGtc    <  1:992855/36‑1 (MQ=37)
                                  aTCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGtct   >  1:5729973/1‑36 (MQ=255)
                                  aTCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGtct   >  1:3992030/1‑36 (MQ=255)
                                  aTCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGtct   >  1:2535551/1‑36 (MQ=255)
                                   tCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGtctc  >  1:1871226/1‑36 (MQ=255)
                                   |                                    
GATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCC.GAAGAACTGCCGCTGCTGCCGAAAGAGTCTC  >  REL606/4231584‑4231654

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG