Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,504,957	(T)8→7	coding (900/1218 nt)	yjhB →	putative sialic acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,504,950	0	T	.	100.0%	81.3 / NA	20	coding (893/1218 nt)	yjhB	putative sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (15/5);  total (15/5)

GATAAGATTGGTGTAAAAAAAGCCTTTGTAGTCGGTCTAATAACTTCATTTATTTTCCTTTGTCCTCTTTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCTTTGGATTAATGTTTACAAATTTAGGTATT  >  NC_000913/4504883‑4505026                                                                    |                                                                             gATAAGATTGGTGTAAAAAAAGCCTTTGTAGTCGGTCTAATAACTTCATTTATTTTCCTTTGTCCTC‑TTTTTTTATTTCtg                                                                >  1:862945/1‑81 (MQ=255)     aGATTGGTGTAAAAAAAGCCTTTGTAGTCGGTCTAATAACTTCATTTATTTTCCTTTGTCCTC‑TTTTTTTATTTCTGTGa                                                             >  1:862933/1‑80 (MQ=255)       aTTGGTGTAAAAAAAGCCTTTGTAGTCGGTCTAATAACTTCATTTATTTTCCTTTGTCCTC‑TTTTTTTATTTCTGTGaaa                                                           >  2:862935/1‑80 (MQ=255)       aTTGGTGTAAAAAAAGCCTTTGTAGTCGGTCTAATAACTGCATTTATTTTCCTTTGTCCTC‑TTTTTTTATTTCTGTGaaa                                                           >  2:862934/1‑80 (MQ=255)                  aaaaGCCTTTGTAGTCGGTCTAATAACTTCATTTATTTTCCTTTGTCCTC‑TTTTTTTATTTCTGTGAAAAActcttctct                                                >  1:862949/1‑80 (MQ=255)                   aaaGCCTTTGTAGTCGGTCTAATAACTTCATTTATTTTCCTTTGTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTta                                              >  1:862943/1‑81 (MQ=255)                          ttGTAGTCGGTCTAATAACTTCATTTATTTTCCTTTGTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGAt                                        >  1:862950/1‑80 (MQ=255)                              aGTCGGTCTAATAACTTCATTTATTTTCCTTTGTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGt                                   >  2:862953/1‑81 (MQ=255)                                  ggTCTAATAACTTCATTTATTTTCCTTTGTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGtctc                                <  2:862917/80‑1 (MQ=255)                                   gTCTAATAACTTCATTTATTTTCCTTTGTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCtt                              <  1:862919/81‑1 (MQ=255)                                       aataaCTTCATTTATTTTCCTTTGTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCTTTGGa                          >  2:862952/1‑81 (MQ=255)                                        ataaCTTCATTTATTTTCCTTTGTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCTTTGGAt                         >  2:862938/1‑81 (MQ=255)                                                 tttatttTCCTTTGTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCTTTGGATTAATGTTTa                <  2:862919/81‑1 (MQ=255)                                                 tttatttTCCTTTGTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCTTTGGACTAATGTTTa                >  2:862945/1‑81 (MQ=255)                                                        tCCTTTGTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCTTTGGATTAATGTTTACAAAtt          >  2:862961/1‑80 (MQ=255)                                                          cTTTGTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCTTTGGATTAATGTTTACAAATTTAg       >  1:862951/1‑81 (MQ=255)                                                             tGTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCTTTGGATTAATGTTTACAAATTTAgg      <  2:862920/79‑1 (MQ=255)                                                              gTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCTTTGGATTAATGTTTACAAATTTAGGTAt   >  1:862952/1‑81 (MQ=255)                                                              gTCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCTTTGGATTAATGTTTACAAATTTAGGTAt   >  2:862957/1‑81 (MQ=255)                                                               tCCTC‑TTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCTTTGGATTAATGTTTACAAATTTAGGTAtt  <  2:862921/81‑1 (MQ=255)                                                                    |                                                                             GATAAGATTGGTGTAAAAAAAGCCTTTGTAGTCGGTCTAATAACTTCATTTATTTTCCTTTGTCCTCTTTTTTTTATTTCTGTGAAAAACTCTTCTCTTATAGGATTATGTCTCTTTGGATTAATGTTTACAAATTTAGGTATT  >  NC_000913/4504883‑4505026

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