Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,123,875	G→A	43.7%	R255H (CGC→CAC)	damX →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,123,875	0	G	A	43.7%	2.2 / 15.9	16	R255H (CGC→CAC)	damX	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (3/6);  new base A (3/4);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGCCG  >  W3110S.gb/4123833‑4123939                                           |                                                                 acacGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGa                                       <  1:1483747/70‑1 (MQ=255) acacGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGa                                       <  1:17411/70‑1 (MQ=255)   acGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAAcc                                    >  1:1665168/1‑71 (MQ=255)      gCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTcag                                                     <  1:694835/51‑1 (MQ=255)           aaCGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAAc                            >  1:936797/1‑71 (MQ=255)           aaCGCAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGaaaaaa                             >  1:1235961/1‑70 (MQ=255)               cAAACCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgca                        <  1:1693134/71‑1 (MQ=255)                 aaCCGCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaag                      >  1:1566240/1‑71 (MQ=255)                     gCTGAACGTCCGTCCACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCgcaagcaa                   <  1:142225/70‑1 (MQ=255)                           cgtccgtcCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGaa            <  1:211679/71‑1 (MQ=255)                              ccgtccACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGaaaa          <  1:1387977/70‑1 (MQ=255)                                gtccACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACgg       <  1:826306/71‑1 (MQ=255)                                 tccACTACGCACCCAGCTCGTCAGCAGGCGGTGATTGAACCGaa                                 >  1:1016823/1‑44 (MQ=255)                                     cTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGaa            <  1:1319531/61‑1 (MQ=255)                                     cTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGcc   >  1:992663/1‑70 (MQ=255)                                     cTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGCCg  <  1:1316024/71‑1 (MQ=255)                                           |                                                                 ACACGGCGAAAACGCAAACCGCTGAACGTCCGTCCACTACGCGCCCAGCTCGTCAGCAGGCGGTGATTGAACCGAAAAAACCGCAAGCAACCGTGAAAACGGAGCCG  >  W3110S.gb/4123833‑4123939

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