Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	928,085	G→A	S511S (AGC→AGT)	cydC ←	fused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	928,085	0	G	A	78.6%	117.2 / 28.3	56	S511S (AGC→AGT)	cydC	fused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (12/0);  new base A (0/44);  total (12/44)
Fisher's exact test for biased strand distribution p-value = 1.79e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

ATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCG  >  W3110S.gb/928056‑928140                              |                                                        atcatTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTc                      >  1:96048/1‑65 (MQ=255) atcatTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTc                      >  1:730939/1‑65 (MQ=255) atcatTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTc                      >  1:598649/1‑65 (MQ=255) atcatTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTc                      >  1:2417332/1‑65 (MQ=255) atcatTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTc                      >  1:2400189/1‑65 (MQ=255) atcatTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTc                      >  1:2173694/1‑65 (MQ=255) atcatTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTc                      >  1:1808711/1‑65 (MQ=255) atcatTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTc                      >  1:176836/1‑65 (MQ=255) atcatTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTc                      >  1:1178659/1‑65 (MQ=255) atcatTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTc                      >  1:1441997/1‑65 (MQ=255) atcatTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTc                      >  1:1221782/1‑65 (MQ=255) atcatTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTc                      >  1:1286895/1‑65 (MQ=255)                      aGGATCTGACTTTCGGTTGTGTCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:1752395/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:330845/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2701719/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2710101/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2767977/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2793932/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2799432/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2937107/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:88847/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:3028100/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:3032576/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2509206/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:349065/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:403028/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:440410/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:453102/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:510934/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:532278/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:586414/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:808568/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:1945640/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:1581645/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:1591943/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:1672259/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:176699/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:1774567/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:1405922/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:1821858/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:187026/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:1927053/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2662693/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2092394/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2152158/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2227160/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:1185581/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:244681/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:1085991/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2479649/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:956168/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2634427/64‑1 (MQ=255)                      aGGATCTGACTTTCGGTTGTGGCATCTAACCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2467802/64‑1 (MQ=255)                      aGGATCTGACATTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:1135274/64‑1 (MQ=255)                       ggATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTATGTTCATCCAGCAATACCAgcggcg  <  1:1549933/63‑1 (MQ=255)                       ggATCTGACTTTCGGTTGTGGCATCTAAGCCTTCTGTAGGTTCATCCAGCAATACCAgcggcg  <  1:2957300/63‑1 (MQ=255)                              |                                                        ATCATTTCTGCAAGCAATTCAAGGATCTGGCTTTCGGTTGTGGCATCTAAGCCTTCGGTAGGTTCATCCAGCAACACCAGTGGCG  >  W3110S.gb/928056‑928140

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