Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	364,867	T→C	T221T (ACA→ACG)	lacZ ←	beta‑D‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	364,867	0	T	C	94.2%	128.6 / ‑2.1	52	T221T (ACA→ACG)	lacZ	beta‑D‑galactosidase
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base C (0/49);  total (3/49)
Fisher's exact test for biased strand distribution p-value = 4.52e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCATTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAA  >  W3110S.gb/364829‑364903                                       |                                     tcatTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACga               >  1:256238/1‑62 (MQ=255) tcatTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACga               >  1:188175/1‑62 (MQ=255) tcatTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACga               >  1:537473/1‑62 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:491931/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:300833/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:318965/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:324635/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:348754/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:376047/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:39212/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:443576/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:451573/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:473122/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:473611/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:284379/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:516664/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:624932/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:670152/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:679758/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:695825/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:747488/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:747541/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:76553/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:775059/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:978902/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:999434/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1230391/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1019069/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1021230/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1041387/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1104466/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1105224/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:11119/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1177459/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:122481/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1229148/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1230000/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1012529/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1246602/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1248126/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1271942/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1414937/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:1442171/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:148790/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:160064/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:189323/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:206204/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:230204/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:247934/51‑1 (MQ=255)                         aaaTCGCTGATTTGCGTGGTCGGGTTATGCAGCAACGAGACGTCACGGaaa  <  1:593846/51‑1 (MQ=255)                          aaTCGCTGATTTGCGTGGTCGGTTTATGCAGCAACGAGACGTCACGGaaa  <  1:74772/50‑1 (MQ=255)                                   tttGCGTGGTCGGTTTATGCAGCAACGAGACGTCACTGaaa  <  1:1272766/41‑1 (MQ=37)                                       |                                     TCATTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAA  >  W3110S.gb/364829‑364903

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