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%	201.7 / NA	58	G370G (GGC→GGA)	yjhB	putative sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (43/15);  total (43/15)

GGTTCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCG  >  NC_000913/4505035‑4505302                                                                                                                                     |                                                                                                                                        ggttCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAgggg                                                                                                                                  >  1:196435/1‑138 (MQ=255) ggttCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGgtg                                                                                                                                  <  2:701057/140‑1 (MQ=255) ggttCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGgtg                                                                                                                                  <  1:249227/140‑1 (MQ=255)  gttCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGgtgt                                                                                                                                 >  1:554280/1‑140 (MQ=255)  gttCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGgtgt                                                                                                                                 >  2:70758/1‑140 (MQ=255)       aaaaTTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCat                                                                                                                            >  2:111496/1‑140 (MQ=255)         aaTTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCatta                                                                                                                          >  2:585085/1‑140 (MQ=255)                   tGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTa                                                                                                                <  2:730486/140‑1 (MQ=255)                   tGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTa                                                                                                                <  1:475198/140‑1 (MQ=255)                     aTTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACg                                                                                                              >  1:477435/1‑140 (MQ=255)                               aacaaAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTCCATTATTCATTGTTACGGTTGCATtct                                                                                                    >  2:38984/1‑140 (MQ=255)                                acaaAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATtctc                                                                                                   >  2:334723/1‑140 (MQ=255)                                           agGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTaa                                                                                        >  1:74580/1‑140 (MQ=255)                                                     aCCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGt                                                                              <  1:37278/140‑1 (MQ=255)                                                       cGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTg                                                                            >  2:264673/1‑140 (MQ=255)                                                        ggTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAAtt                                                                                      >  2:709756/1‑129 (MQ=255)                                                        ggTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAAtt                                                                                      <  1:709756/129‑1 (MQ=255)                                                        ggTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTgg                                                                           >  1:87027/1‑140 (MQ=255)                                                        ggTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTgg                                                                           >  2:624457/1‑140 (MQ=255)                                                                  ataACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATAtt                                                                 <  2:296195/140‑1 (MQ=255)                                                                    aaCTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTcc                                                               >  2:685962/1‑140 (MQ=255)                                                                    aaCTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTcc                                                               >  1:112604/1‑140 (MQ=255)                                                                         aGGGGCAACTGGTGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTa                                                          >  1:177956/1‑140 (MQ=255)                                                                         aGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTa                                                          >  2:689796/1‑140 (MQ=255)                                                                         aGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTa                                                          >  1:603405/1‑140 (MQ=255)                                                                              cAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAAtt                                                     <  2:74580/140‑1 (MQ=255)                                                                                aCTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTta                                                   <  1:626446/140‑1 (MQ=255)                                                                                           ggCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCg                                        >  1:667923/1‑140 (MQ=255)                                                                                            gCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGt                                       >  2:16521/1‑140 (MQ=255)                                                                                            gCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGt                                       >  2:755655/1‑140 (MQ=255)                                                                                            gCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGt                                       >  2:224405/1‑140 (MQ=255)                                                                                            gCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGt                                       >  2:432598/1‑140 (MQ=255)                                                                                              cGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTgg                                     <  2:685282/140‑1 (MQ=255)                                                                                                   cTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTaaa                                <  1:111496/140‑1 (MQ=255)                                                                                                             tacatacaTTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAgg                      >  1:93314/1‑140 (MQ=255)                                                                                                             tacatacaTTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAgg                      >  2:574764/1‑140 (MQ=255)                                                                                                             tacatacaTTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAgg                      >  1:504480/1‑140 (MQ=255)                                                                                                               catacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCt                    >  2:511221/1‑140 (MQ=255)                                                                                                               catacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCt                    >  2:554537/1‑140 (MQ=255)                                                                                                                atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAATTGATTGGAGGCtt                   >  2:97927/1‑140 (MQ=255)                                                                                                                atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                   >  2:467074/1‑140 (MQ=255)                                                                                                                atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                   >  2:94929/1‑140 (MQ=255)                                                                                                                atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                   <  2:76981/140‑1 (MQ=255)                                                                                                                atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                   >  2:214300/1‑140 (MQ=255)                                                                                                                atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                   >  2:223208/1‑140 (MQ=255)                                                                                                                atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                   >  2:718902/1‑140 (MQ=255)                                                                                                                atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                   >  1:645554/1‑140 (MQ=255)                                                                                                                atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                   >  1:215058/1‑140 (MQ=255)                                                                                                                atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                   >  2:312328/1‑140 (MQ=255)                                                                                                                atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                   >  1:44393/1‑140 (MQ=255)                                                                                                                    atTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATg               >  1:211391/1‑140 (MQ=255)                                                                                                                      ttCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGat             >  2:7870/1‑140 (MQ=255)                                                                                                                           gATCTTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGattaatt        >  1:76784/1‑140 (MQ=255)                                                                                                                            atatTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGattaatta       <  1:633938/140‑1 (MQ=255)                                                                                                                             tattatGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTAt      <  2:44393/140‑1 (MQ=255)                                                                                                                              attatGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATg     >  2:20408/1‑140 (MQ=255)                                                                                                                              attatGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATg     >  2:494538/1‑140 (MQ=255)                                                                                                                                 atGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCg  <  2:620401/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        GGTTCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCG  >  NC_000913/4505035‑4505302

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