Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,250,940	(T)6→7	coding (353/522 nt)	ruvC ←	component of RuvABC resolvasome, endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,250,934	1	.	T	85.0%	105.4 / 17.1	40	coding (359/522 nt)	ruvC	component of RuvABC resolvasome, endonuclease
Reads supporting (aligned to +/- strand):  ref base . (0/6);  new base T (0/34);  total (0/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

ATATGCTGCACCTGGC‑TTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACG  >  minE/1250919‑1250980                 |                                               atatGCTGCACCTGGC‑TTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:1815207/62‑1 (MQ=255) atatGCTGCACCTGGC‑TTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2621967/62‑1 (MQ=255) atatGCTGCACCTGGC‑TTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:279772/62‑1 (MQ=255) atatGCTGCACCTGGC‑TTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2958050/62‑1 (MQ=255) atatGCTGCACCTGGC‑TTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:317690/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:528259/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:33264/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:3186690/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:797552/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:3136046/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:91668/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:3120814/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2959206/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:806157/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2927703/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2865027/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2837350/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:815414/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2726212/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2665555/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:1156130/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:1813545/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:1255895/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:1322800/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:1459699/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:1523549/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:1564944/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:1788273/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:1802602/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2603521/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2022824/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2146863/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:215918/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2308421/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2384746/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2390600/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2444947/62‑1 (MQ=255)  tatGCTGCACCTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:2520606/62‑1 (MQ=255)   atGCTGCACCTGGC‑TTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:1519562/60‑1 (MQ=255)            cTGGCTTTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACg  <  1:3192163/52‑1 (MQ=255)                 |                                               ATATGCTGCACCTGGC‑TTTTTTCGGCACTGCCAATACCTACCACCGTTTGCTTTACCTGACG  >  minE/1250919‑1250980

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