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%	64.9 / NA	22	D147D (GAT→GAC)	bglB	6‑phospho‑beta‑glucosidase BglB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/6);  total (16/6)

CGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTT  >  NZ_CP009273/3896555‑3896707                                                                   |                                                                                       cGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCa                                                                 >  2:223447/1‑90 (MQ=255)    cAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                                             >  2:475866/1‑75 (MQ=255)    cAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGc                                                                             <  1:475866/75‑1 (MQ=255)    cAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAAcc                                                              >  1:317744/1‑90 (MQ=255)          tAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATCAGCCCAACCGCCGTa                                                        >  2:310676/1‑90 (MQ=255)                 aCTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTc                                                 >  1:311136/1‑90 (MQ=255)                 aCTTTATGTTGGTAGCGAGTAAAGACCGCGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTc                                                 >  1:196219/1‑90 (MQ=255)                           ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  1:156392/1‑90 (MQ=255)                           ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  1:195404/1‑90 (MQ=255)                           ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  2:276687/1‑90 (MQ=255)                              aGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATg                                    >  1:104950/1‑90 (MQ=255)                                   gaaaaGCCGGTGGGGGCATAATGCTCGAAGTGGTCGATCCCCCCTCGATTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCAtt                               <  2:311136/88‑1 (MQ=255)                                          ccGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAAt                        >  2:96437/1‑90 (MQ=255)                                           cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  1:198569/1‑90 (MQ=255)                                            gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:173640/1‑90 (MQ=255)                                            gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:255349/1‑90 (MQ=255)                                            gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:31953/1‑90 (MQ=255)                                            gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:330140/1‑90 (MQ=255)                                                gggCGTAGTGCCCGAGGTGGTGAATCCCCGCTCTATTACCCAAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATa                  <  1:348311/90‑1 (MQ=255)                                                gggCGTAATGCTCGAAGTGGTGGATGCCCGCTCGATTAGCCAAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATa                  <  1:223080/90‑1 (MQ=255)                                                 ggCGTATTGCCGGAAGGGGTCGATGCCCGCTTGATAACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATaa                 <  2:422438/90‑1 (MQ=255)                                                                gTGGTCGATGCCCGCTCTATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCtt  <  2:127919/90‑1 (MQ=255)                                                                   |                                                                                       CGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTT  >  NZ_CP009273/3896555‑3896707

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