Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,174,7460T→C 3.6% 0.4 ‑2.6 ‑0.784965thiCthiamine biosynthesis protein ThiC
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (33/48):  total (36/48)
Fisher's exact test strand distribution p-value = 7.49e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.96e-01

CATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTA  >  REL606/4174711‑4174781
                                   |                                   
cATGGATAGTGAAGTANTCCAATCCTTGCTCGGCCt                                     <  1:2576493/36‑1 (MQ=37)
 aTGGATAGTGAAGTAATCCACACCTTGCTCGGCCTg                                    >  1:5988341/1‑36 (MQ=255)
  tGGTTAGTGAAGTAATCCACACCTTGCTCGGCCTGt                                   >  1:4732082/1‑36 (MQ=255)
  tGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGt                                   >  1:1099074/1‑36 (MQ=255)
  tGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGt                                   <  1:1263631/36‑1 (MQ=255)
  tGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGt                                   >  1:663732/1‑36 (MQ=255)
   ggATAGTGAAGTAATCCACACCTTGCTCGGCCTGtt                                  <  1:1772956/36‑1 (MQ=255)
   ggATAGTGAAGTAATCCACACCTTGCTCGGCCTGtt                                  <  1:7145594/36‑1 (MQ=255)
    gATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTc                                 >  1:5508102/1‑36 (MQ=255)
    gATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTc                                 >  1:6235074/1‑36 (MQ=255)
    gATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTc                                 <  1:6242359/36‑1 (MQ=255)
    gATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTc                                 <  1:920194/36‑1 (MQ=255)
    gATAGTGAAGTAATCCACACCTCGCTCGGCCTGTTc                                 >  1:1414023/1‑36 (MQ=255)
     aTAGTGAAGTAATCCACACCTTGCTCGGCCTGTTcc                                <  1:2119551/36‑1 (MQ=255)
      tAGTGAAGTAATCCACGCCTTGCTCGGCCTGTTCca                               <  1:2175473/36‑1 (MQ=255)
      tAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCca                               <  1:3909882/36‑1 (MQ=255)
      tAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCca                               <  1:6190786/36‑1 (MQ=255)
       aGTGAAGTAATCCACACCTTGCTCGGCCCGTTCcag                              >  1:7459893/1‑36 (MQ=255)
        gTGAAGTAATCCACACCTTGCTCGGCCTGTTCcagc                             <  1:508034/36‑1 (MQ=255)
        gTGAAGTAATCCACACCTTGCTCGGCCTGTTCcagc                             <  1:3581788/36‑1 (MQ=255)
        gTGAAGTAATCCACACCTTGCTCGGCCTGTTCcagc                             <  1:1425594/36‑1 (MQ=255)
        gTGAAGTAATCCACACCTTGCTCGGCCTGTTCcagc                             <  1:7337900/36‑1 (MQ=255)
         tGAAGTAATCCACACCTTGCTCGGCCTGTTCcagca                            <  1:6748203/36‑1 (MQ=255)
         tGAAGTAATCCACACCTTGCTCGGCCTGTTCcagca                            <  1:673278/36‑1 (MQ=255)
         tGAAGTAATCCACACCTTGCTCGGCCTGTTCcagca                            <  1:5948281/36‑1 (MQ=255)
         tGAAGTAATCCACACCTTGCTCGGCCTGTTCcagca                            <  1:1508977/36‑1 (MQ=255)
          gAAGTAATCCACACCTTGCTCGGCCTGTTCcagcag                           <  1:4431112/36‑1 (MQ=255)
          gAAGTAATCCACACCTTGCTCGGCCTGTTCcagcag                           >  1:4015061/1‑36 (MQ=255)
          gAAGTAATCCACACCTTGCTCGGCCTGTTCcagcag                           >  1:1564444/1‑36 (MQ=255)
           aaGTAATCCACACCTTGCTCGGCCTGTTCCagcagc                          <  1:5194800/36‑1 (MQ=255)
           aaGTAATCCACACCTTGCTCGGCCTGTTCCagcagc                          >  1:2573022/1‑36 (MQ=255)
              tAATCCACACCTTGCTCGGCCTGTTCCAGCAGCgtg                       >  1:6179085/1‑36 (MQ=255)
               aaTCCACACCTTGCTCGGCCTGTTCCAGCAGCgtgt                      >  1:3958826/1‑36 (MQ=255)
                aTCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTc                     >  1:5486011/1‑36 (MQ=255)
                aTCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTc                     <  1:3291445/36‑1 (MQ=255)
                aTCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTc                     >  1:5749902/1‑36 (MQ=255)
                aTCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTc                     >  1:4390225/1‑36 (MQ=255)
                aTCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTc                     <  1:6855495/36‑1 (MQ=255)
                aTCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTc                     <  1:6011647/36‑1 (MQ=255)
                aTCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTc                     >  1:2005560/1‑36 (MQ=255)
                 tCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTcg                    <  1:397695/36‑1 (MQ=255)
                 tCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTcg                    >  1:4652675/1‑36 (MQ=255)
                 tCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTcg                    <  1:5238379/36‑1 (MQ=255)
                   cacaCCTTGCTCGGCCTGTTCCAGCAGCGTGTcgcg                  >  1:4506512/1‑36 (MQ=255)
                     cacCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGa                <  1:2944825/36‑1 (MQ=255)
                     cacCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGa                <  1:2741050/36‑1 (MQ=255)
                     cacCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGa                <  1:6897338/36‑1 (MQ=255)
                      acCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGaa               >  1:850571/1‑36 (MQ=255)
                      acCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGaa               <  1:5791369/36‑1 (MQ=255)
                      acCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGaa               >  1:6998231/1‑36 (MQ=255)
                      acCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGaa               >  1:2422725/1‑36 (MQ=255)
                      acCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGaa               <  1:3535314/36‑1 (MQ=255)
                      acCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGaa               <  1:439474/36‑1 (MQ=255)
                      acCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGaa               >  1:101885/1‑36 (MQ=255)
                       ccTTGCTCGTCCTGTTCCAGCAGCGTGTCGCCGAAc              <  1:5522795/36‑1 (MQ=37)
                       ccTTGCTCGGCCTGTTCCAGCAGCGTGTTGCGGAAc              <  1:3150920/36‑1 (MQ=255)
                       ccTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAAc              <  1:2911474/36‑1 (MQ=255)
                        cTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACg             >  1:3305087/1‑36 (MQ=255)
                        cTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACg             <  1:5883878/36‑1 (MQ=255)
                        cTTGCTCGGCCTGCTCCAGNAGCGTGTCGCGGAACg             >  1:2697094/1‑36 (MQ=38)
                         ttGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGc            >  1:3319916/1‑36 (MQ=255)
                         ttGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGc            >  1:7493289/1‑36 (MQ=255)
                          tGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCt           >  1:753391/1‑36 (MQ=255)
                          tGCTCGGCCCGTTCCAGCAGCGTGTCGCGGAACGCt           >  1:2892478/1‑36 (MQ=255)
                           gCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCtt          <  1:1843895/36‑1 (MQ=255)
                           gCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCtt          >  1:3551244/1‑36 (MQ=255)
                           gCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCtt          <  1:4941035/36‑1 (MQ=255)
                           gCTCGGCCCGTTCCAGCAGCGTGTCGCGGAACGCtt          >  1:641999/1‑36 (MQ=255)
                            cTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTc         <  1:5195168/36‑1 (MQ=255)
                            cTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTc         >  1:4090453/1‑36 (MQ=255)
                            cTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTc         >  1:3927384/1‑36 (MQ=255)
                             tCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTcc        >  1:7424250/1‑36 (MQ=255)
                             tCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTcc        <  1:4155304/36‑1 (MQ=255)
                             tCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTcc        <  1:4392056/36‑1 (MQ=255)
                             tCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTcc        <  1:4514349/36‑1 (MQ=255)
                              cGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTccc       <  1:3867860/36‑1 (MQ=255)
                              cGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTccc       >  1:1023696/1‑36 (MQ=255)
                              cGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTccc       <  1:577893/36‑1 (MQ=255)
                              cGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTccc       <  1:4556015/36‑1 (MQ=255)
                              cGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTccc       <  1:4635830/36‑1 (MQ=255)
                               ggCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCa      <  1:6180634/36‑1 (MQ=255)
                               ggCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCa      <  1:594021/36‑1 (MQ=255)
                               ggCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCa      >  1:5052423/1‑36 (MQ=255)
                                gCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAg     <  1:7359732/36‑1 (MQ=255)
                                  cTGTTCCAGCAGGGTGTCGCGGAACGCTTCCCAGGt   <  1:2451293/36‑1 (MQ=255)
                                   tGTTCCAGCAGCGTGTCGCGGAACGCTTCCCCGGTa  >  1:3597060/1‑36 (MQ=255)
                                   tGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTa  >  1:6177474/1‑36 (MQ=255)
                                   |                                   
CATGGATAGTGAAGTAATCCACACCTTGCTCGGCCTGTTCCAGCAGCGTGTCGCGGAACGCTTCCCAGGTA  >  REL606/4174711‑4174781

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG