Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,131,742	(G)5→6	coding (1138/1227 nt)	yhfW →	predicted mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,131,737	1	.	G	100.0%	95.9 / NA	28	V378G (GTG→GGG)	yhfW	predicted mutase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (21/7);  total (21/7)

TGCCAGTGCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGT‑GGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATGGTCGCTCTTTTCT  >  W3110S.gb/4131668‑4131797                                                                       |                                                             tgCCAGTGCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGAtgt‑g                                                             <  1:3742151/71‑1 (MQ=255) tgCCAGTGCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGAtgt‑g                                                             <  1:193483/71‑1 (MQ=255) tgCCAGTGCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGAtgt‑g                                                             <  1:498551/71‑1 (MQ=255) tgCCAGTGCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGAtgt‑g                                                             <  1:300839/71‑1 (MQ=255)     aGTGCTGGTTTATCAGCAAGTGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGT‑ggggg                                                         <  1:62387/71‑1 (MQ=255)     aGTGCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGT‑ggggg                                                         <  1:3476431/71‑1 (MQ=255)     aGTGCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGT‑ggggg                                                         <  1:2626635/71‑1 (MQ=255)     aGTGCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGGGTGCGCACCACGCTTTCTGATGT‑ggggg                                                         <  1:2513880/71‑1 (MQ=255)        gCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGT‑ggggggc                                                       >  1:2863205/1‑68 (MQ=255)        gCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGT‑ggggggc                                                       >  1:2837950/1‑68 (MQ=255)        gCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGT‑ggggggc                                                       >  1:1739653/1‑68 (MQ=255)        gCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGT‑ggggggc                                                       >  1:3170277/1‑68 (MQ=255)        gCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGT‑ggggggc                                                       >  1:2771431/1‑68 (MQ=255)        gCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGT‑ggggggc                                                       >  1:2383324/1‑68 (MQ=255)        gCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGT‑ggggggc                                                       >  1:3053587/1‑68 (MQ=255)                     cAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAAt                                          >  1:883324/1‑71 (MQ=255)                     cAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAAt                                          >  1:679045/1‑71 (MQ=255)                     cAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAAt                                          >  1:2131903/1‑71 (MQ=255)                     cAAGGGATGATCGCAACGCAGCTCGGTGTGCGCACCACGCTTACTGATGTGGGGGGCTACCg                                                   >  1:3678343/1‑62 (MQ=255)                     cAAGGGATGATAGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAAt                                          >  1:411999/1‑71 (MQ=255)                                 gCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCgcgcgc                              >  1:2773392/1‑71 (MQ=255)                                 gCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCgcgcgc                              >  1:3934316/1‑71 (MQ=255)                                 gCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCgcgcgc                              >  1:552667/1‑71 (MQ=255)                                 gCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCgcgcgc                              >  1:2600239/1‑71 (MQ=255)                                 gCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCgcgcgc                              >  1:888055/1‑71 (MQ=255)                                 gCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCgcgcgc                              >  1:1585802/1‑71 (MQ=255)                                   tACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCgcgcgc                              >  1:20822/1‑69 (MQ=255)                                   tACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCgcgcgc                              >  1:4099016/1‑69 (MQ=255)                                                 tgCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATg               <  1:3334370/70‑1 (MQ=255)                                                 tgCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATg               <  1:2031910/70‑1 (MQ=255)                                                 tgCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATg               <  1:3670737/70‑1 (MQ=255)                                                 tgCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATg               <  1:1616601/70‑1 (MQ=255)                                                 tgCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATg               <  1:1602692/70‑1 (MQ=255)                                                 tgCGCACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATg               <  1:4009371/70‑1 (MQ=255)                                                   cgcACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATGGTc            >  1:91462/1‑71 (MQ=255)                                                   cgcACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATGGTc            >  1:1396720/1‑71 (MQ=255)                                                   cgcACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATGGTc            >  1:3452084/1‑71 (MQ=255)                                                   cgcACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATGGTc            >  1:294951/1‑71 (MQ=255)                                                   cgcACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATGGTc            >  1:2867197/1‑71 (MQ=255)                                                   cgcACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATGGTc            >  1:2360351/1‑71 (MQ=255)                                                   cgcACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATGGTc            >  1:2166031/1‑71 (MQ=255)                                                   cgcACCACGCTTTCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATGGTc            >  1:1720028/1‑71 (MQ=255)                                                               tCTGATGTGGGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATGGTCGCtcctttct  <  1:1770003/69‑1 (MQ=255)                                                                       |                                                             TGCCAGTGCTGGTTTATCAGCAAGGGATGATCGCTACGCAGCTCGGTGTGCGCACCACGCTTTCTGATGT‑GGGGGCTACCGTGTGTGAATTTTTCCGCGCGCCACCGCCACAAAATGGTCGCTCTTTTCT  >  W3110S.gb/4131668‑4131797

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: ‑