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	88.0%	49.7 / 1.2	25	D147D (GAT→GAC)	bglB	6‑phospho‑beta‑glucosidase BglB
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (14/8);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

ATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGC  >  NZ_CP009273/3896538‑3896702                                                                                    |                                                                                  aTGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATg                                                                             >  1:823186/1‑90 (MQ=255)  tGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGa                                                                            >  1:507785/1‑90 (MQ=255)       aTTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCt                                                                       <  1:825822/90‑1 (MQ=255)                                  aCTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTc                                            >  1:308095/1‑90 (MQ=255)                                      tttGTTGGTGGCGAAAAAAGCCCGGGGGGGCGTAATGCTCGAAGGGGTCGATGCCCCCTCGATTACCCCACCCGCCGTAGTTTTTCACCa                                        <  2:280019/88‑1 (MQ=255)                                      tATGTTGGTAGCGAGAAAAGCCCGTGCGGGCGTAATGCTCGAAGTGGTCGATGCCCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCa                                        <  1:230379/90‑1 (MQ=255)                                       aTGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAg                                       >  2:545594/1‑90 (MQ=255)                                            ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                  >  2:78876/1‑90 (MQ=255)                                            ggTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCat                                  >  1:362832/1‑90 (MQ=255)                                               aGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATg                               >  2:382378/1‑90 (MQ=255)                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                 >  1:552088/1‑90 (MQ=255)                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                 >  1:733175/1‑90 (MQ=255)                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                 >  1:778404/1‑90 (MQ=255)                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                 >  2:58084/1‑90 (MQ=255)                                                             gTGCGGGCGTAATGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                 >  1:470903/1‑90 (MQ=255)                                                             gTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATgg                 >  2:364104/1‑90 (MQ=255)                                                               gggggCGAAATGCTCGAAGTGGTCGATGCCCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGa               <  1:5570/88‑1 (MQ=255)                                                               gggggCGAAATCCTCGAAGTGGTCGAGGCCCGCTCAATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGa               <  1:798195/88‑1 (MQ=255)                                                               gCGGGCGTAATGCTGGAGGGGGCCGATGCCCGCCCTATTACCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGa               <  2:788683/90‑1 (MQ=255)                                                                 gtgCGTAATGCTCGGCGTGGGCGATCCCCCCTGGATCACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATa             <  2:733175/88‑1 (MQ=255)                                                                 gggCGTAATGCTCGAAGTGATCGATGCCCGCTCGATTACCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATa             <  1:194602/90‑1 (MQ=255)                                                                 gggCGAAATGCGGGAGGTGGTGGTTGCCCGCTCGATTAGCCCACCCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATa             <  1:433762/90‑1 (MQ=255)                                                                   gCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAAc           >  2:683577/1‑90 (MQ=255)                                                                        aTGCTCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTAc      >  2:95730/1‑90 (MQ=255)                                                                            tCGAAGTGGTCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGc  >  2:746954/1‑90 (MQ=255)                                                                                    |                                                                                  ATGTTGATTTCATTAAACGTCAGCCATAACGCCACTTTATGTTGGTAGCGAGTAAAGACCGTGCGGGCGTAATGCTCGAAGTGATCGATGACCGCTCGATTAGCCCAACCGCCGTAGTTTTTCACCAGCCCATATGGCATTTCGTAATGGGATAACGTTACCAGC  >  NZ_CP009273/3896538‑3896702

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