Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,345,414	T→G	H147P (CAC→CCC)	glpT ←	glycerol‑3‑phosphate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,345,414	0	T	G	100.0%	101.0 / NA	31	H147P (CAC→CCC)	glpT	glycerol‑3‑phosphate transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (15/16);  total (15/16)

AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAA  >  NZ_CP009273/2345327‑2345483                                                                                        |                                                                      aaCAGCAGCGGCGGAATACCACCCCCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                     <  2:648737/90‑1 (MQ=255) aaCAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                     <  2:651145/90‑1 (MQ=255) aaCAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                     <  1:40582/90‑1 (MQ=255)    agcagcGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCAcc                                                                  <  2:628176/90‑1 (MQ=255)      cagcGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCAt                                                                <  2:191014/90‑1 (MQ=255)           gcggAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTAc                                                           <  1:440335/90‑1 (MQ=255)                         ccGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACggcggc                                             <  1:206001/90‑1 (MQ=255)                          cGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGcc                                            <  1:142188/90‑1 (MQ=255)                          cGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGcc                                            <  2:147519/90‑1 (MQ=255)                                   gcgcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATcc                                   >  2:207808/1‑90 (MQ=255)                                          gTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaa                            >  1:248526/1‑90 (MQ=255)                                          gTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaa                            >  1:289963/1‑90 (MQ=255)                                               acacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaaccaac                       >  1:124644/1‑90 (MQ=255)                                                  cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTAc                                                           <  1:460040/51‑1 (MQ=255)                                                  cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTAc                                                           >  2:460040/1‑51 (MQ=255)                                                  cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGc                    >  1:598888/1‑90 (MQ=255)                                                    gccacaATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                  >  2:229861/3‑90 (MQ=255)                                                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                  >  2:108891/1‑90 (MQ=255)                                                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                  >  2:511153/1‑90 (MQ=255)                                                    gACACAATGCC‑CCGCGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCaga                 <  2:465787/90‑1 (MQ=255)                                                     acacAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCaga                 <  1:108891/90‑1 (MQ=255)                                                     acacAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCaga                 <  2:61527/90‑1 (MQ=255)                                                        caATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAATCGCAGAGGa              <  2:235242/90‑1 (MQ=255)                                                          aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCAt                                     <  1:370361/65‑1 (MQ=255)                                                          aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCAt                                     >  2:370361/1‑65 (MQ=255)                                                          aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac            >  1:297500/1‑90 (MQ=255)                                                              cgccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAacaaca        >  1:51562/1‑90 (MQ=255)                                                              cgccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAacaaca        >  1:16920/1‑90 (MQ=255)                                                                    gTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCAccc                                        >  1:292568/1‑52 (MQ=255)                                                                    gTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCAccc                                        <  2:292568/52‑1 (MQ=255)                                                                    gTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACaa  >  1:407671/1‑90 (MQ=255)                                                                                        |                                                                      AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAA  >  NZ_CP009273/2345327‑2345483

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