Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,896,621	A→G	D147D (GAT→GAC)	bglB ←	6‑phospho‑beta‑glucosidase BglB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,896,621	0	A	G	100.0%	30.5 / NA	12	D147D (GAT→GAC)	bglB	6‑phospho‑beta‑glucosidase BglB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/4);  total (8/4)

ACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGG  >  NZ_CP009273/3896536‑3896687                                                                                      |                                                                   aCATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGa                                                                >  1:83415/1‑90 (MQ=255)            tCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGa                                                     <  1:237397/90‑1 (MQ=255)                                   cACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGttttt                              >  1:104673/1‑90 (MQ=255)                                       ttATGTTGGTAGCGAGAAAAGCCCGGGGGGGCGTAATGCTCGAAGTGGTCGAGGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCAcc                          <  2:104673/90‑1 (MQ=255)                                                     aaaaaaGCCCGGGGGGGCTTATTCCTGGAAGGGGTCGTTGCCCGTTCGTTTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCAt            <  1:303039/87‑1 (MQ=255)                                                             ccGGGCGGGCGTAATGCCCGAAGTGGTCGATGCCCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAAt    <  2:58694/90‑1 (MQ=255)                                                              cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg   >  1:488599/1‑90 (MQ=255)                                                              cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg   >  2:189608/1‑90 (MQ=255)                                                               gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg  >  1:339002/1‑90 (MQ=255)                                                               gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg  >  1:457816/1‑90 (MQ=255)                                                               gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg  >  2:16840/1‑90 (MQ=255)                                                               gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg  >  2:303039/1‑90 (MQ=255)                                                                                      |                                                                   ACATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGG  >  NZ_CP009273/3896536‑3896687

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑