Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,505,167	C→A	G370G (GGC→GGA)	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,167	0	C	A	100.0%	118.0 / NA	34	G370G (GGC→GGA)	yjhB	putative sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (22/12);  total (22/12)

TTGGTTCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGC  >  NC_000913/4505033‑4505301                                                                                                                                       |                                                                                                                                       ttggttCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAgg                                                                                                                                   >  2:81555/1‑140 (MQ=255) ttggttCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAgg                                                                                                                                   <  1:91053/140‑1 (MQ=255)    gttCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGgtgt                                                                                                                                >  2:180139/1‑140 (MQ=255)                       aTTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACg                                                                                                             >  2:128197/1‑140 (MQ=255)                        ttACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACgg                                                                                                            <  1:239777/140‑1 (MQ=255)                                   caaAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTctct                                                                                                 >  2:126524/1‑140 (MQ=255)                                                aTTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAAttt                                                                                    <  2:515919/140‑1 (MQ=255)                                                      gACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAg                                                                              <  1:126524/140‑1 (MQ=255)                                                       aCCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGt                                                                             <  1:488020/140‑1 (MQ=255)                                                         cGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTg                                                                           >  2:329144/1‑140 (MQ=255)                                                          ggTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTgg                                                                          >  1:559452/1‑140 (MQ=255)                                                           gTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGt                                                                         >  2:6276/1‑140 (MQ=255)                                                                           aGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTa                                                         >  2:254568/1‑140 (MQ=255)                                                                           aGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTa                                                         >  2:481958/1‑140 (MQ=255)                                                                            ggggCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTaa                                                        <  2:544111/140‑1 (MQ=255)                                                                                   cTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTaa                                                        >  1:218017/1‑133 (MQ=255)                                                                                   cTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTaa                                                        <  2:218017/133‑1 (MQ=255)                                                                                                    ccTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTaa                                <  1:6276/140‑1 (MQ=255)                                                                                                             gCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGa                       <  1:254568/140‑1 (MQ=255)                                                                                                              cTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAg                      >  2:259/1‑140 (MQ=255)                                                                                                                  atacatTTCAGGATATTATGGATTGGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGc                                   >  2:19768/1‑123 (MQ=255)                                                                                                                  atacatTTCAGGATATTATGGATTGGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGc                                   <  1:19768/123‑1 (MQ=255)                                                                                                                  atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                  >  2:284366/1‑140 (MQ=255)                                                                                                                  atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                  >  2:46535/1‑140 (MQ=255)                                                                                                                  atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                  >  1:425883/1‑140 (MQ=255)                                                                                                                  atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                  >  2:532545/1‑140 (MQ=255)                                                                                                                   tacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTAGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCttt                 <  1:81555/140‑1 (MQ=255)                                                                                                                        ttCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGat            <  1:128197/140‑1 (MQ=255)                                                                                                                           aGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGattaa         >  1:390088/1‑140 (MQ=255)                                                                                                                            ggATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTAATAAACTATCCGTGGCTAAATGATTGGAGGCTTAATGattaat        >  2:103136/1‑140 (MQ=255)                                                                                                                                  tatGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGc  >  2:388170/1‑140 (MQ=255)                                                                                                                                  tatGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGc  >  2:444934/1‑140 (MQ=255)                                                                                                                                  tatGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGc  >  2:455452/1‑140 (MQ=255)                                                                                                                                  tatGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGc  >  1:400753/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                       TTGGTTCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGC  >  NC_000913/4505033‑4505301

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