Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,893,550 +G 89.3% intergenic (+5/‑51) kup → / → insJ‑5 potassium transporter/IS150 hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,893,5501.→G89.3% 6.4 0.0 ‑0.328intergenic (+5/‑51)kup/insJ‑5potassium transporter/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  new base (12/13):  ref base (2/1):  total (14/14)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.49e-01

TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..G.ACTGTACTGCACCCATTTTGTTGGACGATGAA  >  REL606/3893515‑3893583
                                    |                                   
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCa                                      <  1:5891403/36‑1 (MQ=255)
 aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..g                                   <  1:3586645/36‑1 (MQ=255)
 aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..g                                   >  1:1875698/1‑36 (MQ=255)
 aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..g                                   >  1:1159314/1‑36 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCCA..ggac                                <  1:5594422/36‑4 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCCA..ggac                                >  1:5534411/1‑33 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCCA..ggac                                >  1:1189449/1‑33 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCCA..ggac                                >  1:3914195/1‑33 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCCA..ggac                                >  1:3362433/1‑33 (MQ=255)
       cTGGGTACTCAGGTCGAAATCTAACGCCAG.G.ACTGt                             <  1:3320652/36‑1 (MQ=37)
         gggTACTCAGGTCGAAATCTAACGCCAG.G.ACTGTAc                           >  1:4095428/1‑36 (MQ=37)
          ggTACTCAGGTCGAAATCTAACGCCAG.G.ACTGTACt                          <  1:1503102/36‑1 (MQ=37)
             aCTCAGGGCGAAATCTAACGCCAG.G.ACTGTACTGCa                       >  1:3093069/1‑36 (MQ=25)
              cTCAGGTCGAACTCTAACGCCAG.G.ACTGTACTGCAc                      <  1:1265000/36‑1 (MQ=25)
               tCAGGTCGAAATCTAACGCCAG.G.ACTGTACTGCAcc                     <  1:4576341/36‑1 (MQ=37)
               tCAGGTCGAAATCTAACGCCAG.G.ACTGTACTGCAcc                     <  1:2311700/36‑1 (MQ=37)
               tCAGGTCGAAATCTAACGCCAG.G.ACTGTACTGCAcc                     <  1:3504491/36‑1 (MQ=37)
                 aGGTCGAAATCTAACGCCAG.G.ACTGTACTGCACCCa                   <  1:3520246/36‑1 (MQ=37)
                  ggTCGAAGTCTGACGCCA..GTACTGTACTGCACCCAt                  <  1:1114119/36‑1 (MQ=16)
                  ggTCGAAATCTAACGCCAG.G.ACTGTACTGCACCCAt                  >  1:667394/1‑36 (MQ=37)
                    tCGAAATCTAACGCCAG.G.ACTGTACTGCACCCAttt                <  1:5284910/36‑1 (MQ=25)
                     cGAAATCTAACGCCAGGG.ACTGTACTGCACCCAttt                >  1:1365924/1‑36 (MQ=18)
                     cGAAATCTAACGCCAG.G.ACTGTACTGCACCCAtttt               <  1:299801/36‑1 (MQ=25)
                     cGAAATCTAACGCCAG.G.ACTGTACTGCACCCAtttt               >  1:645821/1‑36 (MQ=25)
                      gAAATCTAACGCCAG.G.ACTGTACTGCACCCATTttg              >  1:2173644/1‑36 (MQ=21)
                      gAAATCTAACGCCAG.G.ACTGTACTGCACCCATTttg              >  1:174083/1‑36 (MQ=21)
                      gAAATCTAACGCCA..GTACTGTACTGCACCCAtttcg              >  1:715336/1‑34 (MQ=21)
                         aTCTAACGCCAG.G.ACTGTACTGCACCCATTttgttg           >  1:1260686/1‑36 (MQ=18)
                         aTCTAACGCCAG.G.ACTGTACTGCACCCATTttgttg           >  1:2770401/1‑36 (MQ=18)
                         aTCTAACGCCAG.G.ACTGTACTGCACCCATTtcgtgg           >  1:452704/1‑34 (MQ=18)
                          tttAACGCCAG.G.ACTGTACTGCACCCATTTTGTTgg          <  1:4016714/34‑1 (MQ=17)
                          tCTAACGCCAG.G.ACTGTACTGCACCCATTTTGTTgg          >  1:585928/1‑36 (MQ=18)
                          tCTAACGCCAG.G.ACTGTACTGCACCCATTTTGTTgg          <  1:2801576/36‑1 (MQ=18)
                            tAACGCCAG.G.ACTGTACTGCACCCATTTTGTTGGAc        <  1:5379827/36‑1 (MQ=17)
                            tAACGCCAG.G.ACTGTACTGCACCCATTTTGTTGGAc        <  1:550156/36‑1 (MQ=17)
                              aCGCCAG.G.ACTGTACTGCACCCATTTTGTTGGACGa      >  1:1325916/1‑36 (MQ=14)
                                 acA..G.ACTGTACTGCACCCATTTTGTTGGACGATGaa  >  1:5371539/2‑36 (MQ=255)
                                    |                                   
TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..G.ACTGTACTGCACCCATTTTGTTGGACGATGAA  >  REL606/3893515‑3893583

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG