Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,382,164	G→T	A173D (GCT→GAT)	folE ←	GTP cyclohydrolase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,382,164	0	G	T	76.3%	50.0 / 27.2	38	A173D (GCT→GAT)	folE	GTP cyclohydrolase I
Reads supporting (aligned to +/- strand):  ref base G (9/0);  new base T (0/29);  total (9/29)
Fisher's exact test for biased strand distribution p-value = 6.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

TGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTG  >  minE/1382143‑1382213                      |                                                  tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:1312267/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:915058/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:765863/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:68970/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:644643/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:634967/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:49879/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:3140542/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:3071695/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:2917779/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:2738266/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:2726886/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:2468000/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:234672/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:2283020/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:2251426/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:2143269/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:2141155/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:1650531/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:1614110/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:159732/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:1298248/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:1582101/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:2221780/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:1496771/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:1561489/71‑1 (MQ=255) tGCACCGCGTCGATCGAGACATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:1512051/71‑1 (MQ=255)          tcgatcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgc    >  1:3200020/1‑60 (MQ=255)          tcgatcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgc    >  1:982125/1‑60 (MQ=255)          tcgatcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgc    >  1:953208/1‑60 (MQ=255)          tcgatcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgc    >  1:1485306/1‑60 (MQ=255)          tcgatcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgc    >  1:594638/1‑60 (MQ=255)          tcgatcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgc    >  1:338654/1‑60 (MQ=255)          tcgatcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgc    >  1:3017330/1‑60 (MQ=255)          tcgatcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgc    >  1:1764357/1‑60 (MQ=255)          tcgatcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgc    >  1:2152189/1‑60 (MQ=255)            gatcgaGACATCCACGTTATTGGTGGCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:1338257/60‑1 (MQ=255)                gagaCATCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTtgctg  <  1:1582871/56‑1 (MQ=255)                      |                                                  TGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTG  >  minE/1382143‑1382213

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