Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,271,926	T→C	D38G (GAC→GGC)	yrbA ←	predicted DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,271,926	0	T	C	100.0%	339.4 / NA	88	D38G (GAC→GGC)	yrbA	predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (50/38);  total (50/38)

TCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCC  >  REL606/3271879‑3271972                                                |                                               tCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCaa                                             <  1:2505957/51‑1 (MQ=255) tCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCaa                                             <  1:2664893/51‑1 (MQ=255)  cAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCaaa                                            >  1:4426917/1‑51 (MQ=255)  cAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCaaa                                            <  1:3892053/51‑1 (MQ=255)  cAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCaaa                                            <  1:6523752/51‑1 (MQ=255)   aGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAaac                                           <  1:4145237/51‑1 (MQ=255)   aGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAaac                                           <  1:5661448/51‑1 (MQ=255)    gCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAaaca                                          <  1:1929741/51‑1 (MQ=255)     cGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAaacaa                                         >  1:1865255/1‑51 (MQ=255)      ggACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAaacaac                                        >  1:974488/1‑51 (MQ=255)         ccATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCa                                     >  1:4740880/1‑51 (MQ=255)         ccATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCa                                     <  1:393521/51‑1 (MQ=255)         ccATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCa                                     <  1:677691/51‑1 (MQ=255)          cATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCAc                                    >  1:4841305/1‑51 (MQ=255)          cATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCAc                                    >  1:962643/1‑51 (MQ=255)          cATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCAc                                    >  1:4264094/1‑51 (MQ=255)          cATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCAc                                    >  1:7114243/1‑51 (MQ=255)          cATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCAc                                    >  1:86051/1‑51 (MQ=255)           aTAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCAcc                                   >  1:6983778/1‑51 (MQ=255)              gACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCAc                                >  1:2476887/1‑51 (MQ=255)                ccGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACgg                              <  1:1092956/51‑1 (MQ=255)                ccGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACgg                              >  1:6292201/1‑51 (MQ=255)                ccGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACgg                              <  1:6019721/51‑1 (MQ=255)                ccGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACg                               <  1:53653/50‑1 (MQ=255)                 cGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGc                             >  1:5812016/1‑51 (MQ=255)                  gTCTGCTGTTTTTTAACCCGACTCATGTCGCCAAACAACTCACCCACGGCa                            >  1:5877958/1‑51 (MQ=255)                  gTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCa                            <  1:4305330/51‑1 (MQ=255)                  gTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCa                            >  1:2869875/1‑51 (MQ=255)                  gTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCa                            <  1:5735507/51‑1 (MQ=255)                  gTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCa                            <  1:4645377/51‑1 (MQ=255)                  gTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCa                            >  1:2183982/1‑51 (MQ=255)                  gTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCa                            <  1:4442465/51‑1 (MQ=255)                   tCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCaa                           >  1:2814865/1‑51 (MQ=255)                   tCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCaa                           >  1:7059253/1‑51 (MQ=255)                   tCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCaa                           >  1:1632700/1‑51 (MQ=255)                   tCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCaa                           >  1:1156638/1‑51 (MQ=255)                   tCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCaa                           <  1:7173515/51‑1 (MQ=255)                   tCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCaa                           <  1:3795284/51‑1 (MQ=255)                    ctgctgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAAt                          >  1:7214779/1‑51 (MQ=255)                    ctgctgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAAt                          >  1:5302507/1‑51 (MQ=255)                    ctgctgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAAt                          >  1:5182799/1‑51 (MQ=255)                     tgctgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATa                         <  1:2195672/51‑1 (MQ=255)                     tgctgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATa                         >  1:3546307/1‑51 (MQ=255)                     tgctgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATa                         <  1:1393931/51‑1 (MQ=255)                        tgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAAcc                      >  1:3788399/1‑51 (MQ=255)                        tgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAAcc                      >  1:1945196/1‑51 (MQ=255)                         gTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCt                     <  1:6514456/51‑1 (MQ=255)                           tttttAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGa                   <  1:4783372/51‑1 (MQ=255)                            ttttAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGaa                  >  1:657205/1‑51 (MQ=255)                             tttAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGaaa                 <  1:4442981/51‑1 (MQ=255)                             tttAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGaaa                 <  1:1559947/51‑1 (MQ=255)                             tttAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGaaa                 >  1:3902197/1‑51 (MQ=255)                               tAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGt               <  1:5560791/51‑1 (MQ=255)                                aaCCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTg              <  1:1222798/51‑1 (MQ=255)                                aaCCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTg              <  1:7063542/51‑1 (MQ=255)                                 aCCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTgg             <  1:4729192/51‑1 (MQ=255)                                 aCCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTgg             >  1:844452/1‑51 (MQ=255)                                   ccGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCt           >  1:3232181/1‑51 (MQ=255)                                   ccGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCt           >  1:2659515/1‑51 (MQ=255)                                    cGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTg          <  1:1785693/51‑1 (MQ=255)                                    cGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTg          >  1:3151248/1‑51 (MQ=255)                                    cGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTg          <  1:1384205/51‑1 (MQ=255)                                     gACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGc         <  1:424742/51‑1 (MQ=255)                                     gACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGc         >  1:1186556/1‑51 (MQ=255)                                     gACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGc         >  1:747184/1‑51 (MQ=255)                                     gACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGc         <  1:2119939/51‑1 (MQ=255)                                     gACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGc         >  1:6033472/1‑51 (MQ=255)                                      aCTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGcc        >  1:5231832/1‑51 (MQ=255)                                      aCTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGcc        <  1:6797151/51‑1 (MQ=255)                                      aCTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGcc        >  1:311257/1‑51 (MQ=255)                                       cTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCa       <  1:5247191/51‑1 (MQ=255)                                       cTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCa       >  1:2013154/1‑51 (MQ=255)                                       cTCATGCCGCCAAACAACTCACCCACGCCAATAACCTGAAAGTGGCTGCCa       >  1:2863050/1‑51 (MQ=255)                                        tCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCAt      >  1:5894751/1‑51 (MQ=255)                                        tCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCAt      >  1:1272213/1‑51 (MQ=255)                                        tCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCAt      >  1:3494499/1‑51 (MQ=255)                                        tCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCAt      >  1:3742838/1‑51 (MQ=255)                                        tCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCAt      <  1:6259092/51‑1 (MQ=255)                                         cATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATc     >  1:5478226/1‑51 (MQ=255)                                         cATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATc     <  1:6793350/51‑1 (MQ=255)                                         cATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATc     >  1:792572/1‑51 (MQ=255)                                         cATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATc     >  1:356548/1‑51 (MQ=255)                                          aTGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCg    >  1:7059529/1‑51 (MQ=255)                                          aTGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCg    >  1:1970726/1‑51 (MQ=255)                                          aTGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCg    <  1:2762036/51‑1 (MQ=255)                                           tGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGc   >  1:7152253/1‑51 (MQ=255)                                           tGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGc   <  1:3337399/51‑1 (MQ=255)                                           tGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGc   >  1:30790/1‑51 (MQ=255)                                            gccgccAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGcc  <  1:1288517/51‑1 (MQ=255)                                            gccgccAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGcc  >  1:1115225/1‑51 (MQ=255)                                            gccgccAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGcc  <  1:5597231/51‑1 (MQ=255)                                                |                                               TCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCC  >  REL606/3271879‑3271972

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑