Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	219,917	T→C	64.6%	P737P (CCA→CCG)	fadE ←	acyl coenzyme A dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	219,917	0	T	C	64.6%	18.0 / 14.4	31	P737P (CCA→CCG)	fadE	acyl coenzyme A dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (5/6);  new base C (16/4);  total (21/10)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTAAACGGCAGGTTTTTACCCAGCTCTTTACAGATCCGCTGATGAATTGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGG  >  minE/219868‑219982                                                  |                                                                  gggTAAACGGCAGGTTTTTACCCAGCTCTTTACAGATCCGCTGATGAATTGGGTCGGCGGCAATCACAt                                                <  1:845680/69‑1 (MQ=255) gggTAAACGGCAGGTTTTTACCCAGCTCTTTACAGATCCGCTGATGAATTGGGTCGGCGGCAATCACAt                                                <  1:1127943/69‑1 (MQ=255) gggTAAACGGCAGGTTTTTACCCAGCTCTTTACAGATCCGCTGATGAATTGGGTCGGCGGCAATCACAt                                                <  1:3397680/69‑1 (MQ=255)    tAAACGGCAGGTTTTTACCCAGCTCTTTACAGATCCGCTGATGAATTGGGTCGGCGGCAATCACATcc                                              >  1:1943226/1‑68 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCTGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:1869926/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:104499/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:983382/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:943138/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:892634/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:700054/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:517537/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:3177033/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:3025272/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:2890615/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:2749094/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:2503310/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:2043143/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:1124850/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:1066803/1‑67 (MQ=255)                      ccAACTCTTTACAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAgta                           >  1:1047230/1‑67 (MQ=255)                               tACAGATCCGCTGATGAATTGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTc                                >  1:462926/1‑55 (MQ=255)                               tACAGATCCGCTGATGAATTGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCg                  >  1:3165136/1‑69 (MQ=255)                               tACAGATCCGCTGATGAATTGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCg                  >  1:2161877/1‑69 (MQ=255)                                 cAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAGTAAACCAACCgg                 <  1:2557885/68‑1 (MQ=255)                                 cAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAGTAAACCAACCgg                 <  1:3198976/68‑1 (MQ=255)                                 cAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAGTAAACCAACCgg                 <  1:1660060/68‑1 (MQ=255)                                 cAGATCCGCTGATGAATCGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAGTAAACCAACCgg                 <  1:1014694/68‑1 (MQ=255)                                      ccGCTGATGAATTGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCg                  <  1:2351282/62‑1 (MQ=255)                                      ccGCTGATGAATTGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCg                  <  1:3090598/62‑1 (MQ=255)                                      ccGCTGATGAATTGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCg                  <  1:107180/62‑1 (MQ=255)                                               aaTTGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCgg  >  1:2409679/1‑69 (MQ=255)                                                  |                                                                  GGGTAAACGGCAGGTTTTTACCCAGCTCTTTACAGATCCGCTGATGAATTGGGTCGGCGGCAATCACATCCACCAGCGCCTCTTCCAGCAAGCCAACCGGATTATGCTCGCTCGG  >  minE/219868‑219982

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