Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	3,035,425	T→C	Q446Q (CAA→CAG)	recJ ←	ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,035,425	0	T	C	95.6%	73.2 / ‑5.7	45	Q446Q (CAA→CAG)	recJ	ssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (0/43);  total (0/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGC  >  W3110S.gb/3035404‑3035469                      |                                             aaCCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACccg      <  1:45682/62‑1 (MQ=255)  aCCAGTTCGCCAAACCGTTGTTGAAAGAGGTTGAATTTATCCTCTTCCAGCGACAAACccg      <  1:45488/61‑1 (MQ=255)     aGTTCGCCAAAGCTTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:71122/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:64338/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:47136/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:49461/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:52705/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:53906/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:54901/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:55338/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:59781/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:60056/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:61847/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:46102/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:66779/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:67629/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:68773/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:75144/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:79044/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:79103/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:79776/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:81561/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:81876/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:30144/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:11911/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:13097/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:18141/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:18331/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:19119/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:19779/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:2086/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:22305/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:27628/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:46519/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:35288/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:37777/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:39286/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:39482/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:41195/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:41347/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:45858/62‑1 (MQ=255)     aGTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:10180/62‑1 (MQ=255)      gTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:35421/61‑1 (MQ=255)      gTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:79214/61‑1 (MQ=255)      gTTCGCCAAAGCGTTGCTGAAAGAGTTCGAATTTATCCTCTTCCAGCGACAATCCtgccgc  <  1:4149/61‑1 (MQ=255)                      |                                             AACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACAAACCCGCCGC  >  W3110S.gb/3035404‑3035469

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