Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	573,489	C→T	100%	G106G (GGC→GGT)	ybjD →	conserved hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	573,489	0	C	T	82.9%	40.7 / 2.1	29	G106G (GGC→GGT)	ybjD	conserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base C (5/0);  new base T (0/24);  total (5/24)
Fisher's exact test for biased strand distribution p-value = 8.42e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGCCAGGCCGACATCGGGTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAGCGTGA  >  minE/573422‑573544                                                                    |                                                        tGCCAGGCCGACATCGGGTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGc                                                         <  1:2646626/68‑1 (MQ=255)   ccAGGCCGACATCGGGTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTAt                                                      >  1:3174863/1‑69 (MQ=255)                            aTCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAg                             >  1:3430793/1‑69 (MQ=255)                            aTCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAg                             >  1:3276631/1‑69 (MQ=255)                            aTCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAg                             >  1:3112924/1‑69 (MQ=255)                            aTCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAg                             >  1:2482764/1‑69 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGTAgtgtca  <  1:2294363/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:2747692/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:894221/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:873789/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:3251478/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:3197767/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:3147365/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:3117207/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:3106076/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:2863185/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:2844542/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:270799/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:2445519/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:22686/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:2263740/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:2048155/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:1967843/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:1770850/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:1527134/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:1349312/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:1314988/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:1229022/69‑3 (MQ=255)                                                       cATGCACCGATGGTTATCACCGTACTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:545935/69‑3 (MQ=255)                                                       cATGCAACGATGGTTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAgtgtca  <  1:1097073/69‑3 (MQ=255)                                                                    |                                                        TGCCAGGCCGACATCGGGTTCGCCGTTATCGGCCGCTGGAAGCGTGCTGGACGCCATGCACCGATGGCTATCACCGTATTTTTTATCGTCTGGAAGGGGAGAGTGCGGAAGACGGCAGCGTGA  >  minE/573422‑573544

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