Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,893,550 +G 52.5% 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.→G52.5% 122.4 ‑3.9 0.059intergenic (+5/‑51)kup/insJ‑5potassium transporter/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  new base (10/21):  ref base (18/10):  total (28/31)
Fisher's exact test strand distribution p-value = 1.94e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTACTGCACCCATTTTGTTGGACGATGAAAT  >  REL606/3893515‑3893585
                                    |                                    
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCa                                       <  1:3959775/36‑1 (MQ=255)
 aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..g                                    <  1:6432068/36‑1 (MQ=255)
  tCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..gg                                   >  1:1169334/1‑35 (MQ=255)
  tCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..gg                                   <  1:566442/36‑2 (MQ=255)
   cGAACTGGGTACTCAGGTCGAAATCTAACGCCA..ggg                                  <  1:2284526/36‑3 (MQ=255)
   cGAACTGGGTACTCAGGTCGAAATCTAACGCCA..ggg                                  <  1:1568981/36‑3 (MQ=255)
   cGAACTGGGTACTCAGGTCGAAATCTAACGCCA..GAc                                  >  1:5062583/1‑36 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCCA..ggga                                 <  1:865649/36‑4 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCCA..GACt                                 >  1:4669058/1‑36 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCCA..GACt                                 >  1:5910371/1‑36 (MQ=255)
    gAACTGGGTACTCAGGGCGAAATCTAACGCCA..ggac                                 >  1:1373942/1‑33 (MQ=37)
     aaCTGGGTACTCAGGTCGAAATCTAACGCCA..GACTg                                >  1:4082048/1‑36 (MQ=255)
      aCTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGt                               >  1:5392904/1‑36 (MQ=255)
       cTGGGTACTCAGGTCGAAATCTAACGCCAG.GACTGt                               <  1:3921093/36‑1 (MQ=37)
       cTGGGTACTCAGGTCGAAATCTAACGCCAG.GACTGt                               <  1:3345186/36‑1 (MQ=37)
        tGGGTACTCAGGTCGAAATCTAACGCCAG.GACTGTa                              >  1:1488490/1‑36 (MQ=37)
        tGGGTACTCAGGTCGAAATCTAACGCCAG.GACTGTa                              >  1:5857391/1‑36 (MQ=37)
        tGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTAc                             <  1:174289/36‑1 (MQ=255)
          ggTCCTCAGGTCGAAAACTAACGCCAG.GACTGTACt                            <  1:5164588/36‑1 (MQ=16)
          ggTACTCAGGTCGAAATCTAACGCCAG.GACTGTACt                            <  1:7885861/36‑1 (MQ=37)
          ggTACTCAGGTCGAAATCTAACGCCAG.GACTGTACt                            <  1:1435460/36‑1 (MQ=37)
          ggTACTCAGGTCGAAATCTAACGCCA..GACTGTACTg                           >  1:1229646/1‑36 (MQ=255)
          ggTACTCAGGTCGAAATCTAACGCCA..GACTGTACTg                           <  1:7623076/36‑1 (MQ=255)
           gTACTCAGGACGAAATCTAACGCCAG.GACTGTACTg                           <  1:3235028/36‑1 (MQ=25)
            tACTCAGGTCGAAATCTAACGCCAG.GACTGTACTGc                          <  1:5682266/36‑1 (MQ=37)
              cTCAGGTCGAAATCTAACGCCAG.GACTGTACTGCAc                        <  1:3132448/36‑1 (MQ=37)
               tCAGGTCGAAATCTAACGCCAGGGACTGTACTGCAc                        <  1:1220785/36‑1 (MQ=25)
               tCAGGTCGAAATCTAACGCCAGGGACTGTACTGCAc                        <  1:4426115/36‑1 (MQ=25)
               tCAGGTCGAAATCTAACGCCA..GACTGTACTGCAccc                      >  1:7750356/1‑36 (MQ=255)
               tCAGGTCGAAATCTAACGCCA..GACTGTACTGCAccc                      >  1:5687016/1‑36 (MQ=255)
                cAGGTCGAAATCTAACGCCAG.GACTGTACTGCAccc                      >  1:1338767/1‑36 (MQ=37)
                cAGGTCGAAATCTAACGCCA..GACTGTACTGCACCCa                     <  1:3541793/36‑1 (MQ=255)
                cAGGTCGAAATCTAACGCCA..GACTGTACTGCACCCa                     <  1:3884614/36‑1 (MQ=255)
                cAGGGCGAAATCTAACGCCAG.GACTGTACTGCAccc                      >  1:4221981/1‑36 (MQ=25)
                 aGGTCGAAATCTAACGCCA..GACTGTACTGCACCCAt                    >  1:1668159/1‑36 (MQ=255)
                  ggTCGAAATCTAACGCCA..GACTGTACTGCACCCAtt                   <  1:1436699/36‑1 (MQ=255)
                  ggTCGAAATCTAACGCCA..GACTGTACTGCACCCAtt                   >  1:464716/1‑36 (MQ=255)
                    tCGAAATCTAACGCCAGGGACTGTACTGCACCCAtt                   <  1:3501556/36‑1 (MQ=21)
                    tCGAAATCTAACGCCAG.GACTGTACTGCACCCAttt                  <  1:2539491/36‑1 (MQ=25)
                    tCGAAATCTAACGCCAG.GACTGTACTGCACCCAttt                  <  1:597603/36‑1 (MQ=25)
                    tCGAAATCTAACGCCA..GACTGTACTGCACCCAtttt                 <  1:7091624/36‑1 (MQ=255)
                      gAAATCTAACGCCAGGGACTGTACTGCACCCAtttt                 <  1:4414649/36‑1 (MQ=18)
                      gAAATCTAACGCCA..GACTGTACTGCACCCATTttgt               >  1:5988658/1‑36 (MQ=255)
                        aaTCTAACGCCAGGGACTGTACTGCACCCATTttgt               >  1:4000488/1‑36 (MQ=17)
                        aaTCTAACGCCAG.GACTGTACTGCACCCATTttgtt              <  1:1601477/36‑1 (MQ=18)
                        aaTCTAACGCCA..GACTGTACTGCACCCATTttgttg             <  1:1638731/36‑1 (MQ=255)
                        aaTCTAACGCCA..GACAGGACTGCACACATTttgttg             >  1:2066881/1‑36 (MQ=25)
                         aTCTAACGCCA..GACTGTACTGCACCCATTTTGTTgg            >  1:1193553/1‑36 (MQ=255)
                         aTCTAACGCCA..GACTGTACTGCACCCATTTTGTTgg            >  1:6129105/1‑36 (MQ=255)
                         aTCTAACGCCA..GACTGTACTGCACCCATTTTGTTgg            <  1:2037049/36‑1 (MQ=255)
                          tCTAACGCCAGGGACTGTACTGCACCCATTttgttg             <  1:5294812/36‑1 (MQ=14)
                          tCTAACGCCAG.GACTGTACTGCACCCATTTTGTTgg            >  1:3797144/1‑36 (MQ=18)
                           cTAACGCCAGGGACTGTACTGCACCCATTTTGTTgg            <  1:57519/36‑1 (MQ=14)
                           cTAACGCCAGGGACTGTACTGCACCCATTTTGTTgg            >  1:6925174/1‑36 (MQ=14)
                           cTAACGCCA..GACTGTACTGCACCCATTTTGTTGGAc          <  1:1695972/36‑1 (MQ=255)
                            tAACGCCAG.GACTGTACTGCACCCATTTTGTTgggg          >  1:2175791/1‑34 (MQ=17)
                            tAACGCCAG.GACTGTACTGCACCCATTTTGTTGGAc          <  1:3352411/36‑1 (MQ=17)
                             caCGCCAG.GACTGTACTGCACCCATTTTGTTGGACg         <  1:6278421/35‑1 (MQ=14)
                             aaCGCCAG.GACTGTACTGCACCCATTTTGTTGGACg         >  1:4151082/1‑36 (MQ=14)
                             aaCGCCAG.GACTGTACTGCACCCATTTTGTTGGACg         <  1:3975722/36‑1 (MQ=14)
                             aaCGCCA..GACTGTACTGCACCCATTTTGTTGGACGa        >  1:296834/1‑36 (MQ=255)
                              aCGCCAG.GACTGTACTGCACCCATTTTGTTGGACGa        >  1:5554334/1‑36 (MQ=14)
                               cGCCAG.GACTGTACTGCACCCATTTTGTTGGACGAt       <  1:2613995/36‑1 (MQ=11)
                                gCCA..GACTGTACTGCACCCATTTTGTTGGACGATGa     >  1:5497695/1‑36 (MQ=255)
                                gCCA..GACTGTACTGCACCCAGTTTGTTGGACGATGa     >  1:2400344/1‑36 (MQ=255)
                                 ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa    <  1:4751092/36‑1 (MQ=255)
                                 ccA..GACTGTACTGCACCCATTTTGTTGGACGATGaa    >  1:7238705/1‑36 (MQ=255)
                                   a..GACTGTACTGCACCCATTTTGTTGGACGATGAAAt  >  1:1978697/1‑36 (MQ=38)
                                   a..GACTGTACTGCACCCATTTTGTTGGACGATGAAAt  <  1:7419286/36‑1 (MQ=38)
                                    |                                    
TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA..GACTGTACTGCACCCATTTTGTTGGACGATGAAAT  >  REL606/3893515‑3893585

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