Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,604,716	C→A	A186E (GCA→GAA)	bglJ →	DNA‑binding transcriptional regulator BglJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,716	0	C	A	100.0%	12.1 / NA	7	A186E (GCA→GAA)	bglJ	DNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (3/4);  total (3/4)

ATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGC  >  NC_000913/4604649‑4604826                                                                    |                                                                                                               aTGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCCATATCAAAACGATCCGTGAACATAAATTTAAAGTATGGTCGAAACTGGGCGTCCATTCTCAC‑CGGGGTTGTTTGGGGccc                                                  >  1:80124/1‑128 (MQ=255)                     gACACAAATTGCCGGGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatgtat                               <  1:304201/129‑1 (MQ=255)                           aaTTGCCGAACA‑CTTAAAACCAAATACAAAACGGTCCCTGAACATAAATTTAATGGGAGTTCGAAACTGGGCGGCATTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCAt                        <  1:472216/129‑1 (MQ=255)                            aTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGCCAGTTCTACCCCAGGGTGTTTGGGGCCCCCAGAT‑TTCCGTTTTTTATGGGGCAt                        >  1:213255/1‑128 (MQ=255)                             ttGCCGGGCAGCTTAAACGCAAAATCAAAACGATCCGTGAACATAAATTTAATGGGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCAtt                       <  1:490243/129‑1 (MQ=255)                                                 aaTATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTTGTTATGCGGCATTGCGAAACAAGAAATGTGttc   >  1:279198/1‑128 (MQ=255)                                                  aaaacAAAACGGTCCGTGAACATAAATTTAATGGGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGc  <  1:330908/125‑1 (MQ=255)                                                                    |                                                                                                               ATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGC  >  NC_000913/4604649‑4604826

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑