Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,426,573:1	+A	100%	coding (558/864 nt)	napH ←	ferredoxin‑type protein essential for electron transfer from ubiquinol to periplasmic nitrate reductase (NapAB)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,426,573	1	.	A	75.0%	14.9 / 2.0	8	coding (558/864 nt)	napH	ferredoxin‑type protein essential for electron transfer from ubiquinol to periplasmic nitrate reductase (NapAB)
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (0/6);  total (0/8)
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
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

ACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAAT‑AAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCA  >  minE/1426523‑1426616                                                    |                                            aCCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAAT‑AAAAACAGCGCGAGaa                             <  1:2101684/67‑1 (MQ=255) aCCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAAT‑AAAAACAGCGCGAGaa                             <  1:860915/67‑1 (MQ=255)                             gtTCAACGACCAGTAGATCAAATAAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCa  <  1:1138969/67‑1 (MQ=255)                             gtTCAACGACCAGTAGATCAAATAAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCa  <  1:1423143/67‑1 (MQ=255)                             gtTCAACGACCAGTAGATCAAATAAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCa  <  1:1979288/67‑1 (MQ=255)                             gtTCAACGACCAGTAGATCAAATAAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCa  <  1:2163700/67‑1 (MQ=255)                             gtTCAACGACCAGTAGATCAAATAAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCa  <  1:2183482/67‑1 (MQ=255)                             gtTCAACGACCAGTAGATCAAATAAAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCa  <  1:316295/67‑1 (MQ=255)                                                    |                                            ACCGGACAAATGTGCCCGCACCAGCCGTGTTCAACGACCAGTAGATCAAAT‑AAAAACAGCGCGAGAATAAGCAGCGCGCCGCTGCCGAAGCCCA  >  minE/1426523‑1426616

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