Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,841,531	G→A	R57* (CGA→TGA)	hycI ←	protease involved in processing C‑terminal end of HycE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,841,531	0	G	A	100.0%	70.1 / NA	22	R57* (CGA→TGA)	hycI	protease involved in processing C‑terminal end of HycE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/15);  total (7/15)

GATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACC  >  W3110S.gb/2841472‑2841582                                                            |                                                    gatgatGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCatg                                                   <  1:524683/62‑1 (MQ=255) gatgatGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCatg                                                   <  1:4263602/62‑1 (MQ=255) gatgatGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCatg                                                   <  1:3872077/62‑1 (MQ=255) gatgatGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCatg                                                   <  1:3859029/62‑1 (MQ=255) gatgatGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCatg                                                   <  1:3834498/62‑1 (MQ=255) gatgatGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCatg                                                   <  1:3248900/62‑1 (MQ=255) gatgatGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCatg                                                   <  1:2637117/62‑1 (MQ=255) gatgatGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCatg                                                   <  1:2631873/62‑1 (MQ=255) gatgatGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCatg                                                   <  1:2089159/62‑1 (MQ=255)                      gTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCATGGCGGTCGCAGTTCACGGATAg                              <  1:2501710/62‑1 (MQ=255)                             cccATATCCGTGGCGTCGACAATCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGAt                       >  1:4379841/1‑62 (MQ=255)                                     cGTGGCGTCGACAATCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGtttt               <  1:3533321/62‑1 (MQ=255)                                     cGTGGCGTCGACAATCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGtttt               <  1:3538754/62‑1 (MQ=255)                                     cGTGGCGTCGACAATCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGtttt               <  1:3828254/62‑1 (MQ=255)                                     cGTGGCGTCGACAATCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGtttt               <  1:1745819/62‑1 (MQ=255)                                          cGTCGACAATCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGAt                       >  1:2695884/1‑49 (MQ=255)                                          cGTCGACAATCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGt          >  1:4481191/1‑62 (MQ=255)                                          cGTCGACAATCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGt          >  1:2729022/1‑62 (MQ=255)                                          cGTCGACAATCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGt          >  1:725905/1‑62 (MQ=255)                                          cGTCGACAATCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGt          >  1:4451264/1‑62 (MQ=255)                                          cGTCGACAATCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGt          >  1:67849/1‑62 (MQ=255)                                                  aTCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAcc  <  1:706990/62‑1 (MQ=255)                                                  aTCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAcc  <  1:3785808/62‑1 (MQ=255)                                                  aTCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAcc  <  1:749908/62‑1 (MQ=255)                                                  aTCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAcc  <  1:972825/62‑1 (MQ=255)                                                  aTCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAcc  <  1:1180911/62‑1 (MQ=255)                                                  aTCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAcc  <  1:3259106/62‑1 (MQ=255)                                                  aTCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAcc  <  1:2909888/62‑1 (MQ=255)                                                  aTCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAcc  <  1:269106/62‑1 (MQ=255)                                                  aTCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAcc  <  1:1442896/62‑1 (MQ=255)                                                   tCAGCAGTCATGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTAcc  <  1:540279/61‑1 (MQ=255)                                                            |                                                    GATGATGCGGATCTCGCCGGGGTTTAGCCCCATATCCGTGGCGTCGACAATCAGCAGTCGTGTCGGGCGCAGTTCACGGATAGCGACGATGTCGTTTTCCGGTGCGCTACC  >  W3110S.gb/2841472‑2841582

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