Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	259,132:1	+A	100%	coding (59/348 nt)	yajD →	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	259,132	1	.	A	100.0%	125.1 / NA	36	L20* (TTA→TAA)	yajD	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (36/0);  total (36/0)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATT‑AAAAATCTATCCGTGGGTCTGCGGTCGCT  >  minE/259096‑259161                                      |                              aagcgcGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:584267/3‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTggg              >  1:242744/1‑55 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTg                >  1:1746284/1‑53 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCt           >  1:1015689/1‑58 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCg        >  1:1581795/1‑61 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:510556/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:330356/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:353393/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:353842/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:364318/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:488812/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:2048737/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:535582/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:577934/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:646532/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:666921/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:671900/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:746889/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:82119/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:892693/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:209852/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:2023576/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:1954967/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:1810026/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:1571609/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:142001/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:1351605/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:1346061/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:1321392/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:1195474/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:1095783/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:1095358/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:1076102/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:1073274/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:106344/1‑67 (MQ=255) aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAAATCTATCCGTGGGTCTGCGGTCGCt  >  1:1033536/1‑67 (MQ=255)                                      |                              ATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATT‑AAAAATCTATCCGTGGGTCTGCGGTCGCT  >  minE/259096‑259161

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