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

AGAGCATTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTG  >  NC_000913/1903350‑1903640                                                                                                                                                    |                                                                                                                                                agagCATTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCa                                                                                                                                              <  1:131890/151‑1 (MQ=255)             cAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTctct                                                                                                                                  >  2:362663/1‑151 (MQ=255)             cAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTctct                                                                                                                                  <  1:379788/151‑1 (MQ=255)                      cACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgc                                                                                                                          >  1:542648/1‑150 (MQ=255)                          ggCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTgcaagc                                                                                                                      <  1:111103/150‑1 (MQ=255)                                cccTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCgt                                                                                                               <  1:290702/151‑1 (MQ=255)                                 ccTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGGTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCgt                                                                                                               >  1:324137/1‑150 (MQ=255)                                      ggCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCt                                                                                                         >  2:509734/1‑151 (MQ=255)                                        ccGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTAt                                                                                                       >  1:359561/1‑151 (MQ=255)                                                  ccAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTc                                                                                             <  2:71975/151‑1 (MQ=255)                                                  ccAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTc                                                                                             >  2:44795/1‑151 (MQ=255)                                                   cAGGTACTGACCATCCTCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGCTTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTc                                                                                             >  1:440360/1‑150 (MQ=255)                                                             ccATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGct                                                                                  >  1:480140/1‑151 (MQ=255)                                                             ccATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGct                                                                                  >  2:152174/1‑151 (MQ=255)                                                                          tACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAcc                                                                     >  1:90392/1‑151 (MQ=255)                                                                                            ccAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCCATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATg                                                   >  2:515254/1‑151 (MQ=255)                                                                                                cACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGa                                               <  2:17818/151‑1 (MQ=255)                                                                                                   gCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAAt                                             >  1:386725/1‑150 (MQ=255)                                                                                                       cGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGAt                                        >  2:494401/1‑151 (MQ=255)                                                                                                        ggATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGAtt                                       <  1:28803/151‑1 (MQ=255)                                                                                                            aaGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCgg                                   <  1:267539/151‑1 (MQ=255)                                                                                                                            gCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGGGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATgg                   <  1:255723/151‑1 (MQ=255)                                                                                                                            gCAACCTGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATgg                   >  2:370088/1‑151 (MQ=255)                                                                                                                                  tGACAGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTAGCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGa              >  1:273198/1‑150 (MQ=255)                                                                                                                                      aGCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATc         <  2:57527/151‑1 (MQ=255)                                                                                                                                       gCGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCg        <  2:542791/151‑1 (MQ=255)                                                                                                                                        cGATTTCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGc       >  1:412020/1‑151 (MQ=255)                                                                                                                                           tttCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCgg    >  1:39292/1‑151 (MQ=255)                                                                                                                                            ttCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCgg    <  2:142480/150‑1 (MQ=255)                                                                                                                                             tCCTGGACCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTg  >  2:480794/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                AGAGCATTGGTGCAGGTATCGCACTGGCAATCCCTCTGGCCGCTGCGGGCCAGGTACTGACCATCATCGTTCGTACTATTACCGTTGCTTTCCAGCACGCTGCGGATAAGGCTGCTGATAACGGCAACCTGACAGCGATTTCCTGGATCCACGTTTCTTCTCTGTTCCTGCAAGCAATGCGTGTGGCTATTCCGGCCGTCATCGTTGCGCTGTCTGTTGGTACCAGCGAAGTACAGAACATGCTGAATGCGATTCCGGAAGTGGTGACCAATGGTCTGAATATCGCCGGTG  >  NC_000913/1903350‑1903640

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