Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,903,497	T→C	I139T (ATC→ACC)	manY →	mannose‑specific enzyme IIC component of PTS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,903,497	0	T	C	100.0%	189.7 / NA	50	I139T (ATC→ACC)	manY	mannose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/28);  total (22/28)

TTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTG  >  NC_000913/1903356‑1903640                                                                                                                                              |                                                                                                                                                ttGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTttc                                                                                                                                        >  2:386635/1‑151 (MQ=255)  tGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTttc                                                                                                                                        <  1:306959/150‑1 (MQ=255)   ggTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTttctt                                                                                                                                      <  2:613967/151‑1 (MQ=255)   ggTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTttct                                                                                                                                       <  2:289480/150‑1 (MQ=255)    gTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGTCCCACGTttct                                                                                                                                       >  2:439324/1‑149 (MQ=255)           tATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGtt                                                                                                                               <  2:142326/150‑1 (MQ=255)              cGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCt                                                                                                                            <  1:165874/150‑1 (MQ=255)              cGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTg                                                                                                                           <  2:328828/151‑1 (MQ=255)              cGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTg                                                                                                                           <  2:44090/151‑1 (MQ=255)               gCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgc                                                                                                                          >  2:402495/1‑151 (MQ=255)               gCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgc                                                                                                                          <  1:560330/151‑1 (MQ=255)               gCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgc                                                                                                                          <  1:552591/151‑1 (MQ=255)                    ggCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgcaagca                                                                                                                     <  1:468895/151‑1 (MQ=255)                    ggCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgcaag                                                                                                                       <  2:135139/149‑1 (MQ=255)                     gCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgcaagca                                                                                                                     >  2:171624/1‑150 (MQ=255)                      cAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAAt                                                                                                                   >  1:516295/1‑151 (MQ=255)                             tctGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGtgtg                                                                                                            >  2:322219/1‑151 (MQ=255)                                 gCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGc                                                                                                          >  2:308656/1‑149 (MQ=255)                                    gCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTc                                                                                                     <  2:770782/151‑1 (MQ=255)                                       gCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTcc                                                                                                    <  2:779476/149‑1 (MQ=255)                                       gCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCgg                                                                                                  >  2:272537/1‑151 (MQ=255)                                          ggCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCg                                                                                               >  1:476948/1‑151 (MQ=255)                                               ggTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATc                                                                                          <  2:187680/151‑1 (MQ=255)                                                     gACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAGGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTgcg                                                                                    <  1:683421/151‑1 (MQ=255)                                                           catcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgtc                                                                               <  2:550889/150‑1 (MQ=255)                                                                ttcgtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTgg                                                                         >  1:133763/1‑151 (MQ=255)                                                                    tACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAcc                                                                     >  2:617783/1‑151 (MQ=255)                                                                     aCTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCa                                                                    <  1:7336/151‑1 (MQ=255)                                                                       tATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGc                                                                  >  1:575350/1‑151 (MQ=255)                                                                        aTTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCg                                                                 <  2:516234/151‑1 (MQ=255)                                                                         ttACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCg                                                                 >  1:649117/1‑150 (MQ=255)                                                                           aCCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGa                                                                >  2:444776/1‑149 (MQ=255)                                                                             cGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGt                                                             >  1:550779/1‑150 (MQ=255)                                                                             cGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGt                                                             >  2:368033/1‑150 (MQ=255)                                                                             cGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTa                                                            <  1:646669/151‑1 (MQ=255)                                                                               ttGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTAc                                                           >  1:431149/1‑150 (MQ=255)                                                                                            cGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAAt                                             >  1:152753/1‑151 (MQ=255)                                                                                                 cGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGAt                                        >  1:530304/1‑151 (MQ=255)                                                                                                       aGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGa                                  <  2:230820/151‑1 (MQ=255)                                                                                                               gATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggtg                           <  2:295627/150‑1 (MQ=255)                                                                                                                 tAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAcc                        >  2:402123/1‑151 (MQ=255)                                                                                                                     ggCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATg                    <  2:1477/151‑1 (MQ=255)                                                                                                                              aCAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAata           <  1:524329/151‑1 (MQ=255)                                                                                                                                aGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAatat          <  2:271522/150‑1 (MQ=255)                                                                                                                                aGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATc         <  1:587131/151‑1 (MQ=255)                                                                                                                                 gCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCg        <  2:116231/151‑1 (MQ=255)                                                                                                                                  cGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGc       >  1:263046/1‑151 (MQ=255)                                                                                                                                   gATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGcc      >  2:533112/1‑151 (MQ=255)                                                                                                                                       tCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGt   <  2:400947/150‑1 (MQ=255)                                                                                                                                        ccTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTg  <  2:557898/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                TTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTG  >  NC_000913/1903356‑1903640

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