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%	147.9 / NA	43	H147P (CAC→CCC)	glpT	glycerol‑3‑phosphate transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (22/21);  total (22/21)

AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTC  >  NZ_CP009273/2345327‑2345500                                                                                        |                                                                                       aaCAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  2:1052337/90‑1 (MQ=255) aaCAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  1:387420/90‑1 (MQ=255)      cagcGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCAt                                                                                 <  1:949413/90‑1 (MQ=255)         cggcggAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGt                                                                              <  2:1141971/90‑1 (MQ=255)           gcggAATACCACCCCCGACGTTGTGCGCACAGTTCCACACTGACCCAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTAc                                                                            <  1:57247/90‑1 (MQ=255)                        accGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACAcggcgg                                                               >  1:255815/1‑90 (MQ=255)                           gACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCa                                                            >  2:998138/1‑90 (MQ=255)                               ttGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCAcccc                                                        >  1:381569/1‑90 (MQ=255)                                 gtgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCAt                                                      <  1:995946/90‑1 (MQ=255)                                     gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATcccc                                                  >  2:229343/1‑90 (MQ=255)                                     gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATcccc                                                  <  2:101577/90‑1 (MQ=255)                                         aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGa                                              >  1:183172/1‑90 (MQ=255)                                          gTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaa                                             >  2:529881/1‑90 (MQ=255)                                            tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaacc                                           <  2:475073/90‑1 (MQ=255)                                            tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaacc                                           <  2:1094308/90‑1 (MQ=255)                                                   tGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCa                                    <  2:1054276/90‑1 (MQ=255)                                                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGTACCAACCGCag                                   >  2:560010/1‑90 (MQ=255)                                                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGc                                     <  2:1086632/88‑1 (MQ=255)                                                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGc                                     >  1:1086632/1‑88 (MQ=255)                                                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  1:336910/1‑90 (MQ=255)                                                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  1:979393/1‑90 (MQ=255)                                                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:838678/1‑90 (MQ=255)                                                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:841362/1‑90 (MQ=255)                                                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:397623/1‑90 (MQ=255)                                                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:366141/1‑90 (MQ=255)                                                        caATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGa                               <  2:1150652/90‑1 (MQ=255)                                                          aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCaga                                  <  1:581300/85‑1 (MQ=255)                                                          aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCaga                                  >  2:581300/1‑85 (MQ=255)                                                          aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac                             <  2:212134/90‑1 (MQ=255)                                                          aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac                             >  2:140798/1‑90 (MQ=255)                                                          aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac                             >  1:830113/1‑90 (MQ=255)                                                          aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac                             >  2:811854/1‑90 (MQ=255)                                                          aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac                             <  1:658608/90‑1 (MQ=255)                                                            gccgccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaacaa                           >  1:492709/1‑90 (MQ=255)                                                             ccgccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaacaac                          <  2:1232831/90‑1 (MQ=255)                                                               gccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAg                        >  2:22302/1‑90 (MQ=255)                                                                 cACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTa                      <  2:958537/90‑1 (MQ=255)                                                                    gTTCTTTCTGCGCCCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACaa                   <  1:658056/90‑1 (MQ=255)                                                                    gTTCTTTCTGCGACCCCCAGGGCCCCATAGTACGACCACACGGCGGCCCCCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACaa                   <  2:734634/90‑1 (MQ=255)                                                                         ttctGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATc              <  2:979393/90‑1 (MQ=255)                                                                             gCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCg          >  2:517893/1‑90 (MQ=255)                                                                               gACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCa        >  2:741632/1‑90 (MQ=255)                                                                                     caGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTc  <  2:1082268/90‑1 (MQ=255)                                                                                        |                                                                                       AACAGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTC  >  NZ_CP009273/2345327‑2345500

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