Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	551,138	C→T	P234S (CCG→TCG)	ybcF →	putative carbamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	551,138	0	C	T	100.0%	104.1 / NA	33	P234S (CCG→TCG)	ybcF	putative carbamate kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/24);  total (9/24)

CGGGTGGTGGCGTCTCCGCAACCGCGTAAAATTCTCGACAGCGAAGCCATCGAGTTGTTGCTCAAAGAGGGGCATGTGGTGATTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGCCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTAT  >  NC_000913/550895‑551381                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     cGGGTGGTGGCGTCTCCGCAACCGCGTAAAATTCTCGACAGCGAAGCCATCGAGTTGTTGCTCAAAGAGGGGCATGTGGTGATTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTGACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCgcagca                                                                                                                                                                                                                                               <  1:20442/251‑1 (MQ=255)   ggtggtGGCGTCTCCGCAACCGCGTAAAATTCTCGACAGCGAAGCCATCGAGTTGGTGCTCAAAGAGGGGCATGTGGTGATTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAAc                                                                                                                                                                                                                                             <  2:284153/251‑1 (MQ=255)                   cAACCGCGTAAAATTCTCGACAGCGAAGCCATCGAGTTGTTGCTCAAAGAGGGGCATGTGGTGATTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCAt                                                                                                                                                                                                                               <  1:5150/249‑1 (MQ=255)                         cgTAAAATTCTCGACAGCGAAGCCATCGAGTTGTTGCTCAAAGAGGGGCATGTGGTGATTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCGCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGTGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACAcc                                                                                                                                                                                                                       <  1:310943/251‑1 (MQ=255)                                           gAAGCCATCGAGTTGTTGCTCAAAGAGGGGCATGTGGTGATTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGGGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCCTGCCACACC‑GGATGAGTTAGCGCCAtt                                                                                                                                                                                                    >  2:64450/1‑251 (MQ=255)                                            aaGCCATCGAGTTGTTGCTCAAAGAGGGGCATGTGGTGATTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCAttt                                                                                                                                                                                                   <  1:85039/251‑1 (MQ=255)                                                   cGAGTTGTTGCTCAAAGAGGGGCATGTGGTGATTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCaaa                                                                                                                                                                                             <  2:261667/250‑1 (MQ=255)                                                   cGAGTTGTTGCTCAAAGAGGGGCATGTGGTGATTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCaa                                                                                                                                                                                              >  1:372839/1‑249 (MQ=255)                                                       ttgttgCTCAAAGAGGGGCATGTGGTGATGTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCCGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCCCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCCCCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCAAGCC‑CACCCGGGTGGGTTAGCGCCATTTGCCAAAACCgg                                                                                                                                                                                        >  2:278147/1‑250 (MQ=255)                                                               ctaagagggggtggtggtgAGTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGAGTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTc                                                                                                                                                                                  <  2:377842/236‑1 (MQ=255)                                                               cAAAGAGGGGCATGTGGTGATTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGAGGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTc                                                                                                                                                                                  <  2:377832/249‑1 (MQ=255)                                                                aaaGAGGGGCATGTGGTGATTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGCTTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGAt                                                                                                                                                                               >  2:167027/1‑251 (MQ=255)                                                                    agGGGCATGAGGTGATTTGCAGTAGCGGCGGGGGGGTGCCTGTGGCGGATGACGGAGAAGGGGGGGAAGCCGGGATTGATAAAGATATCGCCGCTGCGGTGCTCGCAGAGCAGAGTAATGAAGATGGACTGGTGATCCTCACCGATGCTGATGAGGTATATGAAAACTGTGGAACGTCGCAGCAACGTGCCACTTTCAATTCCACACC‑GGAGGAGTTAAGGCCAttttcacaatcggaagggtcgatgggg                                                                                                                                                                           >  2:335852/1‑227 (MQ=255)                                                                                             cggcggcggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGtt                                                                                                                                                  <  1:215444/251‑1 (MQ=255)                                                                                                    ggTGTGCCTGTGACGGATGACGGAGCCGGGAGTGAAGCAGTGTTTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCCGATGAACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTTCCATTCGCCATGCCACACC‑GGATGCGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGTATGTAAcggcggggggggggtttgtca                                                                                                                                            >  2:70056/1‑235 (MQ=255)                                                                                                     gtgtGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTGCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGa                                                                                                                                          <  1:367995/251‑1 (MQ=255)                                                                                                         gCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGGGGGTTTGTCAGAAGcc                                                                                                                                      >  1:86426/1‑251 (MQ=255)                                                                                                                         ggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTg                                                                                                                      <  1:344777/251‑1 (MQ=255)                                                                                                                                      aaGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGGGGGTATGTCAGAAGACGGGGGAAACCCGCGGGGGTTgggggggtt                                                                                                         >  1:48747/1‑246 (MQ=255)                                                                                                                                        gCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATc                                                                                                       <  1:222973/251‑1 (MQ=255)                                                                                                                                        gCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCGCGCCGAGCAGATTAATGCAGATGGACGGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATc                                                                                                       <  2:147730/251‑1 (MQ=255)                                                                                                                                                      aGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCCCCGATGCTGAGGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAgaga                                                                                         <  2:359237/251‑1 (MQ=255)                                                                                                                                                                   gCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGa                                                                            <  1:187652/251‑1 (MQ=255)                                                                                                                                                                       tGCTCGCCGAGCAGATTAATGCAGATGGACTGCTGATCCTCACCGATGCTGCTGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGCGGGGTTTTGTCCGAAGCCGTGGGAAACCCGCGTGGGGTGGGGCGTTTTCGCGCATTGAAGAGACGCTggggggcggaggg                                                                        >  1:235204/1‑246 (MQ=255)                                                                                                                                                                                              gatggccgggTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATTCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTTGTAGACCCGCGTGGATTTGGGCCCTATCGCTAAATGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTc                                                 <  1:226019/243‑1 (MQ=255)                                                                                                                                                                                                        ggATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTa                                       <  2:202002/250‑1 (MQ=255)                                                                                                                                                                                                               cACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACAttt                                <  1:210045/251‑1 (MQ=255)                                                                                                                                                                                                               cACCGATGCTGATGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACAtt                                 <  2:276424/250‑1 (MQ=255)                                                                                                                                                                                                                          aTGCGGTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGc                     <  1:183898/251‑1 (MQ=255)                                                                                                                                                                                                                              ggTATATGAAAACTGGGGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACg                 <  2:136905/251‑1 (MQ=255)                                                                                                                                                                                                                                 atatGAAAACTGGTGAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCTGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATTCGACGCt               <  2:233728/250‑1 (MQ=255)                                                                                                                                                                                                                                            ggggAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCt     <  2:84223/249‑1 (MQ=255)                                                                                                                                                                                                                                             gggAACGTCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTAt  <  1:34593/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     CGGGTGGTGGCGTCTCCGCAACCGCGTAAAATTCTCGACAGCGAAGCCATCGAGTTGTTGCTCAAAGAGGGGCATGTGGTGATTTGCAGTGGCGGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAACGCCGCAGCAACGTGCCATTCGCCATGCCACACC‑GGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTAT  >  NC_000913/550895‑551381

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