Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,171,939	T→C	T5A (ACG→GCG)	ycfS ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,171,939	0	T	C	100.0%	118.9 / NA	33	T5A (ACG→GCG)	ycfS	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (7/26);  total (7/26)

CCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCT  >  W3110S.gb/1171885‑1171999                                                       |                                                             ccagcgccaCGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGCTTTGatcatc                                                    >  1:1040303/1‑65 (MQ=255) ccagcgccaCGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGCTTTGatcatc                                                    >  1:856485/1‑65 (MQ=255) ccagcgccaCGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGCTTTGatcatc                                                    >  1:553782/1‑65 (MQ=255) ccagcgccaCGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGCTTTGatcatc                                                    >  1:1470219/1‑65 (MQ=255) ccagcgccaCGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGCTTTGatcatc                                                    >  1:1700699/1‑65 (MQ=255) ccagcgccaCGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGCTTTGatcatc                                                    >  1:2239996/1‑65 (MQ=255) ccagcgccaCGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGCTTTGatcatc                                                    >  1:1880042/1‑65 (MQ=255)                                  ttAGCCAGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGt                   <  1:1093023/65‑1 (MQ=255)                                  ttAGCCAGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGt                   <  1:750784/65‑1 (MQ=255)                                  ttAGCCAGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGt                   <  1:2954268/65‑1 (MQ=255)                                  ttAGCCAGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGt                   <  1:2333214/65‑1 (MQ=255)                                  ttAGCCAGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGt                   <  1:222929/65‑1 (MQ=255)                                  ttAGCCAGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGt                   <  1:1492424/65‑1 (MQ=255)                                  ttAGCCAGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGt                   <  1:1238290/65‑1 (MQ=255)                                  ttAGCCAGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGt                   <  1:116387/65‑1 (MQ=255)                                        aGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATaaa             <  1:2159425/65‑1 (MQ=255)                                        aGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATaaa             <  1:2148106/65‑1 (MQ=255)                                        aGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATaaa             <  1:2354759/65‑1 (MQ=255)                                        aGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATaaa             <  1:2489816/65‑1 (MQ=255)                                        aGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATaaa             <  1:2570198/65‑1 (MQ=255)                                        aGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATaaa             <  1:2616087/65‑1 (MQ=255)                                        aGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATaaa             <  1:2727939/65‑1 (MQ=255)                                        aGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATaaa             <  1:876392/65‑1 (MQ=255)                                        aGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATaaa             <  1:2897825/65‑1 (MQ=255)                                        aGCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATaaa             <  1:1336208/65‑1 (MQ=255)                                         gCGAGAAAAACGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATaaa             <  1:1113306/64‑1 (MQ=255)                                              aaaaaCGCGCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGt                   <  1:929239/53‑1 (MQ=255)                                                   cgcgCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCt  <  1:935334/65‑1 (MQ=255)                                                   cgcgCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCt  <  1:673086/65‑1 (MQ=255)                                                   cgcgCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCt  <  1:3043874/65‑1 (MQ=255)                                                   cgcgCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCt  <  1:2767989/65‑1 (MQ=255)                                                   cgcgCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCt  <  1:1704257/65‑1 (MQ=255)                                                   cgcgCTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCt  <  1:1755130/65‑1 (MQ=255)                                                       |                                                             CCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCT  >  W3110S.gb/1171885‑1171999

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