Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,353,381:1	+T	100%	coding (1034/1431 nt)	yihF ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,353,381	1	.	T	80.0%	19.3 / 2.3	10	coding (1034/1431 nt)	yihF	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (8/0);  total (10/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CGTTCTTCCATCTTACCGGTTGTTTAATGGTCGGCTCAC‑TTTTTTGTAATAACGGTAAATACTCTTTGAACAATGCAGCATTCACTTCTTGCAGAGC  >  minE/2353343‑2353439                                        |                                                           cGTTCTTCCATCTTACCGGTTGTTTAATGGTCGGCTCAC‑TTTTTTGTAATAACGGTAAATACTCTTTGaa                             >  1:1588639/1‑70 (MQ=255) cGTTCTTCCATCTTACCGGTTGTTTAATGGTCGGCTCAC‑TTTTTTGTAATAACGGTAAATACTCTTTGaa                             >  1:794763/1‑70 (MQ=255)                             ggTCGGCTCACTTTTTTTGTAATAACGGTAAATACTCTTTGAACAATGCAGCATTCACTTCTTGCagag   >  1:2920619/1‑69 (MQ=255)                             ggTCGGCTCACTTTTTTTGTAATAACGGTAAATACTCTTTGAACAATGCAGCATTCACTTCTTGCATAGc  >  1:743617/1‑70 (MQ=255)                             ggTCGGCTCACTTTTTTTGTAATAACGGTAAATACTCTTTGAACAATGCAGCATTCACTTCTTGCAGAGc  >  1:2129460/1‑70 (MQ=255)                             ggTCGGCTCACTTTTTTTGTAATAACGGTAAATACTCTTTGAACAATGCAGCATTCACTTCTTGCAGAGc  >  1:226648/1‑70 (MQ=255)                             ggTCGGCTCACTTTTTTTGTAATAACGGTAAATACTCTTTGAACAATGCAGCATTCACTTCTTGCAGAGc  >  1:272876/1‑70 (MQ=255)                             ggTCGGCTCACTTTTTTTGTAATAACGGTAAATACTCTTTGAACAATGCAGCATTCACTTCTTGCAGAGc  >  1:546535/1‑70 (MQ=255)                             ggTCGGCTCACTTTTTTTGTAATAACGGTAAATACTCTTTGAACAATGCAGCATTCACTTCTTGCAGAGc  >  1:669025/1‑70 (MQ=255)                             ggTCGGCTCACTTTTTTTGTAATAACGGTAAATACTCATTGAACAATGCAGCATTCACTTCTTGCAGAGc  >  1:1077887/1‑70 (MQ=255)                                        |                                                           CGTTCTTCCATCTTACCGGTTGTTTAATGGTCGGCTCAC‑TTTTTTGTAATAACGGTAAATACTCTTTGAACAATGCAGCATTCACTTCTTGCAGAGC  >  minE/2353343‑2353439

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