Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,589,355	C→T	V225M (GTG→ATG)	ugpE ←	sn‑glycerol 3‑phosphate ABC transporter membrane subunit UgpE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,589,355	0	C	T	100.0%	156.6 / NA	44	V225M (GTG→ATG)	ugpE	sn‑glycerol 3‑phosphate ABC transporter membrane subunit UgpE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/20);  total (24/20)

TCACTGAGTTCCATTCCGTGGTGCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCACATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACAAACAGCGCCGCCAGAT  >  NC_000913/3589277‑3589432                                                                               |                                                                              tCACTGAGTTCCATTCCGTGGTGCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCatat                                                                             <  1:638828/81‑1 (MQ=255)    cTGAGTTCCATTCCGTGGTGCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCgg                                                                          >  2:638867/1‑81 (MQ=255)       aGTTCCATTCCGTGGTGCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGta                                                                        <  1:638829/80‑1 (MQ=255)        gTTCCATTCCGTGGTGCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGtaat                                                                      <  1:638831/81‑1 (MQ=255)        gTTCCATTCCGTGGTGCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGtaat                                                                      <  1:638830/81‑1 (MQ=255)          tCCATTCCGTGGTGCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGtaataa                                                                    >  2:638901/1‑81 (MQ=255)              ttCCGTGGTGCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATcaa                                                                <  2:638832/81‑1 (MQ=255)                ccGTGGTGCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATcaac                                                               <  2:638833/80‑1 (MQ=255)                  gtggtgCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCaacaac                                                            >  1:638872/1‑81 (MQ=255)                  gtggtgCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCaacaac                                                            <  2:638834/81‑1 (MQ=255)                    ggtgCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACgg                                                          >  1:638870/1‑81 (MQ=255)                       gCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCca                                                       >  2:638886/1‑81 (MQ=255)                        ccTTCGCCTGTAGCGATCATCCCTTTGATCCCCGCCACGGTGGTGCCGAGATCCATATCGGTATTAATCAACAACGGCcac                                                      <  2:638835/81‑1 (MQ=255)                              ccTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATAc                                                >  2:638885/1‑81 (MQ=255)                               cTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACt                                               <  2:638836/81‑1 (MQ=255)                                  tAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTg                                              <  2:638838/79‑1 (MQ=255)                                  tAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTg                                              <  2:638837/79‑1 (MQ=255)                                  tAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGa                                             <  1:638839/80‑1 (MQ=255)                                      gATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTc                                          >  1:638877/1‑79 (MQ=255)                                      gATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTc                                          >  1:638904/1‑79 (MQ=255)                                      gATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCa                                        <  1:638840/81‑1 (MQ=255)                                       atcatcCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAg                                       >  1:638895/1‑81 (MQ=255)                                       atcatcCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAg                                       <  1:638841/81‑1 (MQ=255)                                            cccTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCg                                    >  1:638893/1‑79 (MQ=255)                                            cccTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCg                                    >  2:638897/1‑79 (MQ=255)                                             ccTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAg                                 <  1:638843/81‑1 (MQ=255)                                             ccTTTGATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAg                                 <  2:638842/81‑1 (MQ=255)                                                  gATCCCTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATaa                             <  1:638844/80‑1 (MQ=255)                                                     cccTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGt                         >  2:638878/1‑81 (MQ=255)                                                      ccTGCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGt                         >  1:638907/1‑80 (MQ=255)                                                         gCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGAt                      >  1:638879/1‑80 (MQ=255)                                                         gCCACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGAt                      >  1:638874/1‑80 (MQ=255)                                                           cACGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATcac                   <  2:638845/81‑1 (MQ=255)                                                            aCGGTGGTGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATcac                   >  1:638892/1‑80 (MQ=255)                                                              ggtggtGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACaaa                <  1:638846/81‑1 (MQ=255)                                                              ggtggtGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACaaa                >  1:638882/1‑81 (MQ=255)                                                              ggtggtGCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACaaa                >  2:638891/1‑81 (MQ=255)                                                                tggtgCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACaaa                >  2:638921/1‑79 (MQ=255)                                                                  gtgCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACAAACAgc            >  1:638883/1‑81 (MQ=255)                                                                  gtgCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACAAACAgc            >  2:638896/1‑81 (MQ=255)                                                                    gCCGAGATCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACAAACAgcgc          >  1:638890/1‑81 (MQ=255)                                                                        agaTCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACAAACAGcgccgc       >  1:638881/1‑80 (MQ=255)                                                                        agaTCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACAAACAGCgccgcc      <  1:638848/81‑1 (MQ=255)                                                                        agaTCCATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACAAACAGCgccgcc      <  2:638847/81‑1 (MQ=255)                                                                            ccATATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACAAACAGCGCCGCCAGAt  >  1:638906/1‑81 (MQ=255)                                                                               |                                                                              TCACTGAGTTCCATTCCGTGGTGCCTTCGCCTGTAGCGATCATCCCTTTGATCCCTGCCACGGTGGTGCCGAGATCCACATCGGTAATAATCAACAACGGCCACAAATACTGATTCCAGCCGTAGATAAAGGTGATCACAAACAGCGCCGCCAGAT  >  NC_000913/3589277‑3589432

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