Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	543,970	C→T	50.0%	intergenic (+113/‑130)	yliL → / → mntR	hypothetical protein/transcriptional regulator of mntH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	543,970	0	C	T	50.0%	‑5.7 / 20.1	16	intergenic (+113/‑130)	yliL/mntR	hypothetical protein/transcriptional regulator of mntH
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base T (2/6);  total (4/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGACCTCCTGACTTGCTAATCCCGTCGATCCTTGAGGGATGATTGCATTACATACAGATATAGCACAGGCTATATTATATAGCTA  >  minE/543915‑544000                                                        |                               aaGACCTCCTGACTTGCTAATCCCGTCGATCCTTGAGGGATGATTGCATTACATACAGatat                          <  1:189840/62‑1 (MQ=255) aaGACCTCCTGACTTGCTAATCCCGTCGATCCTTGAGGGATGATTGCATTACATACAGatat                          <  1:465449/62‑1 (MQ=255)      ctcctGACTTGCTAATCCCGTCGATCCTTGAGGGATGATTGCATTACATATAGATATAGcac                     <  1:149289/62‑1 (MQ=255)      ctcctGACTTGCTAATCCCGTCGATCCTTGAGGGATGATTGCATTACATATAGATATAGcac                     <  1:174838/62‑1 (MQ=255)      ctcctGACTTGCTAATCCCGTCGATCCTTGAGGGATGATTGCATTACATATAGATATAGcac                     <  1:413164/62‑1 (MQ=255)       tcctGACTTGCTAATCCCGTCGATCCTTGAGGGATGATTGCATTACATATAGATATAGcac                     <  1:402870/61‑1 (MQ=255)       tcctGACTTGCTAATCCCGTCGATCCTTGAGGGATGATTGCATTACATACAGATATAGcaca                    <  1:329975/62‑1 (MQ=255)          tGACTTGCTAATCCCGTCGATCCTTGAGGGATGATTGCATTACATATAGATATAGCACAGGc                 >  1:175614/1‑62 (MQ=255)          tGACTTGCTAATCCCGTCGATCCTTGAGGGATGATTGCATTACATATAGATATAGCACAGGc                 >  1:176386/1‑62 (MQ=255)           gACTTGCTAATCCCGTCGATCCTTGAGGGATGATTGCATTACATATAGATATAGCACAGGc                 <  1:208696/61‑1 (MQ=255)                     tCCCGTCGATCCTTGAGGGATGATTGCATTACATACAGATATAGCACAGGCtatatta          >  1:243143/1‑58 (MQ=255)                     tCCCGTCGATCCTTGAGGGATGATTGCATTACATACAGATATAGCACAGGCtatatta          >  1:345317/1‑58 (MQ=255)                      cccGTCGATCCTTGAGGGATGATTGCATTACATATAGATATAGcac                     <  1:120007/46‑1 (MQ=255)                              tCCTTGAGGGATGATTGCATTACATACAGATATAGCACAGGCTATATTATATAGCTa  <  1:392924/57‑1 (MQ=255)                              tCCTTGAGGGATGATTGCATTACATACAGATATAGCACAGGCTATATTATATAGCTa  <  1:536248/57‑1 (MQ=255)                              tCCTTGAGGGATGATTGCATTACATACAGATATAGCACAGGCTATATTATATAGCTa  <  1:550656/57‑1 (MQ=255)                                                        |                               AAGACCTCCTGACTTGCTAATCCCGTCGATCCTTGAGGGATGATTGCATTACATACAGATATAGCACAGGCTATATTATATAGCTA  >  minE/543915‑544000

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