Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	287	C→T	intergenic (+32/‑50)	thrL → / → thrA	thr operon leader peptide/fused aspartokinase I and homoserine dehydrogenase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	287	0	C	T	100.0%	131.8 / NA	40	intergenic (+32/‑50)	thrL/thrA	thr operon leader peptide/fused aspartokinase I and homoserine dehydrogenase I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/13);  total (27/13)

TTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAA  >  W3110S.gb/230‑331                                                          |                                             ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAg                                       >  1:439690/1‑65 (MQ=255) ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAg                                       >  1:2546911/1‑65 (MQ=255) ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAg                                       >  1:1508514/1‑65 (MQ=255) ttACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAg                                       >  1:1430600/1‑65 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:1183347/1‑65 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:953626/1‑65 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:802810/1‑65 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:589232/1‑65 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:2678861/1‑65 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:2469718/1‑65 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:239794/1‑65 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:1295962/1‑65 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:1308522/1‑65 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCggg                                 >  1:1417743/1‑65 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCgg                                  >  1:430212/1‑64 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCg                                   >  1:1125132/1‑63 (MQ=255)       caGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCg                                   >  1:736400/1‑63 (MQ=255)                gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:2142439/1‑65 (MQ=255)                gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:2253127/1‑65 (MQ=255)                gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:832672/1‑65 (MQ=255)                gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:1410595/1‑65 (MQ=255)                gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:1710278/1‑65 (MQ=255)                gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:2717056/1‑65 (MQ=255)                gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:851805/1‑65 (MQ=255)                gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCtttttttt                        >  1:453220/1‑65 (MQ=255)                gTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCttcttttt                        >  1:2516791/1‑65 (MQ=255)                      gCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGAc                   >  1:491470/1‑64 (MQ=255)                                  aGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGa       <  1:1367739/64‑1 (MQ=255)                                  aGGAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGa       <  1:1944291/64‑1 (MQ=255)                                   ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCACAGGTAACGAgg     <  1:2428908/65‑1 (MQ=39)                                   ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAATGTAACGAgg     <  1:2902680/65‑1 (MQ=39)                                   ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg     <  1:2204621/65‑1 (MQ=255)                                   ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg     <  1:1286013/65‑1 (MQ=255)                                   ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg     <  1:1599546/65‑1 (MQ=255)                                   ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAgg     <  1:2804994/65‑1 (MQ=255)                                   ggAAACACAGAAAAAAGCCCGCATCTGACAGTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAgg     <  1:1326534/64‑1 (MQ=38)                                      aaCACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaa  <  1:1265998/65‑1 (MQ=255)                                      aaCACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaa  <  1:2302765/65‑1 (MQ=255)                                      aaCACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaa  <  1:1645211/65‑1 (MQ=255)                                      aaCACAGAAAAAAGCCCGCATCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaa  <  1:1510272/65‑1 (MQ=255)                                                          |                                             TTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACAGTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAA  >  W3110S.gb/230‑331

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