Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	599,692	G→A	20.0%	A586T (GCG→ACG)	cusA →	copper/silver efflux system, membrane component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	599,692	0	G	A	20.0%	51.3 / 8.0	25	A586T (GCG→ACG)	cusA	copper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base G (11/9);  new base A (3/2);  total (14/11)
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.23e-01

ATCAATGAAGGCGACTTGTTGTATATGCCATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAA  >  W3110S.gb/599629‑599760                                                                |                                                                     aTCAATGAAGGCGACTTGTTGTATATGCCATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAg                                                                  <  1:2084736/68‑1 (MQ=255)   cAATGAAGGCGACTTGTTGTATATGCCATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTAtgc                                                             >  1:2077495/1‑71 (MQ=255)   cAATGAAGGCGACTTGTTGTATATGCCATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTAtgc                                                             >  1:3301779/1‑71 (MQ=255)        aaGGCGACTTGTTGTATATGCCATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGACGAGTATGCTGCa                                                         >  1:969365/1‑70 (MQ=255)              aCTTGTTGTATATGCCATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGACGAGTATGCTGCAAAAAAcc                                                  <  1:3862333/71‑1 (MQ=255)                ttgttgTATATGCCATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGa                                                >  1:1336326/1‑71 (MQ=255)                ttgttgTATATGCCATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGa                                                >  1:1398174/1‑71 (MQ=255)                  gttgtATATGCCATCGACGCTGCCGGTGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGAc                                               >  1:3795894/1‑70 (MQ=255)                          tGCCATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGACAAGCTAAt                                       <  1:1162424/70‑1 (MQ=255)                           gCCATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGACAAGCTAATTa                                     >  1:2127315/1‑71 (MQ=255)                           gCCATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGACAAGCTAATTa                                     >  1:2480212/1‑71 (MQ=255)                            ccATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGACAAGCTAATTa                                     <  1:1882551/70‑1 (MQ=255)                            ccATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGACAAGCTAATTa                                     <  1:2120945/70‑1 (MQ=255)                            ccATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGACAAGCTAATTa                                     <  1:2164137/70‑1 (MQ=255)                               tCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCaaaa                                                      <  1:975694/50‑1 (MQ=255)                               tCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCaaaa                                                      <  1:2819750/50‑1 (MQ=255)                                cGACGCTGCCGGGGATTTCCGCAGCAGAGGCGACGAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGc                                >  1:3129063/1‑71 (MQ=255)                                cGACGCTGCCGGGGATTTCCGCAGCAGAGGCGACGAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGc                                >  1:498798/1‑71 (MQ=255)                                 gACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCaaaaaa                                                    >  1:928851/1‑50 (MQ=255)                                 gACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGCg                               >  1:3734474/1‑71 (MQ=255)                                 gACGCTGCCGGGGATTTCCGCAGCAGAGACGGCGAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGCg                               >  1:1074085/1‑71 (MQ=255)                                        ccGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGCGTACCTGa                        <  1:580375/71‑1 (MQ=255)                                           gggATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGCGTACCTGa                        <  1:3010346/68‑1 (MQ=255)                                            ggATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGCGTACCTGAAGTg                    >  1:1631845/1‑71 (MQ=255)                                                   cGCAGCAGAGGCGACGAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGCGTACCTGAAGTGGCGCggg             <  1:926374/71‑1 (MQ=255)                                                             gcggcgAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGCGTACCTGAAGt                     >  1:169332/1‑53 (MQ=255)                                                             gcggcgAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCaa   <  1:722406/71‑1 (MQ=255)                                                              cggcgAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCaaa  <  1:949930/71‑1 (MQ=255)                                                                |                                                                     ATCAATGAAGGCGACTTGTTGTATATGCCATCGACGCTGCCGGGGATTTCCGCAGCAGAGGCGGCGAGTATGCTGCAAAAAACCGACAAGCTAATTATGAGCGTACCTGAAGTGGCGCGGGTATTTGGCAAA  >  W3110S.gb/599629‑599760

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