Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,929,9200T→. 2.5% 0.2 0.0 ‑0.081coding (319/567 nt)yecDpredicted hydrolase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (1/1):  ref base (36/43):  total (37/44)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.92e-01

GCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGC  >  REL606/1929885‑1929955
                                   |                                   
gCGATATCGAAATCATCAAACGTCAATGGGGTGCGt                                     >  1:7360966/1‑36 (MQ=255)
gCGATATCGAAATCATCAAACGTCAATGGGGTGCGt                                     <  1:519912/36‑1 (MQ=255)
gCGATATCGAAATCATCAAACGTCAATGGGGTGCGt                                     <  1:6291014/36‑1 (MQ=255)
  gATATCGAAATCATCAAACGTCACTGGGGTGCGttt                                   <  1:7554259/36‑1 (MQ=255)
  gATATCGAAATCATCAAACGTCAATGGGGTGCGttt                                   >  1:7878967/1‑36 (MQ=255)
   atatCGAAATCATCAAACGTCAATGGGGTGCGtttt                                  <  1:2944965/36‑1 (MQ=255)
   atatCGAAATCATCAAACGTCAATGGGGTGCGtttt                                  <  1:80335/36‑1 (MQ=255)
   atatCGAAATCATCAAACGTCAATGGGGTGCGtttt                                  >  1:3170871/1‑36 (MQ=255)
   atatCGAAATCATCAAACGTCAATGGGGTGCGtttt                                  <  1:1764385/36‑1 (MQ=255)
    tatCGAAATCATCAAACGTCAATGGGGTGCGTTTta                                 <  1:4597808/36‑1 (MQ=255)
    tatCGAAATCATCAAACGTCAATGGGGTGCGTTTta                                 >  1:2355011/1‑36 (MQ=255)
    tatCGAAATCATCAAACGTCAATGGGGTGCGTTTta                                 <  1:6022962/36‑1 (MQ=255)
    tatCGAAATCATCAAACGTCAATGGGGTGCGTTTta                                 >  1:6109213/1‑36 (MQ=255)
     atCGAAATCATCAAACGTCAATGGGGTGCGTTTtac                                <  1:6642186/36‑1 (MQ=255)
     atCGAAATCATCAAACGTCAATGGGGTGCGTTTtac                                >  1:2707586/1‑36 (MQ=255)
     atCGAAATCATCAAACGTCAATGGGGTGCGTTTtac                                <  1:7341701/36‑1 (MQ=255)
     atCGAAATCATCAAACGTCAATGGGGTGCGTTTtac                                <  1:908553/36‑1 (MQ=255)
      tCGAAATCATCAAACGTCAATGGGGTGCGTTTtacg                               <  1:6529291/36‑1 (MQ=255)
      tCGAAATCATCAAACGTCAATGGGGTGCGTTTtacg                               <  1:6993989/36‑1 (MQ=255)
      tCGAAATCATCAAACGTCAATGGGGTGCGTTTtacg                               <  1:4040970/36‑1 (MQ=255)
       cGAAATCATCAAACGTCAATGGGGTGCGTTTtacgg                              <  1:5037145/36‑1 (MQ=255)
        gAAATCATCAAACGTCACTGGGGTGCGTTTtacggt                             <  1:2858678/36‑1 (MQ=255)
        gAAATCATCAAACGTCAATGGGGTGCGTTTtacggt                             <  1:5317654/36‑1 (MQ=255)
        gAAATCATCAAACGTCAATGGGGTGCGTTTtacggt                             <  1:3579184/36‑1 (MQ=255)
        gAAATCATCAAACGTCAATGGGGTGCGTTTtacggt                             <  1:1043471/36‑1 (MQ=255)
        gAAATCATCAAACGTCAATGGGGTGCGTTTtacggt                             >  1:1456714/1‑36 (MQ=255)
        gAAATCATCAAACGTCAATGGGGTGCGTTTtacgga                             >  1:1215047/1‑35 (MQ=255)
         aaaTCATCAAACGTCAATGGGGTGCGTTTtacggta                            <  1:6186845/36‑1 (MQ=255)
         aaaTCATCAAACGTCAATGGGGTGCGTTTtacggta                            >  1:2055752/1‑36 (MQ=255)
         aaaTCATCAAACGACAATGGGGTGCGTTTtacggta                            >  1:815412/1‑36 (MQ=255)
          aaTCATCAAACGTCAATGGGGTGCGTTTtacggtac                           <  1:1302449/36‑1 (MQ=255)
          aaTCATCAAACGTCAATGGGGTGCGTTTtacggtac                           <  1:4012414/36‑1 (MQ=255)
          aaTCATCAAACGTCAATGGGGTGCGTTTtacggtac                           >  1:2783301/1‑36 (MQ=255)
          aaTCATCAAACGTCAATGGGGTGCGTTTtacggtac                           <  1:7820223/36‑1 (MQ=255)
          aaTCATCAAACGTCAATGGGGTGCGTTTtacggtac                           <  1:5750459/36‑1 (MQ=255)
           atcatcAAACGTCAATGGGGTGCTTTTtacggtacg                          >  1:7752684/1‑36 (MQ=255)
           atcatcAAACGTCAATGGGGTGCGTTTtacggtacg                          <  1:3157498/36‑1 (MQ=255)
           atcatcAAACGTCAATGGGGTGCGTTTtacggtacg                          >  1:2782373/1‑36 (MQ=255)
            tcatcaAACGTCAATGGGGTGAGTTTtacggtacgg                         >  1:1725859/1‑36 (MQ=255)
             catcaAACGTCAATGGGGTGCGTTTTACGGTACGGa                        <  1:3714397/36‑1 (MQ=255)
             catcaAACGTCAATGGGGTGCGTTTTACGGTACGGa                        >  1:3228240/1‑36 (MQ=255)
              atcaAACGTCAATGGGGTGCGTTTtaccgtacggtt                       >  1:7786721/1‑34 (MQ=37)
              atcaAACGTCAATGGGGTGCGTTTTACGGTACGGAt                       <  1:6624247/36‑1 (MQ=255)
              atcaAACGTCAATGGGGTGCGTTTTACGGTACGGAt                       <  1:3828495/36‑1 (MQ=255)
              atcaAACGTCAATGGGGTGCGTTTTACGGTACGGAt                       <  1:4822687/36‑1 (MQ=255)
              atcaAACGTCAATGGGGTGCGTTTTACGGTACGGAt                       <  1:1504794/36‑1 (MQ=255)
              atcaAACGTCAATGGGGTGCGTTTTACGGCACGGAt                       >  1:3185379/1‑36 (MQ=255)
               tcaAACGTCAATGGGGTGCGTTTTACGGTACGGATc                      <  1:4933992/36‑1 (MQ=255)
               tcaAACGTCAATGGGGTGCGTTTTACGGTACGGATc                      <  1:2874981/36‑1 (MQ=255)
               tcaAACGTCAATGGGGTGCGTTTTACGGTACGGATc                      <  1:5761259/36‑1 (MQ=255)
                caAACGTCAATGGGGTGCGTTTTACGGTACGGATCt                     <  1:6375734/36‑1 (MQ=255)
                 aaaCGTCAATGGGGTGCGTTTTACGGTACGGATCTg                    <  1:7159859/36‑1 (MQ=255)
                 aaaCGTCAATGGGGTGCGTTTTACGGTACGGATCTg                    <  1:4641855/36‑1 (MQ=255)
                  aaCGTCAATGGGGTGCGTTTTACGGTACGGATCTgg                   >  1:18095/1‑36 (MQ=255)
                  aaCGTCAATGGGGTGCGTTTTACGGTACGGATCTgg                   >  1:7250979/1‑36 (MQ=255)
                  aaCGTCAATGGGGTGCGTTTTACGGTACGGATCTgg                   >  1:3020862/1‑36 (MQ=255)
                  aaCGTCAATGGGGTGCGTTTTACGGTACGGATCTgg                   >  1:5067060/1‑36 (MQ=255)
                  aaCGTCAATGGGGTGCG.TTTACGGTACGGATCTGGa                  >  1:5490468/1‑36 (MQ=37)
                   aCGTCAATGGGGTGCGTTTTACGGTACGGATCTGGa                  >  1:1430078/1‑36 (MQ=255)
                   aCGTCAATGGGGTGCGTTTTACGGTACGGATCTGGa                  <  1:4687402/36‑1 (MQ=255)
                    cGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAg                 >  1:1342164/1‑36 (MQ=255)
                    cGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAg                 >  1:7005961/1‑36 (MQ=255)
                     gTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGt                <  1:4497199/36‑1 (MQ=255)
                     gTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGt                >  1:2101624/1‑36 (MQ=255)
                      tCAATGGGGTGCGTTTTACGGTACGGATCTGGAGtt               >  1:7080672/1‑36 (MQ=255)
                       cAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTg              <  1:6430289/36‑1 (MQ=255)
                       cAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTg              >  1:6614021/1‑36 (MQ=255)
                        aaTGGGGTGCGTTTTACGTTACGGATCTGGAGTTGc             >  1:6314548/1‑36 (MQ=255)
                        aaTGGGGTGCGTTTTACGGTACGGATCTGGAGTTGc             >  1:4920617/1‑36 (MQ=255)
                        aaTGGGGTGCGTTTTACGGTACGGATCTGGAGTTGc             >  1:5730454/1‑36 (MQ=255)
                        aaTGGGGTGCGTTTTACGGTACGGATCTGGAGTTGc             >  1:601076/1‑36 (MQ=255)
                        aaTGGGGTGCGTTTTACGGTACGGATCTGGAGTTGc             >  1:5886261/1‑36 (MQ=255)
                         aTGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCa            <  1:1649384/36‑1 (MQ=255)
                         aTGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCa            <  1:786110/36‑1 (MQ=255)
                         aTGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCa            >  1:449611/1‑36 (MQ=255)
                         aTGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCa            <  1:7615532/36‑1 (MQ=255)
                           ggggTGCGTTTTACGGTACGGATCTGGAGTTGCAAt          >  1:5023291/1‑36 (MQ=255)
                           ggggTGCGTTTTACGGTACGGATCTGGAGTTGCAAt          >  1:7430700/1‑36 (MQ=255)
                            gggTGCGTTTTACGGTACGGATCTGGAGTTGCAAtt         >  1:2557647/1‑36 (MQ=255)
                            gggTGCGTTTTACGGTACGGATCTGGAGTTGCAAtt         <  1:4320190/36‑1 (MQ=255)
                            gggTGCGTTTTACGGTACGGATCTGGAGTTGCAAtt         <  1:749938/36‑1 (MQ=255)
                             ggTGCG.TTTACGGTACGGATCTGGAGTTGCAAtttc       <  1:2484235/36‑3 (MQ=37)
                              gTGCGTTTTACGGTACGGATCTGGAGTTGCAATTAc       <  1:454685/36‑1 (MQ=255)
                                gCGTTTTACGGTACGGATCTGGAGTTGCAATTAcgc     >  1:3444015/1‑36 (MQ=255)
                                 cGTTTTACGGTACGGATCTGGAGTTGCAATTAcgcc    <  1:2808327/36‑1 (MQ=255)
                                 cGTTTTACGGTACGGATCTGGAGTTGCAATTAcgcc    >  1:7540606/1‑36 (MQ=255)
                                 cGTTTTACGGTACGGATCTGGAGTTGCAATTAcgcc    <  1:7431970/36‑1 (MQ=255)
                                 cGTTTTACGGTACGGATCTGGAGTTGCAATTAcgcc    <  1:1701305/36‑1 (MQ=255)
                                  gTTTTACGGTACGGATCTGGAGTTGCAATTAcgccg   <  1:3963432/36‑1 (MQ=255)
                                  gTTTTACGGTACGGATCTGGAGTTGCAATTAcgccg   >  1:7043775/1‑36 (MQ=255)
                                   ttttACGGTACGGATCTGGAGTTGCAATTAcgccgc  <  1:6641921/36‑1 (MQ=255)
                                   ttttACGGTACGGATCTGGAGTTGCAATTAcgccgc  <  1:3284816/36‑1 (MQ=255)
                                   ttttACGGTACGGAGCTGGAGTTGCAATTAcgccgc  >  1:7829444/1‑36 (MQ=255)
                                   |                                   
GCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGC  >  REL606/1929885‑1929955

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG