Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,893,550 +G 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.→G100.0% 115.1 42intergenic (+5/‑51)kup/insJ‑5potassium transporter/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  new base (20/21):  ref base (0/1):  total (20/22)

TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA.G.ACTGTACTGCACCCATTTTGTTGGACGAT  >  REL606/3893515‑3893580
                                    |                               
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCa                                  >  1:1452688/1‑36 (MQ=255)
tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCa                                  >  1:149672/1‑36 (MQ=255)
 aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCA.g                                <  2:1416990/36‑1 (MQ=255)
 aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCA.g                                <  1:1374618/36‑1 (MQ=255)
 aTCGAACTGGGTACTCAAGTCGAAATCTAACGCCA.g                                <  1:188005/36‑1 (MQ=255)
  tCGAACTGGGTACTCAGGTCGAAATCTAACGCCA.gg                               >  2:2564717/1‑35 (MQ=255)
   cGAACTGGGTACTCAGGTCGAAATCTAACGCCA.gga                              >  2:1036625/1‑34 (MQ=255)
   cGAACTGGGTACTCAGGTCGAAATCTAACGCCA.gga                              >  2:2594534/1‑34 (MQ=255)
   cGAACGGGGTACTCAGGTCGAAATCTAACGCCA.gga                              <  2:3299090/36‑3 (MQ=38)
    gAACTGGGTACTCAGGTTGAAATTTAACGTCA.ggac                             >  1:204224/1‑33 (MQ=16)
       cTGGGTACTCAGGTCGAAATCTAACGCCAGG.ACTGt                          <  2:3175880/36‑1 (MQ=37)
       cTGGGTACTCAGGTCGAAATCTAACGCCAGG.ACTGt                          <  1:2238065/36‑1 (MQ=37)
       cTGGGTACTCAGGTCGAAATCTAACGCCAGG.ACTGt                          >  1:985841/1‑36 (MQ=37)
        tGGGTAGTCAGGTCGAGATCTAGCGCCAGG.ACTGTa                         <  1:861215/36‑1 (MQ=16)
        tGGGTACGCAGGTCGAAATCTAACGCCAGG.ACTGTa                         <  2:1321021/36‑1 (MQ=25)
         gggTACTCAGGTCGAAATCTAACGCCAGG.ACTGTAc                        <  2:1917827/36‑1 (MQ=37)
          ggTACTCAGGTCGAAATCTAACGCCAGG.ACTGTACt                       <  2:2149596/36‑1 (MQ=37)
           gTACTCAGGTCGAAATCTAACGCCAGG.ACTGTACTg                      >  1:3490016/1‑36 (MQ=37)
            tACTCAGGTCGAAATCTAACGCCAGG.ACTGTACtac                     >  1:2535997/1‑34 (MQ=37)
            tACTCAGGTCGAAATCTAACGCCAGG.ACTGTACTGc                     <  2:2427016/36‑1 (MQ=37)
            tACTCAGGTCGAAATCTAACGCCAGG.ACTGTACTGc                     >  1:1497438/1‑36 (MQ=37)
              cgcAGATCGAAATCTAACGCCAGG.ACTGTACTGCAc                   <  2:1103444/34‑1 (MQ=21)
              cTCAGGTCGAAATCTAACGCCAGG.ACTGTACTGCAc                   <  2:2559903/36‑1 (MQ=37)
                cAGGTCGAAATCTAACGCCAGG.ACTGTACTGCAccc                 <  1:2006929/36‑1 (MQ=37)
                cAGGTCGAAATCTAACGCCAGG.ACTGTACTGCAccc                 <  2:2936310/36‑1 (MQ=37)
                 aGGTAGAAATCTAACGCCAGG.ACTGTACTGCACCCa                <  2:2808363/36‑1 (MQ=25)
                   gTCGAAATCTAACGCCAGG.ACTGTACTGCACCCAtt              <  2:1223936/36‑1 (MQ=37)
                   gTCGAAATCTAACGCCAGG.ACTGTACTGCACCCAtt              >  2:790420/1‑36 (MQ=37)
                   gTCGAAATCTAACACCAGG.CCTGTACTGCACCCAtt              >  1:833964/1‑36 (MQ=16)
                    tCGAATTCTAACGCCAGG.ACTGTACTGCACCCAttt             <  2:3011493/36‑1 (MQ=21)
                    tCGAAATCTAACGCCAGG.ACTGTACTGCACCCAttt             >  1:1765962/1‑36 (MQ=25)
                     cGAAATCTAACGCCAGG.ACTGTACTGCACCCAtttt            <  2:2889975/36‑1 (MQ=25)
                     cGAAATCTAACGCCAGG.ACTGAACTGCACCCAtttt            >  1:55673/1‑36 (MQ=21)
                      gAAATCTAACGCCAGG.ACTGTACTGCACCCATTttg           >  1:1095798/1‑36 (MQ=21)
                      gAAATCTAACGCCAGG.ACTGTACTGCACCCATTttg           >  1:3468000/1‑36 (MQ=21)
                      gAAATCTAACGCCAGG.ACTGTACTGCACCCATAttg           >  2:2343004/1‑36 (MQ=21)
                      aaaaTCTAACGCCAGG.ACTGTACTGCACCCATTttg           <  1:943748/35‑1 (MQ=21)
                       aaaTCTAACGCCAGG.ACTGTACTGCACCCATTttgt          >  2:1497339/1‑36 (MQ=21)
                        aaTCTAACGCCAGG.ACTGTACTGCACCCATTttgtt         <  2:3360392/36‑1 (MQ=18)
                        aaTCTAACGCCAGG.ACTGTACTGCACCCATTttgtt         >  1:2212070/1‑36 (MQ=18)
                        aaTCTAACGCCAGG.ACTGTACTGCACCCATTttgtt         <  2:188578/36‑1 (MQ=18)
                         aTCTAACGCCAGG.ACTGTACTGCACCCATTttgttg        <  2:1555118/36‑1 (MQ=18)
                          tCTAACGCCAGG.ACTGTACTGCACCCATTTTGTTgg       >  2:3640017/1‑36 (MQ=18)
                           cTAACGCCAGG.ACTGTACTGCACCCATTTTGTTGGa      >  1:3083527/1‑36 (MQ=17)
                           cTAACGCCAGG.ACTGTACTGCACCCATTTTGTTGGa      <  1:2951286/36‑1 (MQ=17)
                            tAACGCCAGG.ACTGTACTGCACCCATTTTGTTGGAc     >  1:1043061/1‑36 (MQ=17)
                            tAACGCCA.GTACTGTACTGCACCCATTTTGTTGGAc     <  2:318651/36‑1 (MQ=17)
                             aaCGCCAGG.ACTGTACTGCACCCATTTTGTTGGACg    >  1:339039/1‑36 (MQ=14)
                             aaCGCCAGG.ACTGTACTGCACCCATTTTGTTGGACg    >  2:799598/1‑36 (MQ=14)
                               cGCCAGG.ACTGTACTGCACCCATTTTGTTGGACGAt  >  1:261264/1‑36 (MQ=11)
                               cGCCAGG.ACTGTACTGCACCCATTTTGTTGGACGAt  >  1:2395945/1‑36 (MQ=11)
                               cGCCAGG.ACTGTACTGCACCCATTTTGTTGGACGAt  <  1:2325459/36‑1 (MQ=11)
                                    |                               
TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA.G.ACTGTACTGCACCCATTTTGTTGGACGAT  >  REL606/3893515‑3893580

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG