Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,201,911	C→T	A217T (GCG→ACG)	iclR	DNA‑binding transcriptional repressor

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	REL606	4,201,911	0	C→T	100.0%	175.8	59	A217T (GCG→ACG)	iclR	DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  new base (29/28):  ref base (0/1):  total (29/30)

AAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTT  >  REL606/4201878‑4201943                                  |                                 aaaTACGCGCGGCGAGGCAACGTAGCCCCAGCGTAt                                <  1:1178323/36‑1 (MQ=18) aaaTACACGCTGCAAGGCAACGTAGCCCCAGCGTAt                                <  1:2252761/36‑1 (MQ=255) aaaTACACGCTGCAAGGCAACGTAGCCCCAGCGTAt                                >  1:5630703/1‑36 (MQ=255) aaaAACACGCTGCAAGGCAACGTAGCCCCAGCGTAt                                <  1:1738620/36‑1 (MQ=25)  aaTACACGCTGCAAGGCAACGTAGCCCCCGCGTATg                               >  1:3964999/1‑36 (MQ=25)  aaTACACGCTGCAAGGCAACGTAGCCCCAGCGTATg                               >  1:6434321/1‑36 (MQ=255)   aTACACGCTGCAAGGCAACGTAGCCCCAGCGTATGt                              >  1:2156437/1‑36 (MQ=255)   aTACACGCTGCAAGGCAACGTAGCCCCAGCGTATGt                              >  1:4179876/1‑36 (MQ=255)    tACACGCTGCAAGGCAACGTAGCCCCAGCGTATGtt                             <  1:4127550/36‑1 (MQ=255)     acacGCTGCAAGGCAACGTAGCCCCAGCGTATGtta                            >  1:539173/1‑35 (MQ=255)     acacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTc                            >  1:3308928/1‑36 (MQ=255)     acacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTc                            <  1:2285619/36‑1 (MQ=255)      cacGCTGCAAGGCAACGTAGCTCCAGCGTATGTTcc                           <  1:2838817/36‑1 (MQ=25)      cacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTcc                           <  1:4933504/36‑1 (MQ=255)      cacGCTGCAAGGCAACGTAGCCCCAGCGTATGTTcc                           >  1:322309/1‑36 (MQ=255)       acGCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCt                          >  1:2833591/1‑36 (MQ=255)        cGCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTc                         <  1:3124998/36‑1 (MQ=255)         gCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCa                        >  1:5072320/1‑36 (MQ=255)         gCTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCa                        >  1:5111065/1‑36 (MQ=255)          cTGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCAt                       >  1:863483/1‑36 (MQ=255)           tGCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATc                      <  1:2353411/36‑1 (MQ=255)            gaaaGGCAACGTAGCCCCAGCGTATGTTCCTCATCg                     <  1:118627/34‑1 (MQ=25)            gCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATCg                     >  1:7225668/1‑36 (MQ=255)            gCAAGGCAACGTAGCCCCAGCGTATGTTCCTCATCg                     >  1:813299/1‑36 (MQ=255)             cAAGGCAGCGTAGCCCCAGCGTATGTTCCTCATCGt                    <  1:5380188/36‑1 (MQ=25)              aaGGCAACGTAGCCCCAGCGTATGTTCCTCCTCGTc                   >  1:3460714/1‑36 (MQ=25)              aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                   >  1:4070607/1‑36 (MQ=255)              aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                   <  1:837876/36‑1 (MQ=255)              aaGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTc                   >  1:6981796/1‑36 (MQ=255)               tggCAACGTAGCCCCAGCGTATGTTCCTCATCGTCa                  >  1:2198307/2‑36 (MQ=37)               atgCAACGTAGCCCCAGCGTATGTTCCTCATCGTCa                  <  1:4512492/34‑1 (MQ=25)               aGGCAACGTAGCCCCAGCGTATGTTCCTCATCGTCa                  <  1:4362032/36‑1 (MQ=255)                ggCTACGTAGCCCCAGCGTATGTTCCTCATCGTCaa                 <  1:4367480/36‑1 (MQ=25)                ggCAACGTAGCCCCAGCGTATGTTCCTCATCGTCaa                 <  1:2205778/36‑1 (MQ=255)                gcgaaCGTAGCCCCAGCGCATGTTCCTCATCGTCaa                 <  1:7328647/33‑1 (MQ=255)                 gaaaCGAAGCCCCAGCGAATGTTCCTCATCGTCaaa                <  1:4343362/34‑1 (MQ=21)                 gCAACGTAGCCCCAGCGTATGTTCCTCATCGTCaaa                >  1:1837897/1‑36 (MQ=255)                  cAACGTAGCCCCAGCGTATGTTCCTCATCGTCAAAt               <  1:3220794/36‑1 (MQ=255)                   aaCGTAGCCCCAGCGTATGTTCCTCATCGTCATATg              >  1:4429167/1‑36 (MQ=25)                   aaCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATg              >  1:6494826/1‑36 (MQ=255)                   aaCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATg              <  1:134193/36‑1 (MQ=255)                    aCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATGa             <  1:6028841/36‑1 (MQ=255)                    aCGTAGCCCCAGCGTATGTTCCTCATCGTCAAATGa             <  1:5819855/36‑1 (MQ=255)                      gTAGCCCCAGCGTATGTTCCTCATCGTCAAATGAAt           >  1:7376098/1‑36 (MQ=255)                      gTAGCCCCAGCGTATGTTCCTCATCGTCAAATGAAt           >  1:6454864/1‑36 (MQ=255)                       tAGTCCCAGCGTATGTTCCTCATCGTCAAATGAATa          <  1:5417027/36‑1 (MQ=25)                       tAGCCCCAGCGTATGTTCCTCATCGTCAAATGAATa          <  1:5971182/36‑1 (MQ=255)                       tAGCCCCAGCGTATGTTCCTCATCGTCAAATGAATa          <  1:4568354/36‑1 (MQ=255)                        aGCCCCAGCGTATGTTCCTCATCGTCAAATGAATaa         >  1:1322666/1‑36 (MQ=255)                        aGCCCCAGCGTATGTTCCTCATCGTAAAAGGACTaa         >  1:6275042/1‑36 (MQ=255)                         gCCCCAGCGTATGTTCCTCATCGTCAAATGAATAAc        >  1:4643924/1‑36 (MQ=255)                         gCCCCAGCGTATGTTCCTCATCGTCAAATGAATAAc        <  1:393560/36‑1 (MQ=255)                          ccctAGCGTATGTTCCTCATCGTCAAATGAATAAcc       <  1:5693696/36‑1 (MQ=25)                          ccccAGCGTATGTTCCTCATCGTCAAATGAATAAcc       >  1:3286873/1‑36 (MQ=255)                           cccAGCGTATGTTCCTCATCGTCAAATGAATAAccc      <  1:2957970/36‑1 (MQ=255)                            ccAGCGTATGTTCCTCATCGTCAAATGAATAAcccc     <  1:5707190/36‑1 (MQ=255)                             cAGCGTATGTTCCTCATCGTCAAATGAATAACCCCg    >  1:4985690/1‑36 (MQ=255)                             cAGCGTATGTTCCTCATCGTCAAATGAATAACCCCg    >  1:7354050/1‑36 (MQ=255)                              aGCGTATGTTCCTCATCGTCAAATGAATAACCCCGt   <  1:1333136/36‑1 (MQ=255)                               gcgtATGTTCCTCATCGTCAAATGAATAACCCCGtt  >  1:5250507/1‑36 (MQ=255)                                  |                                 AAATACACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTT  >  REL606/4201878‑4201943 Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG