Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,237,239	A→C	I115L (ATC→CTC)	yjbG →	extracellular polysaccharide export OMA protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,237,239	0	A	C	100.0%	188.9 / NA	54	I115L (ATC→CTC)	yjbG	extracellular polysaccharide export OMA protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (35/19);  total (35/19)

CAACGGCGGCAGCATTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAAATCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGCCGTCCTGGCAAGCAGCCATTCACTC  >  NC_000913/4237096‑4237383                                                                                                                                                |                                                                                                                                                 cAACGGCGGCAGCATTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCt                                                                                                                                           >  1:715662/1‑151 (MQ=255) cAACGGCGGCAGCATTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCt                                                                                                                                           <  2:259897/151‑1 (MQ=255)  aaCGGCGGCAGCATTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCtt                                                                                                                                          >  2:665033/1‑151 (MQ=255)  aaCGGCGGCAGCATTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCtt                                                                                                                                          <  1:705739/151‑1 (MQ=255)  aaCGGCGGCAGCATTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCtt                                                                                                                                          <  2:275861/151‑1 (MQ=255)    cggcggCAGCATTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTg                                                                                                                                         >  1:763233/1‑150 (MQ=255)     ggcggcAGCATTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGAt                                                                                                                                       >  1:186508/1‑151 (MQ=255)      gcggcAGCATTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATc                                                                                                                                      >  2:482895/1‑151 (MQ=255)         gcagcaTTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCg                                                                                                                                   >  2:696816/1‑151 (MQ=255)            gcaTTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGata                                                                                                                                <  2:12008/151‑1 (MQ=255)             caTTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGatat                                                                                                                               <  1:321431/151‑1 (MQ=255)             caTTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGatat                                                                                                                               >  1:870410/1‑151 (MQ=255)             caTTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGatat                                                                                                                               <  1:131073/151‑1 (MQ=255)              aTTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGTGCAGGCAAAAACTCAATCTTGATCCCGata                                                                                                                                >  2:617115/1‑149 (MQ=255)               ttGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCg                                                                                                                             >  2:857852/1‑151 (MQ=255)                 gCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGt                                                                                                                            <  1:573493/150‑1 (MQ=255)                  cGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTcc                                                                                                                          >  2:60772/1‑151 (MQ=255)                     cagcagcagGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCcgcg                                                                                                                       >  2:523451/1‑151 (MQ=255)                        cagcagGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTg                                                                                                                    <  1:876526/151‑1 (MQ=255)                         agcagGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGc                                                                                                                   <  1:769031/151‑1 (MQ=255)                         agcagGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGc                                                                                                                   <  1:276555/151‑1 (MQ=255)                            agGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCg                                                                                                                 >  2:889977/1‑150 (MQ=255)                                  tGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAAcgcg                                                                                                           <  2:423944/150‑1 (MQ=255)                                      gACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAAc                                                                                                      >  1:562519/1‑151 (MQ=255)                                       acacGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAAcc                                                                                                     >  2:860255/1‑151 (MQ=255)                                         acGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACccg                                                                                                   >  1:483044/1‑151 (MQ=255)                                          cGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACccgc                                                                                                  >  1:690011/1‑151 (MQ=255)                                              tGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGtt                                                                                              >  2:263798/1‑151 (MQ=255)                                                  aGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGc                                                                                            >  1:906908/1‑149 (MQ=255)                                                    aaCGGGGCGCAGATTCCAGCGCGGACGATGCCGCTGCGATTACCGCTTTACCCCAACAAATTCAGGCGTTGAAAGTGAGGGGCGGGAAAAACCTCATTCTTGATCCCGATATCGTCCGCGTTGCGGAACGCGGTACCCCGCCGTTGc                                                                                            <  2:890458/147‑1 (MQ=255)                                                           cgcAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTac                                                                                 >  2:813235/1‑151 (MQ=255)                                                            gcAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTac                                                                                 >  1:742454/1‑150 (MQ=255)                                                                 ttCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCtg                                                                           >  2:749069/1‑151 (MQ=255)                                                                       cgcCGACGAGGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGGGACGGGCAGGCAAAAACTCAATCTTGATCCCGATACCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGt                                                                      <  1:432502/150‑1 (MQ=255)                                                                          cgacgaTGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCg                                                                    <  1:766594/149‑1 (MQ=255)                                                                                    gCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGt                                                        >  1:890054/1‑151 (MQ=255)                                                                                    gCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACAAccgccg                                                         >  2:293338/1‑150 (MQ=255)                                                                                     cTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTc                                                       >  1:866153/1‑151 (MQ=255)                                                                                       gcgttaaaCGCCTCCCGCCAGGAAATCCAGGCTTTAAAAGTGGCGGGAAGGAAAAACCTCACTCTTGATCCCAATATGTCCCGGGTGGCCAATCCGGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTc                                                       <  1:161888/143‑1 (MQ=255)                                                                                       gCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTcc                                                      >  2:691052/1‑150 (MQ=255)                                                                                          aTTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACgg                                                  >  1:368955/1‑151 (MQ=255)                                                                                          aTTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACgg                                                  >  2:776814/1‑151 (MQ=255)                                                                                           ttACCGCCTCACCCCAGAAAATTCGGGCGTTAAAGGGGCCGGGCAGGAAAAAACTCAATCTTGATCCAGATATCGCCCGCGTTGCCAAACGCGGTACCCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACg                                                   <  2:69236/149‑1 (MQ=255)                                                                                              aCGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCAc                                              >  1:48373/1‑151 (MQ=255)                                                                                               cGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACg                                             >  1:85642/1‑151 (MQ=255)                                                                                                   ttccGCCAGCAAATTCGGGCGTTGAAGGGGCCGGGCAGCAAAAACCTCAATCTTGACCCCGATATCGCCCGCGTTGCCGAACGCGGTACCCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACgt                                            <  2:870391/145‑1 (MQ=255)                                                                                                    tACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACgttgtt                                        >  1:638169/1‑151 (MQ=255)                                                                                                        ccAGCAAATTCGGGCTTTAAAAGTGACGGGCAGGAAAAAACTCAATCTTCACCCCGATATCGTCCGCTTTGCCAAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCg                                      <  2:553750/149‑1 (MQ=255)                                                                                                          aGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCt                                  >  1:400634/1‑151 (MQ=255)                                                                                                            cAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCtt                                 >  1:632637/1‑150 (MQ=255)                                                                                                             aaaTTCAGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTAt                               >  1:861131/1‑151 (MQ=255)                                                                                                                   aGGCGTTGAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGCCg                         >  2:313075/1‑151 (MQ=255)                                                                                                                          gAAAGTGACGGGCAGGCAAAAACTCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGCCGTCCTg                    >  2:485642/1‑149 (MQ=255)                                                                                                                           aaaGTGACGGGCAGGAAAAAACTCAATCTTGACCCCGATATCGTCCGCGTTGCCGACCGGGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGCCGTCCTgg                   <  1:216605/149‑1 (MQ=255)                                                                                                                                           aaaaacTCAATCTTGACCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGCCGTCCTGGCAAGCAGCCATTCACTc  <  1:468841/150‑1 (MQ=255)                                                                                                                                           aaaaacTCAATCTTGACCCCGATATCGCCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGCCGTCCTGGCAAGCAGCCATTCACt   <  1:56736/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 CAACGGCGGCAGCATTGCGTCAGCAGCAGGCGTTGCTGACACGGCTGGCAGAACAGGGCGCAGATTCCAGCGCCGACGATGCCGCTGCGATTAACGCCTTACGCCAGCAAATTCAGGCGTTGAAAGTGACGGGCAGGCAAAAAATCAATCTTGATCCCGATATCGTCCGCGTTGCCGAACGCGGTAACCCGCCGTTGCAGGGCAACTACACGCTGTGGGTCGGACCACCGCCGTCCACGGTCACGTTGTTCGGGCTTATCAGCCGTCCTGGCAAGCAGCCATTCACTC  >  NC_000913/4237096‑4237383

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