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%	212.8 / NA	55	I139T (ATC→ACC)	manY	mannose‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/31);  total (24/31)

AGCATTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCAT  >  NC_000913/1903352‑1903644                                                                                                                                                  |                                                                                                                                                    agCATTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACg                                                                                                                                                >  2:30537/1‑151 (MQ=255)    aTTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGttt                                                                                                                                             <  2:82417/151‑1 (MQ=255)      tGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTttct                                                                                                                                           >  2:477447/1‑151 (MQ=255)      tGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTttct                                                                                                                                           >  1:354674/1‑151 (MQ=255)                aTCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTcc                                                                                                                                 >  1:461630/1‑151 (MQ=255)                aTCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTc                                                                                                                                  <  2:48284/150‑1 (MQ=255)                  cGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCt                                                                                                                                <  2:375280/150‑1 (MQ=255)                  cGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTg                                                                                                                               >  1:442121/1‑151 (MQ=255)                   gCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgc                                                                                                                              >  2:321886/1‑151 (MQ=255)                      cTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgca                                                                                                                             >  2:665785/1‑149 (MQ=255)                        ggCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgcaagc                                                                                                                          <  1:318894/150‑1 (MQ=255)                           aaTCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATg                                                                                                                      <  2:722347/151‑1 (MQ=255)                            aTCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGc                                                                                                                     <  1:292430/151‑1 (MQ=255)                            aTCACTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGCTTCTTCTCTGTTCCTGCAAGCAATg                                                                                                                      >  2:445744/1‑150 (MQ=255)                              cccTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCg                                                                                                                    <  1:261408/150‑1 (MQ=255)                               ccTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCgtg                                                                                                                  <  1:296479/151‑1 (MQ=255)                                     gCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCt                                                                                                             >  1:639340/1‑150 (MQ=255)                                       cGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTAtt                                                                                                          <  2:231620/151‑1 (MQ=255)                                        gCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTc                                                                                                         >  1:28415/1‑151 (MQ=255)                                           gCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCgg                                                                                                      <  1:232972/151‑1 (MQ=255)                                               gCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTTACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCCTGTAGCTATTCCGAcca                                                                                                   <  2:232163/150‑2 (MQ=255)                                                 cAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCa                                                                                                <  1:516347/151‑1 (MQ=255)                                                 cAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCa                                                                                                <  2:696779/151‑1 (MQ=255)                                                         gACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTgcg                                                                                        <  1:399837/151‑1 (MQ=255)                                                         gACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTgcg                                                                                        >  2:749380/1‑151 (MQ=255)                                                         gACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTgcg                                                                                        <  1:305909/151‑1 (MQ=255)                                                           ccATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGct                                                                                      >  1:409279/1‑151 (MQ=255)                                                             atcatcGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGctgt                                                                                    >  1:366933/1‑151 (MQ=255)                                                                     tcgtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGt                                                                            >  1:569676/1‑151 (MQ=255)                                                                      cgtACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGt                                                                            <  2:204845/150‑1 (MQ=255)                                                                                ccGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGt                                                                 <  1:170310/151‑1 (MQ=255)                                                                                    tGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACa                                                              <  1:634419/150‑1 (MQ=255)                                                                                        ttCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGaa                                                           <  2:19577/149‑1 (MQ=255)                                                                                        ttCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAAc                                                          >  1:616602/1‑150 (MQ=255)                                                                                            aGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCt                                                     <  2:44339/151‑1 (MQ=255)                                                                                             gCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTg                                                    >  1:568029/1‑151 (MQ=255)                                                                                              cACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTg                                                    >  2:674817/1‑150 (MQ=255)                                                                                                cGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAAt                                                 <  1:118343/151‑1 (MQ=255)                                                                                                 gCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATg                                                >  2:710628/1‑151 (MQ=255)                                                                                                  cTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGc                                               >  2:495835/1‑151 (MQ=255)                                                                                                      ggATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGAtt                                           <  1:120272/151‑1 (MQ=255)                                                                                                           aGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCgg                                       >  2:640892/1‑150 (MQ=255)                                                                                                                gctgATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtgg                                 <  2:33662/151‑1 (MQ=255)                                                                                                                gctgATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgt                                   <  1:259232/149‑1 (MQ=255)                                                                                                                   gATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAgtggt                                <  2:613929/149‑1 (MQ=255)                                                                                                                     tAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGAcc                            >  1:365686/1‑151 (MQ=255)                                                                                                                             aCCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCt                    <  1:551361/151‑1 (MQ=255)                                                                                                                                 gACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGa                  <  1:581267/149‑1 (MQ=255)                                                                                                                                  aCAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAata               <  2:299112/151‑1 (MQ=255)                                                                                                                                    aGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATc             <  2:442080/151‑1 (MQ=255)                                                                                                                                       gATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGcc          >  2:359341/1‑151 (MQ=255)                                                                                                                                          ttCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCgg        <  2:810564/150‑1 (MQ=255)                                                                                                                                          ttCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGt       >  2:133395/1‑151 (MQ=255)                                                                                                                                              tGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGc    >  2:321593/1‑150 (MQ=255)                                                                                                                                               ggACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCAt  <  1:339260/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AGCATTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTGGCAT  >  NC_000913/1903352‑1903644

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