Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,667,558	T→C	intergenic (‑6/‑130)	hcaR ← / → hcaE	DNA‑binding transcriptional regulator/3‑phenylpropionate dioxygenase, large (alpha) subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,667,558	0	T	C	100.0%	95.8 / NA	27	intergenic (‑6/‑130)	hcaR/hcaE	DNA‑binding transcriptional regulator/3‑phenylpropionate dioxygenase, large (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/13);  total (14/13)

CGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAG  >  W3110S.gb/2667490‑2667615                                                                     |                                                          cgcACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTccc                                                         <  1:4034421/71‑1 (MQ=255) cgcACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTccc                                                         <  1:3380901/71‑1 (MQ=255) cgcACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTccc                                                         <  1:2836466/71‑1 (MQ=255) cgcACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTccc                                                         <  1:2593171/71‑1 (MQ=255) cgcACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTccc                                                         <  1:1553400/71‑1 (MQ=255) cgcACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTccc                                                         <  1:1663560/71‑1 (MQ=255) cgcACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTccc                                                         <  1:2220781/71‑1 (MQ=255)        gTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTCCCCCTTGt                                                   >  1:1390701/1‑70 (MQ=255)        gTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTCCCCCTTGt                                                   >  1:2306051/1‑70 (MQ=255)        gTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTCCCCCTTGt                                                   >  1:2258962/1‑70 (MQ=255)                      cTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGtc                                    <  1:1187700/71‑1 (MQ=255)                      cTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGtc                                    <  1:83205/71‑1 (MQ=255)                      cTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGtc                                    <  1:2397803/71‑1 (MQ=255)                                     gaAATAGCGTAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGt                                     >  1:745835/1‑55 (MQ=255)                                     gaAATAGCGTAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGt                                     >  1:242455/1‑55 (MQ=255)                                     gaAATAGCGTAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGt                                     >  1:2602866/1‑55 (MQ=255)                                     gaAATAGCGTAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGt                                     >  1:3031930/1‑55 (MQ=255)                                         tAGCGTAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAAc                 <  1:3776916/71‑1 (MQ=255)                                              tAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGt                                     >  1:2978320/1‑46 (MQ=255)                                              tAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTCACTAttt            <  1:3256982/71‑1 (MQ=255)                                              tAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAttt            <  1:1667356/71‑1 (MQ=255)                                              tAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAttt            <  1:1292519/71‑1 (MQ=255)                                              tAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAttt            <  1:1306362/71‑1 (MQ=255)                                              tAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAttt            >  1:2878550/1‑71 (MQ=255)                                              tAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAttt            >  1:1373209/1‑71 (MQ=255)                                              tAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAttt            <  1:2149447/71‑1 (MQ=255)                                              tAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAtt             >  1:1587179/1‑70 (MQ=255)                                              tAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAtt             >  1:4079103/1‑70 (MQ=255)                                              tAAATGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTAtt             >  1:900428/1‑70 (MQ=255)                                                  tGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTcaca        <  1:1402080/71‑1 (MQ=255)                                                  tGCCGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTcaca        <  1:1184625/71‑1 (MQ=255)                                                     cGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACatat     <  1:664667/71‑1 (MQ=255)                                                     cGTAGTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACatat     <  1:826350/71‑1 (MQ=255)                                                         gTTCCATCACCTCCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAg  >  1:3380209/1‑70 (MQ=255)                                                                     |                                                          CGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAG  >  W3110S.gb/2667490‑2667615

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑