Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,893,550 +G 2.9% 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.→G 2.9% 3.4 0.0 ‑0.068intergenic (+5/‑51)kup/insJ‑5potassium transporter/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  new base (1/1):  ref base (36/30):  total (37/31)
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.45e-01

TATCGAACTGGGTACTCAGGTCGAAATCTAACGCC..A.GACTGTACTGCACCCATTTTGTTGGACGATGAAAT  >  REL606/3893515‑3893585
                                      |                                   
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCC..a                                      <  1:1399784/36‑1 (MQ=255)
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCC..a                                      <  1:7858509/36‑1 (MQ=255)
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCC..a                                      >  1:4229704/1‑36 (MQ=255)
tATCGAACTGGGTACTAAGGTGGAAATCTAACGCC..a                                      <  1:2905189/36‑1 (MQ=37)
  tCGAACTGGGTACTCAGGTCGAAATCTAACGCC..A.Ga                                   >  1:1958164/1‑36 (MQ=255)
   cGAACTGGGTACTCAGGTCGAAATCTAACGCC..A.GAc                                  <  1:6381399/36‑1 (MQ=255)
   cGAAATGGGTACTCAGGTCGAAATCTAACGCC..A.GAc                                  <  1:995003/36‑1 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCC..A.GACt                                 <  1:3955179/36‑1 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCC..A.GACt                                 >  1:2630911/1‑36 (MQ=255)
    gAACTGGGTACTCAGGTCGAAATCTAACGCC..A.GACt                                 >  1:6626071/1‑36 (MQ=255)
     aaCTGGGTACTCAGGTCGAAATCTAACGCC..A.GACTg                                >  1:2296416/1‑36 (MQ=255)
      aCTGGGTACTCAGGTCGAAATCTAACGCC..A.GTCTGt                               >  1:7015075/1‑36 (MQ=255)
      aCTGGGTACTCAGGTCGAAATCTAACGCC..A.GACTGt                               >  1:4615689/1‑36 (MQ=255)
       cTGGGTACTCAGGTCGAAATCTAACGCC..A.GACTGTa                              >  1:566343/1‑36 (MQ=255)
       cTGGGTACTCAGGTCGAAATCTAACGCC..A.GACTGTa                              >  1:195489/1‑36 (MQ=255)
        tGGGTACTCAGGTCGAAATCTAACGCC..A.GACTGTAc                             >  1:7931614/1‑36 (MQ=255)
         gggTACTCAGGTCGAAATCTAACGCC..A.GACTGTACt                            <  1:6771291/36‑1 (MQ=255)
         gggTACTCAGGTCGAAATCTAACGCC..A.GACTGTACt                            >  1:4274898/1‑36 (MQ=255)
          ggTACTCAGGTCGAAATCTAACGCC..A.GACTGTACTg                           >  1:683973/1‑36 (MQ=255)
           gTACTCAGGTCGAAATCTAACGCC..A.GACTGTACTGc                          <  1:5637108/36‑1 (MQ=255)
            tACTCAGGTCGAAATCTAACGCC..A.GACTGTACTGCa                         <  1:4845119/36‑1 (MQ=255)
            tACTCAGGTCGAAATCTAACGCC..A.GACTGTACTGCa                         <  1:2899220/36‑1 (MQ=255)
            tACTCAGGTCGAAATCTAACGCC..A.GACTGTACTGCa                         >  1:7536128/1‑36 (MQ=255)
            tACTCAGGTCGAAATCTAACGCC..A.GACTGTACTGCa                         <  1:3531016/36‑1 (MQ=255)
              cTCAGGTCGAAATCTAACGCC..AGGACTGTACTGCAc                        >  1:4180558/1‑36 (MQ=37)
              cTCAGGTCGAAATCTAACGCC..A.GACTGTACTGCAcc                       >  1:7228953/1‑36 (MQ=255)
              cTCAGGTCGAAATCTAACGCC..A.GACTGTACTGCAcc                       <  1:4551900/36‑1 (MQ=255)
              cTCAGGTCGAAATCTAACGCC..A.GACTGTACTGCAcc                       <  1:4227202/36‑1 (MQ=255)
               tCAGGTCGAAATCTAACGCC..A.GACTGTACTGCAccc                      >  1:6493059/1‑36 (MQ=255)
               tCAGGTCGAAATCTAACGCC..A.GACTGTACTGCAccc                      <  1:911864/36‑1 (MQ=255)
               tCAGGTCGAAATCTAACGCC..A.GACTGTACTGCAccc                      >  1:857247/1‑36 (MQ=255)
                cAGGTCGAAATCTAACGCC..A.GACTGTACTGCACCCa                     <  1:557977/36‑1 (MQ=255)
                cAGGTCGAAATCTAACGCC..A.GACTGTACTGCACCCa                     >  1:1499139/1‑36 (MQ=255)
                 aGGTCGAAATCTAACGCC..A.GACTGTACTGCACCCAt                    >  1:7260045/1‑36 (MQ=255)
                  ggTCGAAATCTAACGCC..A.GACTGTACTGCACCCAtt                   >  1:6636625/1‑36 (MQ=255)
                  ggTCGAAATCTAACGCC..A.GACTGTACTGCACCCAtt                   >  1:2554023/1‑36 (MQ=255)
                  ggTCGAAATCTAACGCC..A.GACTGTAATGCACCCAtt                   >  1:5643756/1‑36 (MQ=255)
                   gTCGAAATCTAACGCC..A.GACTGTACTGCACCCAttt                  <  1:4135152/36‑1 (MQ=255)
                    tCGAAATCTAACGCC..A.GACTGTACTGCACCCAtttt                 <  1:1279144/36‑1 (MQ=255)
                     cGAAATCTAACGCC..A.GACTGTACTGCACCCATTttg                <  1:5255568/36‑1 (MQ=255)
                      gAAATCTAACGCC..A.GACTGTACTGCACCCATTttgt               >  1:6547010/1‑36 (MQ=255)
                       aaaTCTAACGCC..A.GACTGTACTGCACCCATTttgtt              <  1:4035746/36‑1 (MQ=255)
                       aaaTCTAACGCC..A.GACTGTACTGCACCCATTttgtt              >  1:892995/1‑36 (MQ=255)
                       aaaTCTAACGCC..A.GACTGTACTGCACCCATTttgtt              >  1:235420/1‑36 (MQ=255)
                        aaTCTAACTCCAGA.GACTGTACTGCACCCATTttgt               <  1:3466641/36‑1 (MQ=9)
                        aaTCTAACGCC..A.GACTGTACTGCACCCATTttgttg             >  1:1128868/1‑36 (MQ=255)
                        aaTCTAACGCC..A.GACTGTACTGCACCCATTttgttg             >  1:5608419/1‑36 (MQ=255)
                        aaTCTAACGCC..A.GACGGTACTGCACCCATTttgttg             >  1:3317989/1‑36 (MQ=255)
                         aTCTAACGCC..A.GACTGTACTGCACCCATTTTGTTgg            <  1:936209/36‑1 (MQ=255)
                          tCTAACGCC..A.GACTGTACTGCACCCATTTTGTTGGa           <  1:3035435/36‑1 (MQ=255)
                          tCTAACGCC..A.GACTGTACTGCACCCATTTTGTTGGa           <  1:2422605/36‑1 (MQ=255)
                          tCTAACGCC..A.GACTGTACTGCACCCATTTTGTTGGa           <  1:3358139/36‑1 (MQ=255)
                           cTACCGCC..A.GCCTGTACTGCTCCCATTTTGTTGGAc          <  1:1687327/36‑1 (MQ=14)
                           cTAACGCC..A.GACTGTACTGCACCCATTTTGTTGGAc          <  1:5900857/36‑1 (MQ=255)
                           cTAACGCC..A.GACTGTACTGCACCCATTTTGTTGGAc          <  1:7114689/36‑1 (MQ=255)
                            tAACGCC..AGGACTGTACTGCACCCATTTTGTTGGAc          <  1:2012501/36‑1 (MQ=17)
                            tAACGCC..A.GACTGTACTGCACCCATTTTGTTGGACg         <  1:410175/36‑1 (MQ=255)
                            tAACGCC..A.GACTGTACTGCACCCATTTTGTTGGACg         >  1:1631237/1‑36 (MQ=255)
                              aCGCC..A.GACTGTACTGCACCCATTTTGTTGGACGAt       <  1:3230473/36‑1 (MQ=255)
                              aCGCC..A.GACTGTACTGCACCCATTTTGTTGGACGAt       <  1:3075874/36‑1 (MQ=255)
                              aCGCC..A.GACTGTACTGCACCCAGTTTGTTGGACGAt       >  1:2546061/1‑36 (MQ=255)
                               cGCC..A.GACTGTACTGCACCCATTTTGTTGGACGATg      <  1:3459193/36‑1 (MQ=255)
                               cGCC..A.GACTGTACTGCACCCATTTTGTTGGACGATg      <  1:3443427/36‑1 (MQ=255)
                                gCC..A.GACTGTACTGCCCCCATTTTGTTGGACGATGa     <  1:3925552/36‑1 (MQ=255)
                                gCC..A.GACTGTACTGCACCCATTTTGTTGGACGATGa     >  1:232430/1‑36 (MQ=255)
                                gCC..A.GACTGTACTGCACCCATTTTGTTGGACGATGa     >  1:3816814/1‑36 (MQ=255)
                                gCC..A.GACTGTACTGCACCCATTTTGTTGGACGATGa     >  1:7514632/1‑36 (MQ=255)
                                gCC..A.GACTGTACTGCACCCATTTTGTTGGACGATGa     >  1:1457962/1‑36 (MQ=255)
                                 cc..A.GACTGTACTGCACCCATTTTGTTGGACGATGaa    >  1:1572307/1‑36 (MQ=255)
                                 cc..A.GACTGTACTGCACCCATTTTGTTGGACGATGaa    >  1:4827012/1‑36 (MQ=255)
                                  c..A.GACTGTACTGCACCCATTTTGTTGGACGATGaaa   >  1:4326197/1‑36 (MQ=40)
                                  c..A.GACTGTACTGCACCCATTTTGTTGGACGATGaaa   >  1:3692102/1‑36 (MQ=40)
                                     a.GACTGTACTGCACCCATTTTGTTGGACGATGAAAt  <  1:7890095/36‑1 (MQ=38)
                                     a.GACTGTACTGCACCCATTTTGTTGGACGATGAAAt  <  1:6117298/36‑1 (MQ=38)
                                      |                                   
TATCGAACTGGGTACTCAGGTCGAAATCTAACGCC..A.GACTGTACTGCACCCATTTTGTTGGACGATGAAAT  >  REL606/3893515‑3893585

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG