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
*	REL606	3,893,550	1	.→G	100.0%	101.5	41	intergenic (+5/‑51)	kup/insJ‑5	potassium transporter/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  new base (20/20):  ref base (0/1):  total (20/21)

TATCGAACTGGGTACTCAGGTCGAAATCTAACGCCA.G.ACTGTACTGCACCCATTTTGTTGGACGAT  >  REL606/3893515‑3893580                                     |                                tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCa                                  >  1:1452688/1‑36 (MQ=255) tATCGAACTGGGTACTCAGGTCGAAATCTAACGCCa                                  >  1:149672/1‑36 (MQ=255)  aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCA.g                                <  1:1374618/36‑1 (MQ=255)  aTCGAACTGGGTACTCAGGTCGAAATCTAACGCCA.g                                <  1:5156200/36‑1 (MQ=255)  aTCGAACTGGGTACTCAAGTCGAAATCTAACGCCA.g                                <  1:188005/36‑1 (MQ=255)   tCGAACTGGGTACTCAGGTCGAAATCTAACGCCA.gg                               >  1:6303927/1‑35 (MQ=255)    cGAACTGGGTACTCAGGTCGAAATCTAACGCCA.gga                              >  1:6333744/1‑34 (MQ=255)    cGAACTGGGTACTCAGGTCGAAATCTAACGCCA.gga                              >  1:4775835/1‑34 (MQ=255)    cGAACGGGGTACTCAGGTCGAAATCTAACGCCA.gga                              <  1:7038300/36‑3 (MQ=37)     gAACTGGGTACTCAGGTTGAAATTTAACGTCA.ggac                             >  1:204224/1‑33 (MQ=18)        cTGGGTACTCAGGTCGAAATCTAACGCCAGG.ACTGt                          >  1:985841/1‑36 (MQ=25)        cTGGGTACTCAGGTCGAAATCTAACGCCAGG.ACTGt                          <  1:6915090/36‑1 (MQ=25)        cTGGGTACTCAGGTCGAAATCTAACGCCAGG.ACTGt                          <  1:2238065/36‑1 (MQ=25)         tGGGTACGCAGGTCGAAATCTAACGCCAGG.ACTGTa                         <  1:5060231/36‑1 (MQ=21)          gggTACTCAGGTCGAAATCTAACGCCAGG.ACTGTAc                        <  1:5657037/36‑1 (MQ=25)           ggTACTCAGGTCGAAATCTAACGCCAGG.ACTGTACt                       <  1:5888806/36‑1 (MQ=25)            gTACTCAGGTCGAAATCTAACGCCAGG.ACTGTACTg                      >  1:3490016/1‑36 (MQ=25)             tACTCAGGTCGAAATCTAACGCCAGG.ACTGTACtac                     >  1:2535997/1‑34 (MQ=21)             tACTCAGGTCGAAATCTAACGCCAGG.ACTGTACTGc                     >  1:1497438/1‑36 (MQ=25)             tACTCAGGTCGAAATCTAACGCCAGG.ACTGTACTGc                     <  1:6166226/36‑1 (MQ=25)               cgcAGATCGAAATCTAACGCCAGG.ACTGTACTGCAc                   <  1:4842654/34‑1 (MQ=18)               cTCAGGTCGAAATCTAACGCCAGG.ACTGTACTGCAc                   <  1:6299113/36‑1 (MQ=25)                 cAGGTCGAAATCTAACGCCAGG.ACTGTACTGCAccc                 <  1:6675520/36‑1 (MQ=25)                 cAGGTCGAAATCTAACGCCAGG.ACTGTACTGCAccc                 <  1:2006929/36‑1 (MQ=25)                  aGGTAGAAATCTAACGCCAGG.ACTGTACTGCACCCa                <  1:6547573/36‑1 (MQ=21)                    gTCGAAATCTAACGCCAGG.ACTGTACTGCACCCAtt              >  1:4529630/1‑36 (MQ=21)                    gTCGAAATCTAACGCCAGG.ACTGTACTGCACCCAtt              <  1:4963146/36‑1 (MQ=21)                    gTCGAAATCTAACACCAGG.CCTGTACTGCACCCAtt              >  1:833964/1‑36 (MQ=17)                     tCGAATTCTAACGCCAGG.ACTGTACTGCACCCAttt             <  1:6750703/36‑1 (MQ=18)                     tCGAAATCTAACGCCAGG.ACTGTACTGCACCCAttt             >  1:1765962/1‑36 (MQ=21)                      cGAAATCTAACGCCAGG.ACTGTACTGCACCCAtttt            <  1:6629185/36‑1 (MQ=21)                      cGAAATCTAACGCCAGG.ACTGAACTGCACCCAtttt            >  1:55673/1‑36 (MQ=18)                       gAAATCTAACGCCAGG.ACTGTACTGCACCCATTttg           >  1:1095798/1‑36 (MQ=18)                       gAAATCTAACGCCAGG.ACTGTACTGCACCCATTttg           >  1:3468000/1‑36 (MQ=18)                       gAAATCTAACGCCAGG.ACTGTACTGCACCCATAttg           >  1:6082214/1‑36 (MQ=18)                       aaaaTCTAACGCCAGG.ACTGTACTGCACCCATTttg           <  1:943748/35‑1 (MQ=18)                        aaaTCTAACGCCAGG.ACTGTACTGCACCCATTttgt          >  1:5236549/1‑36 (MQ=18)                         aaTCTAACGCCAGG.ACTGTACTGCACCCATTttgtt         >  1:2212070/1‑36 (MQ=18)                         aaTCTAACGCCAGG.ACTGTACTGCACCCATTttgtt         <  1:7099602/36‑1 (MQ=18)                         aaTCTAACGCCAGG.ACTGTACTGCACCCATTttgtt         <  1:3927788/36‑1 (MQ=18)                          aTCTAACGCCAGG.ACTGTACTGCACCCATTttgttg        <  1:5294328/36‑1 (MQ=17)                           tCTAACGCCAGG.ACTGTACTGCACCCATTTTGTTgg       >  1:7379227/1‑36 (MQ=17)                            cTAACGCCAGG.ACTGTACTGCACCCATTTTGTTGGa      <  1:2951286/36‑1 (MQ=17)                            cTAACGCCAGG.ACTGTACTGCACCCATTTTGTTGGa      >  1:3083527/1‑36 (MQ=17)                             tAACGCCAGG.ACTGTACTGCACCCATTTTGTTGGAc     >  1:1043061/1‑36 (MQ=14)                             tAACGCCA.GTACTGTACTGCACCCATTTTGTTGGAc     <  1:4057861/36‑1 (MQ=14)                              aaCGCCAGG.ACTGTACTGCACCCATTTTGTTGGACg    >  1:4538808/1‑36 (MQ=14)                              aaCGCCAGG.ACTGTACTGCACCCATTTTGTTGGACg    >  1:339039/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  < 21 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG