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

CATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTT  >  NZ_CP009273/3896548‑3896707                                                                          |                                                                                       cATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGAt                                                                        <  1:962616/90‑1 (MQ=255)           cAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAAcc                                                              >  1:1232676/1‑90 (MQ=255)           cAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAAcc                                                              >  1:681210/1‑90 (MQ=255)           cAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAAcc                                                              >  1:799905/1‑90 (MQ=255)           cAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGAGGGCGTAAGGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAAcc                                                              >  1:515153/1‑90 (MQ=255)                 tAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTa                                                        >  2:569787/1‑90 (MQ=255)                               ttgggTAGGGAAAAAAGCCCGGGGGGGCGAAATGCGCGAGGTGGTCGATCCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGcc                                          <  2:1155580/87‑1 (MQ=255)                                  ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  1:512388/1‑90 (MQ=255)                                  ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  2:1135477/1‑90 (MQ=255)                                  ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                       >  2:1052620/1‑90 (MQ=255)                                    tggCGAGAAAAGCCCGTGGGGGGTTAATGCTCGAAGTGGTCGATGCCCCCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCatat                                     <  2:1077172/88‑1 (MQ=255)                                       cGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGc                                  >  1:180542/1‑90 (MQ=255)                                            aaaGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTc                             >  2:928057/1‑90 (MQ=255)                                                cccGGGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTaa                         <  2:436985/89‑1 (MQ=255)                                                  cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                       >  2:238734/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:363515/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:794690/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:652370/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:84796/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:628182/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:424653/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:419672/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:1032886/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:336840/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:264491/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  2:169437/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:421102/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:432490/1‑90 (MQ=255)                                                   gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                      >  1:658738/1‑90 (MQ=255)                                                    tGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATggg                     >  1:88292/1‑90 (MQ=255)                                                       gggCGTAATCCCCGAGGGGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATa                  <  2:432490/90‑1 (MQ=255)                                                                  tCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGc       >  2:1211402/1‑90 (MQ=255)                                                                  tCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGc       >  2:637952/1‑90 (MQ=255)                                                                     aaGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGc    >  2:959184/1‑90 (MQ=255)                                                                       gTGGTCGATCCCCGCTCGTTTCCCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCtt  <  1:652370/90‑1 (MQ=255)                                                                          |                                                                                       CATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGCGGCTT  >  NZ_CP009273/3896548‑3896707

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