Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,733,6490A→C 4.6% 0.2 ‑1.9 ‑0.097685pykFpyruvate kinase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (38/53):  total (39/58)
Fisher's exact test strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.95e-01

Only 100 of 122 total aligned reads displayed.

CATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATC  >  REL606/1733615‑1733684
                                  |                                   
atttcatcgtccaacaaaatgggtgcagtacaacaa                                    >  1:4471506‑M2/33‑36 (MQ=255)
atttcatcgtccaacaaaatgggtgcagtacaacaa                                    >  1:3997097‑M2/33‑36 (MQ=255)
 tttcatcgtccaacaaaatgggtgcagtacaacaac                                   >  1:936180‑M2/32‑36 (MQ=255)
 atcTCCAAAATCGAAAACCAGGAAGGCCTCAACAtg                                   <  1:3749872‑M2/36‑3 (MQ=255)
  ttcatcgtccaacaaaatgcgtgcagtacaACAACt                                  >  1:4190529‑M2/31‑36 (MQ=255)
   tcatcgtccaacaaaatgggtgcagtacaACAACtt                                 >  1:5134107‑M2/30‑36 (MQ=255)
   tcatcgtccaacaaaatgggtgcagtacaACAACtt                                 >  1:4625495‑M2/30‑36 (MQ=255)
   tcatcgtccaacaaaatgggtgcagtacaACAACtt                                 >  1:2531073‑M2/30‑36 (MQ=255)
   ctcCAAAATCGAAAACCAGGAAGGCCTCAACGtgta                                 >  1:1229949‑M2/1‑32 (MQ=255)
   caccaaAATCGAAAACCAGGAAGGCCTCAACAtgta                                 >  1:3820527‑M2/3‑32 (MQ=255)
    catcgtccaacaaaatgggtgcagtacaACAACTTc                                <  1:734171‑M2/8‑1 (MQ=255)
    catcgtccaacaaaatgggtgcagtacaACAACTTc                                <  1:5376590‑M2/8‑1 (MQ=255)
    catcgtccaacaaaatgggtgcagtacaACAACTTc                                <  1:1937326‑M2/8‑1 (MQ=255)
     ccAAAATCGAAAACCAGGAAGGCCTCAACAtgtact                               >  1:2339346‑M2/1‑30 (MQ=255)
     atcgtccaacaaaatgggtgcagtacaACAACTTcg                               >  1:289791‑M2/28‑36 (MQ=255)
      tcgtccaacaaaatgggtgcagtacaACAACTTcga                              <  1:25505‑M2/10‑1 (MQ=255)
        gtccaacaaaatgggtgcagtacaACAACTTcgacg                            >  1:5166559‑M2/25‑36 (MQ=255)
         aaTCGAAAACCAGGAAGGCCTCAACAtgtactgacc                           <  1:2686415‑M2/36‑11 (MQ=255)
          ccaacaaaatgggtgcagtacaACAACTTCGACGaa                          <  1:2992076‑M2/14‑1 (MQ=255)
          ccaacaaaatgggtgcagtacaACAACTTCGACGaa                          <  1:3649576‑M2/14‑1 (MQ=255)
          aTCGAAAACCAGGAAGGCCTCAACAtgtactgaccc                          >  1:1886920‑M2/1‑25 (MQ=255)
          aTCGAAAACCAGGAAGGCCTCAACAtgtactgaccc                          <  1:1447377‑M2/36‑12 (MQ=255)
            cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                        >  1:4555365‑M2/1‑23 (MQ=255)
            aacaaaatgggtgcagtacaACAACTTCGACGAAAt                        >  1:5531616‑M2/21‑36 (MQ=255)
             gAAAACCAGGAAGGCCTCAACAtgtactgaccccaa                       >  1:941291‑M2/1‑22 (MQ=255)
              caaaatgggtgcagtacaACAACTTCGACGAAATcc                      >  1:3213291‑M2/19‑36 (MQ=255)
              aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa                      <  1:4651056‑M2/36‑16 (MQ=255)
              aaaaCCAGGAAGGCCTCAACAtgtactgaccccaaa                      <  1:5221215‑M2/36‑16 (MQ=255)
               aaaatgggtgcagtacaACAACTTTGACGAAATCCt                     >  1:409887‑M2/18‑36 (MQ=255)
               aaaatgggtgcagtacaACAACTTCGACGAAATCCt                     <  1:4492684‑M2/19‑1 (MQ=255)
               aaaatgggtgcagtacaACAACTTCGACGAAATCCt                     <  1:2107502‑M2/19‑1 (MQ=255)
               aaaCCAGGAAGGCCTCAACAtgtactgaccccaaaa                     <  1:675657‑M2/36‑17 (MQ=255)
               aaaCCAGGAAGGCATCAACAtgtactgaccccaaaa                     <  1:3732215‑M2/36‑17 (MQ=255)
               aaaCCAGGAAGGCATCAACAtgtactgaccccaaaa                     >  1:3441423‑M2/1‑20 (MQ=255)
                aaatgggtgcagtacaACAACTTCGACGAAATCCTc                    >  1:3507338‑M2/17‑36 (MQ=255)
                aaatgggtgcagtacaACAACTTCGACGAAATCCTc                    >  1:1580148‑M2/17‑36 (MQ=255)
                aaatgggtgcagtacaACAACTTCGACGAAATCCTc                    >  1:1804980‑M2/17‑36 (MQ=255)
                aaCCAGGAAGGCCTCAACAtgtactgaccccaaaaa                    <  1:678913‑M2/36‑18 (MQ=255)
                 aCCAGGAAGGCCTCAACAtgtactgaccccaaaaag                   <  1:5000256‑M2/36‑19 (MQ=255)
                  atgggtgcagtacaACAACTTCGACGAAATCCTCGa                  >  1:5920655‑M2/15‑36 (MQ=255)
                   tgggtgcagtacaACAACTTCGACGAAATCCTCGaa                 <  1:3201918‑M2/23‑1 (MQ=255)
                    aGGAAGGCCTCAACAtgtactgaccccaaaaagttg                <  1:5912946‑M2/36‑22 (MQ=255)
                    aGGAAGGCCTCAACAtgtactgaccccaaaaagttg                <  1:291479‑M2/36‑22 (MQ=255)
                    aGGAAGGCCTCAACAtgtactgaccccaaaaagttg                <  1:2270580‑M2/36‑22 (MQ=255)
                    aGGAAGGCCTCAACAcgtactgaccccaaaaagttg                <  1:4367657‑M2/36‑22 (MQ=255)
                     ggAAGGCCTCAACAtgtactgaccccaaaaagttgg               <  1:3436889‑M2/36‑23 (MQ=255)
                     ggAAGGCCTCAACAtgtactgaccccaaaaagttgg               <  1:4480476‑M2/36‑23 (MQ=255)
                      gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc              >  1:4214161‑M2/11‑36 (MQ=255)
                      gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc              >  1:1069879‑M2/11‑36 (MQ=255)
                      gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc              <  1:3417889‑M2/26‑1 (MQ=255)
                      gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc              <  1:1892253‑M2/26‑1 (MQ=255)
                      gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc              <  1:2638503‑M2/26‑1 (MQ=255)
                      gAAGGCCTCAACAtgtactgaccccaaaaagttgga              <  1:224995‑M2/36‑24 (MQ=255)
                       aaGGCCTCAACAtgtactgaccccaaaaagttggac             >  1:1755164‑M2/1‑12 (MQ=255)
                       aaGGCCTCAACAtgtactgaccccaaaaagttggac             <  1:5642823‑M2/36‑25 (MQ=255)
                        gcagtacaACAACTTCGACGAAATCCTCGAAGCctc            <  1:4845299‑M2/28‑1 (MQ=255)
                        aGGCCTCAACAtgtactgaccccaaaaagttggaca            <  1:905139‑M2/36‑26 (MQ=255)
                         ggCCTCAACAtgtactgaccccaaaaagttggacag           >  1:1127081‑M2/1‑10 (MQ=255)
                         cagtacaACAACTTCGACGAAATCCTCGAAGCctct           <  1:509864‑M2/29‑1 (MQ=255)
                          gCCTCAACAtgtactgaccccaaaaagttggacagt          <  1:1864177‑M2/36‑28 (MQ=255)
                          gCCTCAACAtgtactgaccccaaaaagttggacagt          <  1:2055804‑M2/36‑28 (MQ=255)
                          agtacaACAACTTCGACGAAATCCTCGAAGCCTCTg          >  1:2195698‑M2/7‑36 (MQ=255)
                           gtacaACAACTTCGACGAAATCCTCGAAGCCTCTGa         >  1:4038898‑M2/6‑36 (MQ=255)
                           ccTCAACAtgtactgaccccaaaaagttggacagtt         >  1:3983917‑M2/1‑8 (MQ=255)
                           ccTCAACAtgtactgaccccaaaaagttggacagtt         >  1:2972698‑M2/1‑8 (MQ=255)
                            tacaACAACTTCGACGAAATCCTCGAAGCCTCTGAc        <  1:57344‑M2/32‑1 (MQ=255)
                            cTCAACAtgtactgaccccaaaaagttggacagtta        <  1:884641‑M2/36‑30 (MQ=255)
                            cTCAACAtgtactgaccccaaaaagttggacagtta        >  1:414481‑M2/1‑7 (MQ=255)
                             tCAACAtgtactgaccccaaaaagttggacagttaa       <  1:5311735‑M2/36‑31 (MQ=255)
                             tCAACAtgtactgaccccaaaaagttggacagttaa       <  1:4891742‑M2/36‑31 (MQ=255)
                             gcaACCtgtactgaccccaaaaagttggacagttaa       <  1:1131‑M2/35‑31 (MQ=255)
                              cAACAtgtactgaccccaaaaagttggacagttaaa      <  1:3132601‑M2/36‑32 (MQ=255)
                              cAACAtgtactgaccccaaaaagtcggacagttaaa      <  1:5144045‑M2/36‑32 (MQ=255)
                              cAACAtgtactgacaccaaaaagttggacagttaaa      <  1:1540506‑M2/36‑32 (MQ=255)
                               cgtatgtactgaccccaaaaagttggacagttaaac     <  1:2549298‑M2/33‑33 (MQ=255)
                               cgtatgtactgaccccaaaaagttggacagttaaac     >  1:2563226‑M2/4‑4 (MQ=255)
                               aaCCtgtactgaccccaaaaagttggacagttaaac     <  1:775052‑M2/36‑33 (MQ=255)
                               aaCCtgtactgaccccaaaaagttggacagttaaac     <  1:1088683‑M2/36‑33 (MQ=255)
                               aACAACTTCGACGAAATCCTCGAAGCCTCTGACGGc     <  1:2083104‑M2/35‑1 (MQ=255)
                                gtatgtactgaccccaaaaagttggacagttaaaca    >  1:5638412‑M2/3‑3 (MQ=255)
                                acaacTTCGACGAAATCCTCGAAGCCTCTGTCggaa    >  1:4462478/1‑34 (MQ=38)
                                aCAtgtactgaccccaaaaagttggacagttaaaca    <  1:5700712‑M2/36‑34 (MQ=255)
                                aCAtgtactgaccccaaaaagttggacagttaaaca    <  1:5642141‑M2/36‑34 (MQ=255)
                                ACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCa    >  1:2263186‑M2/1‑36 (MQ=255)
                                 tatgtactgaccccaaaaagttggacagttaaacac   <  1:4250910‑M2/35‑35 (MQ=255)
                                 tatgtactgaccccaaaaagttggacagttaaacac   <  1:4248795‑M2/35‑35 (MQ=255)
                                 tatgtactgaccccaaaaagttggacagttaaacac   >  1:5070787‑M2/2‑2 (MQ=255)
                                 cAtgtactgaccccaaaaagttggacagttaaacac   >  1:1102539‑M2/1‑2 (MQ=255)
                                 cAtgtactgaccccaaaaagttggacagttaaacac   <  1:4864568‑M2/36‑35 (MQ=255)
                                 CAACTTCGACAAAATCCTCGAAGCCTCTGACGGCAt   <  1:280533‑M2/36‑1 (MQ=255)
                                  atgtactgaccccaaaaagttggacagttaaacacg  <  1:2565483‑M2/36‑36 (MQ=255)
                                  atgtactgaccccaaaaagttggacagttaaacacg  >  1:3066004‑M2/1‑1 (MQ=255)
                                  atgtactgaccccaaaaagttggacagttaaacacg  <  1:5589793‑M2/36‑36 (MQ=255)
                                  atgtactgaccccaaaaagttggacagttaaacacg  <  1:2978426‑M2/36‑36 (MQ=255)
                                  atgtactgaccccaaaaagttggacagttaaacacg  >  1:5308791‑M2/1‑1 (MQ=255)
                                  atgtactgaccccaaaaagttggacagttaaacacg  <  1:852968‑M2/36‑36 (MQ=255)
                                  atgtactgaccccaaaaagttggacagttaaacacg  <  1:2958898‑M2/36‑36 (MQ=255)
                                  aacTTCGACGAAATCCTCGAAGCCTCTGACGGcatc  <  1:5433074/36‑1 (MQ=255)
                                  aacTTCGACGAAATCCTCGAAGCCTCTGACGGcatc  >  1:4652929/1‑36 (MQ=255)
                                  aacTTCGACGAAATCCTCAAAGCCTCTGACGGcatc  <  1:3488054/36‑1 (MQ=255)
                                  |                                   
CATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATC  >  REL606/1733615‑1733684

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