Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,403,465	G→A	100%	L436L (CTG→CTA)	rtn →	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,403,465	0	G	A	75.0%	16.4 / 3.7	12	L436L (CTG→CTA)	rtn	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (9/0);  total (9/3)
Fisher's exact test for biased strand distribution p-value = 4.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCC  >  minE/1403419‑1403517                                               |                                                     ggCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTgg                                        <  1:165839/61‑1 (MQ=255)                                      cGATTATCTAAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTccc  >  1:109747/1‑62 (MQ=255)                                      cGATTATCTAAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTccc  >  1:143594/1‑62 (MQ=255)                                      cGATTATCTAAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTccc  >  1:330580/1‑62 (MQ=255)                                      cGATTATCTAAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTccc  >  1:333594/1‑62 (MQ=255)                                      cGATTATCTAAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTccc  >  1:369670/1‑62 (MQ=255)                                      cGATTATCTAAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTccc  >  1:391592/1‑62 (MQ=255)                                      cGATTATCTAAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTccc  >  1:406566/1‑62 (MQ=255)                                      cGATTATCTAAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTccc  >  1:542670/1‑62 (MQ=255)                                      cGATTATCTAAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTccc  >  1:81518/1‑62 (MQ=255)                                           aTCTGAAAATTGACCGTGGATTTATCAACGCCATCg                       <  1:296600/36‑1 (MQ=255)                                           aTCTGAAAATTGACCGTGGATTTATCAACGCCATCg                       <  1:54257/36‑1 (MQ=255)                                               |                                                     GGCACAGCGCGCTTATCTATCTTGAGCGTTTTACGCTCGATTATCTGAAAATTGACCGTGGATTTATCAACGCCATCGGTACGGAAACGATCACTTCCC  >  minE/1403419‑1403517

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