Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,160,543	C→T	100%	F7F (TTC→TTT)	mdtC →	multidrug efflux system, subunit C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,160,543	0	C	T	80.0%	61.3 / 12.0	30	F7F (TTC→TTT)	mdtC	multidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base C (5/1);  new base T (15/9);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTGGACCAAAAGCCGCTTTGCCCGTCATGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTCATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCTGTGCGGCATA  >  W3110S.gb/2160480‑2160600                                                                |                                                          gtgGACCAAAAGCCGCTTTGCCCGTCATGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTAcc                                                    >  1:693790/1‑71 (MQ=255) gtgGACCAAAAGCCGCTTTGCCCGTCATGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTAcc                                                    >  1:570572/1‑71 (MQ=255) gtgGACCAAAAGCCGCTTTGCCCGTCATGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTAcc                                                    >  1:1210442/1‑71 (MQ=255) gtgGACCAAAAGCCGCTTTGCCCGTCATGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTAcc                                                    >  1:1971088/1‑71 (MQ=255) gtgGACCAAAAGCCGCTTTGCCCGTCATGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTAcc                                                    >  1:1829058/1‑71 (MQ=255) gtgGACCAAAAGCCGCTTTGCCCGTCATGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTAcc                                                    >  1:1758623/1‑71 (MQ=255) gtgGACCAAAAGCCGCTTTGCCCGTCATGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTAcc                                                    >  1:1756023/1‑71 (MQ=255)                           aTGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCg                          >  1:717368/1‑71 (MQ=255)                           aTGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCg                          >  1:335773/1‑71 (MQ=255)                           aTGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCg                          >  1:1301543/1‑71 (MQ=255)                           aTGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCg                          >  1:1521040/1‑71 (MQ=255)                              aaGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCGGt                        <  1:962758/70‑1 (MQ=255)                              aaGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCGGt                        <  1:212286/70‑1 (MQ=255)                                     ggCGTAAGTGAAGTTTTTTGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGcc                    >  1:1122140/1‑67 (MQ=255)                                     ggCGTAAGTGAAGTTTTTTGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCa                   >  1:1666488/1‑68 (MQ=255)                                     ggCGTAAGTGAAGTTTTTTGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCa                   >  1:1845596/1‑68 (MQ=255)                                     ggCGTAAGTGAAGTTTTTTGCCCTCTTTATTTACCGCCCGGTGGCGACGATTATACTGTCGGTAGcc                    >  1:1082518/1‑67 (MQ=255)                                           aGTGAAGTTTTTTGCCCTCTTCATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCtgt          >  1:1714174/1‑71 (MQ=255)                                           aGTGAAGTTTTTTGCCCTCTTCATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCtgt          >  1:1711932/1‑71 (MQ=255)                                           aGTGAAGTTTTTTGCCCTCTTCATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCtgt          >  1:918351/1‑71 (MQ=255)                                           aGTGAAGTTTTTTGCCCTCTTCATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCtgt          >  1:643506/1‑71 (MQ=255)                                           aGTGAAGTTTTTTGCCCTCTTCATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCtgt          >  1:1045380/1‑71 (MQ=255)                                               aaGTTTTTTGCCCTCTTCATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCTGTGCgg      <  1:1615530/71‑1 (MQ=255)                                                  ttttttGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCTGTGCGGCAt   <  1:890026/71‑1 (MQ=255)                                                  ttttttGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCTGTGCGGCAt   <  1:1185295/71‑1 (MQ=255)                                                    ttttGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCTGTGCGGCATa  <  1:1729699/70‑1 (MQ=255)                                                    ttttGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCTGTGCGGCATa  <  1:1060963/70‑1 (MQ=255)                                                    ttttGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCTGTGCGGCATa  <  1:1866607/70‑1 (MQ=255)                                                    ttttGCCCTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCTGTGCGGCATa  <  1:1699381/70‑1 (MQ=255)                                                         cccTCTTTATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCTGTGCGGCATa  <  1:1039050/65‑1 (MQ=255)                                                                |                                                          GTGGACCAAAAGCCGCTTTGCCCGTCATGAAGAGGAGGCGTAAGTGAAGTTTTTTGCCCTCTTCATTTACCGCCCGGTGGCGACGATTTTACTGTCGGTTGCCATTACCCTGTGCGGCATA  >  W3110S.gb/2160480‑2160600

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