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

CAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAG  >  NZ_CP009273/3896558‑3896701                                                                |                                                                                 cAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAAcc                                                        >  1:358039/1‑90 (MQ=255)       tAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTa                                                  >  2:351427/1‑90 (MQ=255)             cACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGttttt                                            >  2:19540/1‑90 (MQ=255)                 ttttGTTGGTAGCGAGAAAAGCCCGGGGGGGCTTATTGCTCGAGGGGGCCGAGGCCCGCTCGTTTACCCCACCCGCCGTAGTTTTTCAcc                                        <  1:470780/87‑1 (MQ=255)                 ttATGTTGGTGGGGAGTAAAGCCGGGGGGGGCGTAAGGCCGGAAGTGGTCGAGGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCAcc                                        <  2:155187/90‑1 (MQ=255)                   aTGTTGGTAGCGAGAAAAGCCCGGGGGGGGGTAATGCCCGAAGTGGTCGATGACCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAg                                      <  2:111947/90‑1 (MQ=255)                    tGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGc                                     >  2:143778/1‑90 (MQ=255)                      ttGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGccc                                   >  1:96398/1‑90 (MQ=255)                         gTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCata                                >  2:233708/1‑90 (MQ=255)                         gTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCata                                >  2:440661/1‑90 (MQ=255)                            gCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATgg                             >  2:439540/1‑90 (MQ=255)                             cGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGc                            >  1:111947/1‑90 (MQ=255)                              gAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCa                           >  1:604569/1‑90 (MQ=255)                                      cccGTGGGGGTGTAATGCTGGAGGTGGTCGATGCCCGCTCGATTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTaa                   <  2:258983/89‑1 (MQ=255)                                        cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                 >  2:140845/1‑90 (MQ=255)                                        cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                 >  2:282654/1‑90 (MQ=255)                                        cGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATg                 >  2:467978/1‑90 (MQ=255)                                         gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                >  1:258983/1‑90 (MQ=255)                                         gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                >  2:33535/1‑90 (MQ=255)                                         gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                >  1:499546/1‑90 (MQ=255)                                          tGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATggg               >  1:652180/1‑90 (MQ=255)                                          tGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATggg               >  2:552152/1‑90 (MQ=255)                                           ggggggGTAATGCCCGAGGTGGTCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGa              <  2:653570/88‑1 (MQ=255)                                                       cccGAAGTGGTCGATGCCCGCTCGATTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAg  <  2:572981/88‑1 (MQ=255)                                                                |                                                                                 CAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAG  >  NZ_CP009273/3896558‑3896701

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