Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,922,670 +G 6.7% coding (681/786 nt) znuB → high‑affinity zinc transporter membrane component

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,922,6701.→G 6.7% 4.1 ‑4.3 ‑0.330coding (681/786 nt)znuBhigh‑affinity zinc transporter membrane component
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (2/26):  total (4/26)
Fisher's exact test strand distribution p-value = 1.38e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.62e-01

TCGCTGTTTTGGTGGGGATGGTGGCAGTGACTGGC.GGTTTAACCTTTTCCGCATTTTACGATACGCCGGC  >  REL606/1922636‑1922705
                                   |                                   
tCGCTGTTTTGGTGGGGATGGTGGCAGTGACTGGC.g                                    <  1:5925219/36‑1 (MQ=255)
tCGCTGTTTTGGTGGGGATGGTGGCAGTGACTGGC.g                                    >  1:5640156/1‑36 (MQ=255)
tCGCTGTTTTGGTGGGGATGGTGGCAGTGACTGGC.g                                    >  1:4613628/1‑36 (MQ=255)
 atcTGTTTTGGTGGGGATGGTGGCAGTGACTGGC.gg                                   <  1:3836756/34‑1 (MQ=255)
  gCTGTTTTGGTGGGGATGGTGGCAGTGACTGGC.ggg                                  >  1:4774359/1‑35 (MQ=255)
  gCTGTTTTGGTGGGGATGGTGGCAGTGACTGGC.GGt                                  >  1:1524402/1‑36 (MQ=255)
   gtattttGGTGGGGATGGTGGCAGTGACTGGC.GGtt                                 <  1:759601/33‑1 (MQ=255)
   cTGTTTTGGTGGGGATGGTGGCAGTGACTGGC.gggt                                 >  1:5068360/1‑34 (MQ=255)
    tGTTTTGGTGGGGATGGTGGCAGTGACTGGC.GGttt                                <  1:88261/36‑1 (MQ=255)
     gTTTTGGGGGGGATGGTGGCAGTGACTGGC.GGTTTa                               <  1:6968173/36‑1 (MQ=255)
      ttttGGTGGGGATGGTGGCAGTGACTGGC.GGttaac                              <  1:6830980/36‑4 (MQ=255)
      ttttGGTGGGGATGGTGGCAGTGACTGGC.GGTTTaa                              <  1:6681362/36‑1 (MQ=255)
       tttGTTGGGGATGGTGGCAGTGACTGGC.GGTTTAAc                             <  1:2044333/36‑1 (MQ=255)
         tggtggGGATGGTGGCAGTGACTGGC.GGTTTAACCt                           <  1:2773445/36‑1 (MQ=255)
          ggtggGGATGGTGGCAGTGACTGGC.GGTTTAACCtt                          <  1:7071175/36‑1 (MQ=255)
          ggtggGGATGGTGGCAGTGACTGGC.GGGTTTAGCtt                          >  1:860168/1‑36 (MQ=25)
           gtggGGATGGTGGCAGTGACTGGC.GGTTTAACCttt                         <  1:3541210/36‑1 (MQ=255)
              gggATGGTGGCAGTGACTGGC.GGTTTAACCTTTTcc                      <  1:2063540/36‑1 (MQ=255)
              gggATGGTGGCAGTGACTGGC.GGTTTAACCTTTTcc                      <  1:6667316/36‑1 (MQ=255)
               ggATGGTGGCAGTGACTGGC.GGTTTAACCTTTTCCg                     <  1:4815918/36‑1 (MQ=255)
                gATGGTGGCAGTGACTGGC.GGTTTAACCTTTTCCGc                    <  1:6516034/36‑1 (MQ=255)
                  tggtggCAGTGACTGGC.GGTTTAACCTTTTCCGCAt                  <  1:4502335/36‑1 (MQ=255)
                   ggtggCAGTGACTGGC.GGTTTAACCTTTTCCGCAtt                 <  1:1513069/36‑1 (MQ=255)
                      ggCAGTGACTGGCGGGTTTAACTTTTTCCCCCtttt               >  1:1027491/1‑36 (MQ=16)
                      ggCAGTGACTGGC.GGTTTAACCTTTTCCGCATTTTa              <  1:5648437/36‑1 (MQ=255)
                      ggCAGTGACTGGC.GGTTTAACCTTTTCCGCATTTTa              <  1:5151595/36‑1 (MQ=255)
                        cAGTGACTGGCGGGTTTAACCTTTTTCCCTtttttt             >  1:8129978/1‑34 (MQ=12)
                           tGACTGGC.GGTTTAACCTTTTCCGCATTTTACGATa         <  1:7861044/36‑1 (MQ=255)
                            gACTGGT.GGTTTAACCTTTTCCGCATTTTACGATAc        <  1:948594/36‑1 (MQ=255)
                            gACTGGC.GGTTTAACCTTTTCCGCATTTTACGATAc        <  1:4088241/36‑1 (MQ=255)
                            gACTGGC.GGTTTAACCTTTTCCGCATTTTACGATAc        <  1:3565954/36‑1 (MQ=255)
                             aCTGGC.GGTTTAACCTTTTCCGCATTTTACGATACg       <  1:1227475/36‑1 (MQ=255)
                              cTGGC.GGTTTAACCTTTTCCGCATTTTACGATACGc      <  1:2815539/36‑1 (MQ=255)
                               tGGC.GGTTTAACCTTTTCCGCATTTTACGATACGcc     <  1:5130716/36‑1 (MQ=255)
                                ggC.GGTTTAACCTTTTCCGCATTTTACGATACGCCg    <  1:1049981/36‑1 (MQ=255)
                                 gC.GGTTTAACCTTTTCCGCATTTTACGATACGCCgg   <  1:7632461/36‑1 (MQ=255)
                                  c.GGTTTAACCTTTTCCGCATTTTACGATACGCCGGc  >  1:3108088/1‑36 (MQ=255)
                                   |                                   
TCGCTGTTTTGGTGGGGATGGTGGCAGTGACTGGC.GGTTTAACCTTTTCCGCATTTTACGATACGCCGGC  >  REL606/1922636‑1922705

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG