Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,842,573	C→T	*157* (TAG→TAA)	hycI ←	protease involved in processing C‑terminal end of HycE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,842,573	0	C	T	100.0%	187.5 / NA	59	*157* (TAG→TAA)	hycI	protease involved in processing C‑terminal end of HycE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/28);  total (31/28)

AACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAA  >  NC_000913/2842438‑2842706                                                                                                                                        |                                                                                                                                      aaCGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAattac                                                                                                                                    <  1:438009/139‑1 (MQ=255) aaCGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAattac                                                                                                                                    >  2:1300635/1‑139 (MQ=255)                   cATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAAc                                                                                                                  >  1:938372/1‑139 (MQ=255)                                    gAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCAtt                                                                                                 <  2:1249696/139‑1 (MQ=255)                                     aaGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCt                                                                                                                     >  2:779835/1‑118 (MQ=255)                                     aaGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCt                                                                                                                     <  1:779835/118‑1 (MQ=255)                                         cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAGCTGCGCGAAGCCGCCATttcctt                                                                                            >  1:630953/1‑139 (MQ=255)                                         cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATttcctt                                                                                            >  2:764913/1‑139 (MQ=255)                                         cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATttcctt                                                                                            >  2:1059645/1‑139 (MQ=255)                                           cGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTcc                                                                                                                           >  1:638457/1‑106 (MQ=255)                                           cGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTcc                                                                                                                           <  2:638457/106‑1 (MQ=255)                                           cGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAGCTGCGCGAAGCCGCCATttccttcc                                                                                          <  2:630953/139‑1 (MQ=255)                                           cGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATttccttcc                                                                                          <  1:764913/139‑1 (MQ=255)                                             aTAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCa                                                                                        >  2:745436/1‑139 (MQ=255)                                                 gcgcagcgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCt                                                                                    <  1:628191/139‑1 (MQ=255)                                                  cgcagcgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCtt                                                                                   >  1:208755/1‑139 (MQ=255)                                                   gcagcgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTc                                                                                  <  1:306772/139‑1 (MQ=255)                                                    cagcgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTcc                                                                                 >  2:792242/1‑139 (MQ=255)                                                       cgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGt                                                                              >  1:167789/1‑139 (MQ=255)                                                       cgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGt                                                                              <  2:208755/139‑1 (MQ=255)                                                              ggCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGAt                                                                       <  1:959703/139‑1 (MQ=255)                                                                             ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATttcct                                                                                             >  1:303615/1‑102 (MQ=255)                                                                             ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATttcct                                                                                             <  2:303615/102‑1 (MQ=255)                                                                             ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                        >  1:981703/1‑139 (MQ=255)                                                                                      gCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTaa                                               <  1:8965/139‑1 (MQ=255)                                                                                        tGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATc                                             >  1:691600/1‑139 (MQ=255)                                                                                          gTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTAcc                                                         >  2:1008081/1‑125 (MQ=255)                                                                                          gTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTAcc                                                         <  1:1008081/125‑1 (MQ=255)                                                                                            tCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCt                                         <  2:691600/139‑1 (MQ=255)                                                                                               ccccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATaa                                                                     >  1:1238929/1‑108 (MQ=255)                                                                                               ccccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATaa                                                                     <  2:1238929/108‑1 (MQ=255)                                                                                                              cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATttccttcc                                                                                          >  1:52699/1‑72 (MQ=255)                                                                                                              cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATttccttcc                                                                                          <  2:52699/72‑1 (MQ=255)                                                                                                              cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAAc                                                                   >  1:1040082/1‑95 (MQ=255)                                                                                                              cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAAc                                                                   <  2:1040082/95‑1 (MQ=255)                                                                                                              cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTaaa                       >  1:291954/1‑139 (MQ=255)                                                                                                              cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTaaa                       >  2:215153/1‑139 (MQ=255)                                                                                                              cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTaaa                       >  1:142215/1‑139 (MQ=255)                                                                                                                 tGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGc                    <  1:1059645/139‑1 (MQ=255)                                                                                                                  gTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGcc                   <  1:1300635/139‑1 (MQ=255)                                                                                                                   tCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAgccg                                                                                                      <  1:1121034/55‑1 (MQ=255)                                                                                                                   tCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAgccg                                                                                                      >  2:1121034/1‑55 (MQ=255)                                                                                                                      tGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCACAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCcgacg               <  1:792242/139‑1 (MQ=255)                                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCAt                                                                                                  >  1:675902/1‑54 (MQ=255)                                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCAt                                                                                                  <  2:675902/54‑1 (MQ=255)                                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCtt                                                                                   >  2:939273/1‑69 (MQ=255)                                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCtt                                                                                   <  1:939273/69‑1 (MQ=255)                                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGtt                                                               <  2:1006525/89‑1 (MQ=255)                                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGtt                                                               >  1:1006525/1‑89 (MQ=255)                                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                        >  2:240704/1‑96 (MQ=255)                                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                        <  1:240704/96‑1 (MQ=255)                                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGc                                          <  2:751223/110‑1 (MQ=255)                                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGc                                          >  1:751223/1‑110 (MQ=255)                                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCcga                 <  2:653734/135‑1 (MQ=255)                                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCAATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCcga                 >  1:653734/1‑135 (MQ=255)                                                                                                                           ccTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATc          <  1:215153/139‑1 (MQ=255)                                                                                                                             ttAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGgtagta                         <  2:708454/122‑1 (MQ=255)                                                                                                                             ttAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGgtagta                         >  1:708454/1‑122 (MQ=255)                                                                                                                             ttAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCg        >  1:412036/1‑139 (MQ=255)                                                                                                                             ttAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCg        >  1:1108822/1‑139 (MQ=255)                                                                                                                                   aaaaaTTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGaa  >  2:581978/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                      AACGCCTTATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAA  >  NC_000913/2842438‑2842706

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