Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,505,148	A→G	Y364C (TAC→TGC)	yjhB →	putative sialic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,505,148	0	A	G	100.0%	32.7 / NA	11	Y364C (TAC→TGC)	yjhB	putative sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (5/6);  total (5/6)

AGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTA  >  NC_000913/4505075‑4505224                                                                          |                                                                             agagGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATGCAttt                                                                         >  1:933037/1‑79 (MQ=255)  gagGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATGCATTTCAg                                                                      >  1:933044/1‑81 (MQ=255)         aGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATGCATTTCAGGAtatt                                                                <  1:933016/80‑1 (MQ=255)           ggACCGGTCTTATTTATAACTTAGGGGCAACAGGAGGAATGGCCGCACCTGTATTAGCTACATGCATTTCAGGATATTATg                                                             <  2:933017/81‑1 (MQ=255)                   cTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATGCATTTCAGGATATTATGGCTTAGgt                                                     >  2:933036/1‑81 (MQ=255)                      atttatAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATGCATTTCAGGATATTATGGCTTAGgt                                                     <  1:933018/78‑1 (MQ=255)                                tAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATGCATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCa                                        >  1:933039/1‑81 (MQ=255)                                                   ggCCGCACCTGTATTAGCTACATGCATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATtctc                     >  1:933041/1‑81 (MQ=255)                                                         aCCTGTATTAGCTACATGCATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCtt               <  1:933019/81‑1 (MQ=255)                                                                   gCTACATGCATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTg     <  1:933020/81‑1 (MQ=255)                                                                    cTACATGCATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGt    <  1:933021/81‑1 (MQ=255)                                                                       catgcatTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTa  <  2:933022/80‑1 (MQ=255)                                                                          |                                                                             AGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTA  >  NC_000913/4505075‑4505224

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: ‑