Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	551,066	G→T	A210S (GCC→TCC)	ybcF →	putative carbamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	551,066	0	G	T	100.0%	125.8 / NA	36	A210S (GCC→TCC)	ybcF	putative carbamate kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (20/16);  total (20/16)

GGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAA  >  NC_000913/550988‑551135                                                                               |                                                                      ggcggcggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTcc                                                                     <  2:106019/81‑1 (MQ=255) ggcggcggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTcc                                                                     >  1:106055/1‑81 (MQ=255) ggcggcggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGCGATTGATAAAGATCTCGCCGCTGCGTTGCTCTcc                                                                     >  1:106062/1‑81 (MQ=255)  gcggcggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCg                                                                    <  1:106021/81‑1 (MQ=255)  gcggcggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCg                                                                    <  1:106022/81‑1 (MQ=255)  gcggcggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCg                                                                    <  1:106020/81‑1 (MQ=255)   cggcggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTcc                                                                     >  1:106050/1‑79 (MQ=255)    ggcggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAg                                                                  <  1:106023/81‑1 (MQ=255)       ggTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAg                                                               <  2:106024/81‑1 (MQ=255)          gtgCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGAt                                                             >  2:106070/1‑80 (MQ=255)            gCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTaa                                                          >  2:106072/1‑81 (MQ=255)                  gacggatgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATg                                                        >  1:106068/1‑77 (MQ=255)                     ggatgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATgg                                                 <  1:106026/81‑1 (MQ=255)                     ggatgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATg                                                  >  2:106069/1‑80 (MQ=255)                      gatgacggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGa                                                >  2:106051/1‑81 (MQ=255)                          acggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTgg                                            <  2:106027/81‑1 (MQ=255)                           cggaGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGt                                           <  1:106028/81‑1 (MQ=255)                              aGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGAt                                        >  2:106059/1‑81 (MQ=255)                                cAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATcc                                      >  2:106081/1‑81 (MQ=255)                                 aGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCt                                     <  2:106029/81‑1 (MQ=255)                                  gggAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTc                                    <  2:106030/81‑1 (MQ=255)                                    gAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCAc                                  <  2:106031/81‑1 (MQ=255)                                     aGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCAcc                                 >  1:106064/1‑81 (MQ=255)                                      gTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCAcc                                 >  1:106088/1‑80 (MQ=255)                                        gAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGAt                              >  1:106057/1‑81 (MQ=255)                                         aaGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGAATAATGCAGATGGACTGGTGATCCTCACCGAt                              <  2:106032/80‑1 (MQ=255)                                                 attgatAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGc                      <  1:106033/80‑1 (MQ=255)                                                     atAAAGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGta                  <  2:106034/80‑1 (MQ=255)                                                        aaGATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTatat               >  1:106076/1‑80 (MQ=255)                                                          gATCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGa             >  2:106058/1‑80 (MQ=255)                                                           aTCTCGCCGCTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGaa            >  1:106067/1‑80 (MQ=255)                                                                gccgcTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTgg      >  2:106087/1‑81 (MQ=255)                                                                 ccgcTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTggg     >  2:585747/1‑81 (MQ=255)                                                                  cgcTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTgggg    <  2:106035/81‑1 (MQ=255)                                                                  cgcTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTggg     >  2:106074/1‑80 (MQ=255)                                                                    cTGCGTTGCTCTCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGaa  >  1:106080/1‑81 (MQ=255)                                                                               |                                                                      GGCGGCGGTGTGCCTGTGACGGATGACGGAGCAGGGAGTGAAGCAGTGATTGATAAAGATCTCGCCGCTGCGTTGCTCGCCGAGCAGATTAATGCAGATGGACTGGTGATCCTCACCGATGCTGATGCGGTATATGAAAACTGGGGAA  >  NC_000913/550988‑551135

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