Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,605,436	G→T	R59S (CGC→AGC)	yhhM ←	DUF2500 domain‑containing protein YhhM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,605,436	0	G	T	100.0%	107.6 / NA	31	R59S (CGC→AGC)	yhhM	DUF2500 domain‑containing protein YhhM
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (19/12);  total (19/12)

CTCCATTCCTCCGCTTTGCGGTTTGAAGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCGCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTCACCACCACCAGCTTTTGCTGGAGCGGAGC  >  NC_000913/3605359‑3605505                                                                              |                                                                      cTCCATTCCTCCGCTTTGCGGTTTGAAGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGctcg                                                                     <  1:630592/80‑1 (MQ=255) cTCCATTCCTCCGCTTTGCGGTTTGAAGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGa                                                                    <  1:630590/81‑1 (MQ=255)  tCCATTCCTCCGCTTTGCGGTTTGAAGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGaa                                                                   >  2:630645/1‑81 (MQ=255)    cATTCCTCCGCTTTGCGGTTTGAAGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAAc                                                                  <  2:630592/80‑1 (MQ=255)         ctccGCTTTGCGGTTTGAAGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGAt                                                            >  1:630711/1‑81 (MQ=255)           ccGCTTTGCGGTTTGAAGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCg                                                          >  2:630635/1‑81 (MQ=255)             gCTTTGCGGTTTGAAGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGt                                                         >  2:630640/1‑80 (MQ=255)                    ggTTTGAAGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGt                                                 <  1:630594/81‑1 (MQ=255)                    ggTTTGAAGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGt                                                 >  2:630653/1‑81 (MQ=255)                     gTTTGAAGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGtt                                                <  2:630594/81‑1 (MQ=255)                          aaGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTc                                            <  1:630596/80‑1 (MQ=255)                              ttGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCt                                       >  2:630636/1‑81 (MQ=255)                                gCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCttg                                     >  1:630650/1‑81 (MQ=255)                                gCCTCATAGCGTATACTGGTGCCTGCCGGAGTAACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCtt                                      <  1:630597/80‑1 (MQ=255)                                  cTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCttgtt                                   >  1:630655/1‑81 (MQ=255)                                  cTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCttgtt                                   >  2:630641/1‑81 (MQ=255)                                    cATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGc                                 >  2:630638/1‑81 (MQ=255)                                    cATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGc                                 >  2:630656/1‑81 (MQ=255)                                    cATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGc                                 >  2:630657/1‑81 (MQ=255)                                      tAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTc                               <  1:630598/81‑1 (MQ=255)                                                tGGTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTcaccaccac                       >  2:630647/1‑79 (MQ=255)                                                  gTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTCACCACCACCAg                    <  1:630600/80‑1 (MQ=255)                                                  gTGCCTGCCGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTCACCACCACCAGc                   <  2:630598/81‑1 (MQ=255)                                                    gcctgccGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTCACCACCACCAGCtt                 >  1:630652/1‑81 (MQ=255)                                                    gcctgccGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCCCACCACCACCAGCtt                 >  2:630646/1‑81 (MQ=255)                                                      ctgccGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTCACCACCACCAGCtttt               <  1:630601/81‑1 (MQ=255)                                                         ccGGAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTCACCACCACCAGCTTTTg              <  2:630601/79‑1 (MQ=255)                                                            gAGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTCACCACCACCAGCTTTTGCTg           >  2:630644/1‑79 (MQ=255)                                                             aGTCACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTCACCACCACCAGCTTTTGCTggag        >  1:630643/1‑81 (MQ=255)                                                                cACTTCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTCACCACCACCAGCTTTTGCTggagcg      >  2:630650/1‑80 (MQ=255)                                                                    tCCTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTCACCACCACCAGCTTTTGCTggagcggagc  >  2:630658/1‑80 (MQ=255)                                                                     ccTGCTGGCTCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTCACCACCACCAGCTTTTGCTggagcggagc  <  2:630602/79‑1 (MQ=255)                                                                              |                                                                      CTCCATTCCTCCGCTTTGCGGTTTGAAGCTTGCCTCATAGCGTATACTGGTGCCTGCCGGAGTCACTTCCTGCTGGCGCGAACGGCGATCGTTAATCGGTTTTTCCCGCTTGTTGCTCACCACCACCAGCTTTTGCTGGAGCGGAGC  >  NC_000913/3605359‑3605505

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