Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,867,299	G→A	71.1%	C18C (TGC→TGT)	ygcS ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,867,299	0	G	A	71.1%	32.9 / 13.6	38	C18C (TGC→TGT)	ygcS	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (9/2);  new base A (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 = 8.24e-01

ACCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTAC  >  minE/1867236‑1867358                                                                |                                                            aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:2094494/1‑68 (MQ=255) aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:268067/1‑68 (MQ=255) aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:2694982/1‑68 (MQ=255) aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:3313739/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) aCCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGt                                                         >  1:989710/1‑68 (MQ=255)           ccGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTa                                              <  1:2294204/69‑1 (MQ=255)                  cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:578419/69‑1 (MQ=255)                  cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTa                                       <  1:1269732/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: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:2673060/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:3014742/1‑69 (MQ=255)                      aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATTCGACAGTGAAAACGGTTAAGCGGTAAatc                                   >  1:2954121/1‑69 (MQ=255)                                                       agcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACAGGTGAAGTGTTCATCTTTAc  <  1:980716/69‑1 (MQ=255)                                                                |                                                            ACCAATGACGCCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTAC  >  minE/1867236‑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: ‑