Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,867,296	G→T	71.0%	R19R (CGC→CGA)	ygcS ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,867,296	0	G	T	71.0%	36.4 / 22.0	38	R19R (CGC→CGA)	ygcS	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (9/2);  new base T (15/12);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

ATGGCGTAACCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTAC  >  minE/1867228‑1867358                                                                     |                                                               atGGCGTAACCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATg                                                                >  1:3135069/1‑69 (MQ=255)         aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:989710/1‑68 (MQ=255)         aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:2694982/1‑68 (MQ=255)         aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:268067/1‑68 (MQ=255)         aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:3313739/1‑68 (MQ=255)         aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:2094494/1‑68 (MQ=255)         aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:1829361/1‑68 (MQ=255)         aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:864337/1‑68 (MQ=255)         aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:1514271/1‑68 (MQ=255)         aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:898133/1‑68 (MQ=255)                   ccGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTa                                              <  1:2294204/69‑1 (MQ=255)                          cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:578419/69‑1 (MQ=255)                          cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:552343/69‑1 (MQ=255)                          cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:525929/69‑1 (MQ=255)                          cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:3281515/69‑1 (MQ=255)                          cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:3190791/69‑1 (MQ=255)                          cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:914773/69‑1 (MQ=255)                          cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:1269732/69‑1 (MQ=255)                          cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:1198411/69‑1 (MQ=255)                          cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:1281800/69‑1 (MQ=255)                          cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:2618776/69‑1 (MQ=255)                          cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:2612878/69‑1 (MQ=255)                          cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:2458770/69‑1 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:465902/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:962168/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:1341497/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:1632231/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:2042484/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:557852/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:2110363/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:2954121/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:450208/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:3399695/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:3357742/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:324395/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:2640211/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:2673060/1‑69 (MQ=255)                              aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:3014742/1‑69 (MQ=255)                                                               agcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACAGGTGAAGTGTTCATCTTTAc  <  1:980716/69‑1 (MQ=255)                                                                     |                                                               ATGGCGTAACCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTAC  >  minE/1867228‑1867358

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