Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,429,248	C→T	E76K (GAA→AAA)	yjfY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,429,248	0	C	T	100.0%	170.0 / NA	51	E76K (GAA→AAA)	yjfY	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/15);  total (36/15)

GATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAA  >  W3110S.gb/4429203‑4429305                                              |                                                          gATAAATATCCGCGCTGGCGTTCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:1056112/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:1088176/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:2268347/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:2205085/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:2176369/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:2088939/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:2037997/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:2444690/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:2536172/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:1930045/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:2868062/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:294740/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:477715/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:821879/61‑1 (MQ=255) gATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTAc                                            <  1:1314980/61‑1 (MQ=255)  aTAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTa                                             >  1:2727215/1‑59 (MQ=255)  aTAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTa                                             >  1:2618737/1‑59 (MQ=255)  aTAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTa                                             >  1:492963/1‑59 (MQ=255)  aTAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTa                                             >  1:94427/1‑59 (MQ=255)  aTAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTa                                             >  1:23854/1‑59 (MQ=255)  aTAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTa                                             >  1:2305947/1‑59 (MQ=255)  aTAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTa                                             >  1:986100/1‑59 (MQ=255)  aTAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTa                                             >  1:1398728/1‑59 (MQ=255)  aTAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTTCCGCATCAGAATTa                                             >  1:1812036/1‑59 (MQ=255)                              ccGTTACTCCCCGGTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTg               >  1:657528/1‑61 (MQ=255)                              ccGTTACTCCCCGGTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTg               >  1:150554/1‑61 (MQ=255)                              ccGTTACTCCCCGGTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGt              >  1:1209312/1‑62 (MQ=255)                              ccGTTACTCCCCGGTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGt              >  1:1354455/1‑62 (MQ=255)                              ccGTTACTCCCCGGTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGt              >  1:2454172/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:2813748/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:842056/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:628494/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:527021/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:508619/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:1477776/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:330916/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:1734818/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:1790536/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:2245458/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:2662894/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:1918489/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:2615584/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:2609188/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:1973983/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:202971/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:2432043/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:2413029/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:2393860/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTaa  >  1:2189824/1‑62 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTa   >  1:2784252/1‑61 (MQ=255)                                          ggTTTCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTAATTTaa  >  1:348297/1‑62 (MQ=255)                                              |                                                          GATAAATATCCGCGCTGGCGTGCATATTGCCGTTACTCCCCGGTTCCCGCATCAGAATTACGTGGTAGTACGTTGCGCCTTGCGCATCTGTTTCTTCATTTAA  >  W3110S.gb/4429203‑4429305

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