Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	minE	1,269,468	(A)7→8	intergenic (‑18/‑282)	uvrC ← / → yedP	excinuclease UvrABC, endonuclease subunit/conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,269,461	1	.	A	100.0%	43.3 / NA	15	intergenic (‑11/‑289)	uvrC/yedP	excinuclease UvrABC, endonuclease subunit/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (14/1);  total (14/1)

GGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGGAATGCACAAT‑AAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCTTAAGCAAGAAGCACTTAAAAAAT  >  minE/1269414‑1269525                                                 |                                                                 ggTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGGAATGCACAATAAAAAAAATCCCGACCCTGAggg                                            >  1:1298614/1‑71 (MQ=255)                     aaCGCGCCCTATCGCCGGAATGCACAATAAAAAAAATCCCGACCCTGAggg                                            >  1:3965595/1‑51 (MQ=255)                     aaCGCGCCCTATCGCCGGAATGCACAATAAAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCt                         >  1:16259/1‑70 (MQ=255)                     aaCGCGCCCTATCGCCGGAATGCACAATAAAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCt                         >  1:2208351/1‑70 (MQ=255)                     aaCGCGCCCTATCGCCGGAATGCACAATAAAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCt                         >  1:2580850/1‑70 (MQ=255)                     aaCGCGCCCTATCGCCGGAATGCACAATAAAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCt                         >  1:793030/1‑70 (MQ=255)                     aaCGCGCCCTATCGCC‑GAATGCACAATAAAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCtt                        >  1:4332274/1‑70 (MQ=255)                          gcCCTATCGCCGGAATGCACAATAAAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCTTAAGc                    <  1:4355639/70‑1 (MQ=255)                                        aTGCACAATAAAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCTTAAGCAAGAAGCACTTaaa      >  1:1269309/1‑70 (MQ=255)                                        aTGCACAATAAAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCTTAAGCAAGAAGCACTTaaa      >  1:3330499/1‑70 (MQ=255)                                        aTGCACAATAAAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCTTAAGCAAGAAGCACTTaaa      >  1:580876/1‑70 (MQ=255)                                        aTGCACAATAAAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCTTAAGCAAGAAGCACTTaaa      >  1:774574/1‑70 (MQ=255)                                        aTGCACAATAAAAAAAATCCAGACCCTGAGGGGGTCGGGATGAAACTTGCTTAAGCAAGAAGCACTTaaa      >  1:2533850/1‑70 (MQ=255)                                           cacaATAAAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCTTAAGCAAGAAGCACTTAAAAAAt  >  1:4087229/1‑71 (MQ=255)                                           cacaATAAAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCTTAAGCAAGAAGCACTTAAAAAAt  >  1:4087232/1‑71 (MQ=255)                                                 |                                                                 GGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGGAATGCACAAT‑AAAAAAATCCCGACCCTGAGGGGGTCGGGATGAAACTTGCTTAAGCAAGAAGCACTTAAAAAAT  >  minE/1269414‑1269525

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