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

GTTGGTTCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGT  >  NC_000913/4505032‑4505303                                                                                                                                        |                                                                                                                                         gttggttCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAg                                                                                                                                      <  1:445704/140‑1 (MQ=255)    ggttCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGgtg                                                                                                                                   <  1:182124/140‑1 (MQ=255)     gttCCAAAATTTATATATGATTACTTTCCAACAAATTTAAGAGGAGTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGgtgt                                                                                                                                  >  2:696542/1‑140 (MQ=255)     gttCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGgtgt                                                                                                                                  >  2:425/1‑140 (MQ=255)     gttCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGACTAGgtgt                                                                                                                                  >  1:191796/1‑140 (MQ=255)          aaaaTTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCat                                                                                                                             >  2:204142/1‑140 (MQ=255)                tatataTGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCa                                                                                                                       <  2:173565/140‑1 (MQ=255)                      tGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTa                                                                                                                 <  2:680884/140‑1 (MQ=255)                            cTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTg                                                                                                           >  1:53624/1‑140 (MQ=255)                            cTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGACTAGGTGTTTCATTATTCATTGTTACGGTTg                                                                                                           <  2:191796/140‑1 (MQ=255)                                  aacaaAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGCTGCATtct                                                                                                     >  2:188937/1‑140 (MQ=255)                                    caaAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTctct                                                                                                   <  1:65573/140‑1 (MQ=255)                                                 aTTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAAttt                                                                                      <  1:268946/140‑1 (MQ=255)                                                       gACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAg                                                                                <  2:486320/140‑1 (MQ=255)                                                       gACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAg                                                                                <  1:672282/140‑1 (MQ=255)                                                          cGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTg                                                                             >  2:178656/1‑140 (MQ=255)                                                           ggTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTgg                                                                            >  2:493178/1‑140 (MQ=255)                                                            gTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGt                                                                           >  1:397159/1‑140 (MQ=255)                                                            gTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGt                                                                           >  1:256972/1‑140 (MQ=255)                                                                       aaCTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTcc                                                                >  1:67200/1‑140 (MQ=255)                                                                       aaCTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTcc                                                                >  1:477619/1‑140 (MQ=255)                                                                            aGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTa                                                           >  2:593786/1‑140 (MQ=255)                                                                            aGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTa                                                           >  2:634750/1‑140 (MQ=255)                                                                                               gCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGt                                        >  2:608210/1‑140 (MQ=255)                                                                                                 cGCACCTGTATTAGCTACATTCATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTTCGGTTGCATTATCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTgg                                      <  2:587719/140‑1 (MQ=255)                                                                                                      cTGTATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTaaa                                 <  2:105538/140‑1 (MQ=255)                                                                                                         tATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGa                              >  2:261128/1‑140 (MQ=255)                                                                                                         tATTAGCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGa                              >  2:54753/1‑140 (MQ=255)                                                                                                              gCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGa                         <  2:442065/140‑1 (MQ=255)                                                                                                              gCTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTCGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGa                         >  2:720446/1‑140 (MQ=255)                                                                                                               cTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAg                        <  1:696542/140‑1 (MQ=255)                                                                                                               cTACATACATTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATT‑CAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAgg                       >  2:679826/1‑140 (MQ=255)                                                                                                                tacatacaTTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAgg                       >  2:378301/1‑140 (MQ=255)                                                                                                                tacatacaTTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAgg                       >  2:697386/1‑140 (MQ=255)                                                                                                                tacatacaTTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAgg                       >  1:523201/1‑140 (MQ=255)                                                                                                                  catacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCt                     <  1:250831/140‑1 (MQ=255)                                                                                                                  catacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATCCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCt                     <  2:15614/140‑1 (MQ=255)                                                                                                                   atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                    >  2:627978/1‑140 (MQ=255)                                                                                                                   atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                    >  2:617410/1‑140 (MQ=255)                                                                                                                   atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                    >  1:206754/1‑140 (MQ=255)                                                                                                                   atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                    >  2:746981/1‑140 (MQ=255)                                                                                                                   atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                    >  2:265764/1‑140 (MQ=255)                                                                                                                   atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                    <  2:232860/140‑1 (MQ=255)                                                                                                                   atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                    >  2:195344/1‑140 (MQ=255)                                                                                                                   atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                    >  1:170726/1‑140 (MQ=255)                                                                                                                   atacatTTCAGGATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCtt                    >  1:446523/1‑140 (MQ=255)                                                                                                                              gATATTATGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGattaatt         >  2:136612/1‑140 (MQ=255)                                                                                                                                 attatGGATTAGGTGTTTCATTATTCATTGTTACGGTT‑CATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATgg     >  2:140277/1‑140 (MQ=255)                                                                                                                                  ttatGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATgg     <  1:265764/140‑1 (MQ=255)                                                                                                                                   tatGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGc    >  2:472479/1‑140 (MQ=255)                                                                                                                                   tatGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGc    >  1:353977/1‑140 (MQ=255)                                                                                                                                   tatGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGc    >  1:416457/1‑140 (MQ=255)                                                                                                                                    atGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCg   <  2:53624/140‑1 (MQ=255)                                                                                                                                     tGGATTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCgt  <  1:54753/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         GTTGGTTCCAAAATTTATATATGATTACTTTCCAACAAAATTAAGAGGATTAGGGACCGGTCTTATTTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGT  >  NC_000913/4505032‑4505303

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