Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,842,591	T→C	L151L (TTA→TTG)	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,591	0	T	C	100.0%	236.2 / NA	63	L151L (TTA→TTG)	hycI	protease involved in processing C‑terminal end of HycE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/27);  total (36/27)

ACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCAC  >  NC_000913/2842455‑2842722                                                                                                                                         |                                                                                                                                     aCATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCaa                                                                                                                                    >  1:390985/1‑139 (MQ=255)       cGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTgcgc                                                                                                                              >  2:472923/1‑139 (MQ=255)       cGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTgcgc                                                                                                                              >  2:1231805/1‑139 (MQ=255)        gCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGcgcg                                                                                                                             >  1:1008494/1‑139 (MQ=255)        gCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGcgcg                                                                                                                             <  2:301877/139‑1 (MQ=255)                        cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATtt‑cctt                                                                                                            >  2:1423621/1‑139 (MQ=255)                          cGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATtt‑ccttcc                                                                                                          <  2:1154079/139‑1 (MQ=255)                            aTAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATttcccttccc                                                                                                         >  2:1358399/1‑139 (MQ=255)                                      cgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGt                                                                                              >  2:448156/1‑139 (MQ=255)                                               cAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTTATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATaa                                                                                     <  1:1215799/139‑1 (MQ=255)                                               cAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATaa                                                                                     <  2:800027/139‑1 (MQ=255)                                               cAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATaa                                                                                     <  2:1430510/139‑1 (MQ=255)                                                       cATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTc                                                                             >  1:327043/1‑139 (MQ=255)                                                        aTTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCt                                                                            >  1:1249872/1‑139 (MQ=255)                                                            ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                                        >  1:190300/1‑139 (MQ=255)                                                            ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                                        <  1:472923/139‑1 (MQ=255)                                                            ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                                        >  1:492465/1‑139 (MQ=255)                                                            ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                                        >  2:283863/1‑139 (MQ=255)                                                              cATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCa                                                                      >  2:200438/1‑139 (MQ=255)                                                               aTCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCAt                                                                     <  2:190300/139‑1 (MQ=255)                                                               aTCAGAGCTGGTTCACCCCATCAAGAACATCCCCGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCAt                                                                     <  2:1464401/139‑1 (MQ=255)                                                                              ccccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTa                                                                <  2:1132476/129‑1 (MQ=255)                                                                              ccccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTa                                                                >  1:1132476/1‑129 (MQ=255)                                                                                 cATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                                   <  1:200438/139‑1 (MQ=255)                                                                                 cATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                                   <  2:492465/139‑1 (MQ=255)                                                                                 cATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                                   <  1:1231805/139‑1 (MQ=255)                                                                                   tCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTTCTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                                   <  1:26729/137‑1 (MQ=255)                                                                                   tCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTTCTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                                   >  2:26729/1‑137 (MQ=255)                                                                                   tCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAgccgcc                                                                                                                     <  1:650382/72‑1 (MQ=255)                                                                                   tCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAgccgcc                                                                                                                     >  2:650382/1‑72 (MQ=255)                                                                                   tCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATc                                                 >  2:963804/1‑139 (MQ=255)                                                                                            tCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTaa                                        >  2:328773/1‑139 (MQ=255)                                                                                             cccTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATtt‑ccttcc                                                                                                          >  1:830779/1‑72 (MQ=255)                                                                                             cccTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATtt‑ccttcc                                                                                                          <  2:830779/72‑1 (MQ=255)                                                                                             cccTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                                   >  2:416433/1‑127 (MQ=255)                                                                                             cccTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                                   <  1:416433/127‑1 (MQ=255)                                                                                             cccTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTaaa                                       >  2:829779/1‑139 (MQ=255)                                                                                               cTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAg                                     >  1:1012195/1‑139 (MQ=255)                                                                                               cTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAg                                     <  2:447945/139‑1 (MQ=255)                                                                                                       aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACg                                                                                  >  2:221122/1‑86 (MQ=255)                                                                                                       aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACg                                                                                  <  1:221122/86‑1 (MQ=255)                                                                                                       aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCAt                                                                     >  1:351146/1‑99 (MQ=255)                                                                                                       aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCAt                                                                     <  2:351146/99‑1 (MQ=255)                                                                                                       aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat                             >  1:1443425/1‑139 (MQ=255)                                                                                                       aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat                             >  2:1355027/1‑139 (MQ=255)                                                                                                       aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat                             >  1:1243575/1‑139 (MQ=255)                                                                                                       aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat                             >  1:141629/1‑139 (MQ=255)                                                                                                        ttCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat                             <  1:1358399/139‑1 (MQ=255)                                                                                                        ttCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGata                            <  1:34210/139‑1 (MQ=255)                                                                                                            ttAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                                   <  1:921303/112‑1 (MQ=255)                                                                                                            ttAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                                   >  2:921303/1‑112 (MQ=255)                                                                                                            ttAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCg                        >  1:51625/1‑139 (MQ=255)                                                                                                            ttAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCg                        >  1:1133721/1‑139 (MQ=255)                                                                                                            ttAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCg                        >  1:744943/1‑139 (MQ=255)                                                                                                            ttAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCg                        >  2:269399/1‑139 (MQ=255)                                                                                                              aaTGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCgg                                               <  1:725110/114‑1 (MQ=255)                                                                                                              aaTGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCgg                                               >  2:142919/1‑114 (MQ=255)                                                                                                              aaTGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCgg                                               >  2:725110/1‑114 (MQ=255)                                                                                                              aaTGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCgg                                               <  1:142919/114‑1 (MQ=255)                                                                                                                       ctactCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCg             <  2:327043/139‑1 (MQ=255)                                                                                                                                ttcCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGgtagt                                          >  2:439408/1‑101 (MQ=255)                                                                                                                                ttcCACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGgtagt                                          <  1:439408/101‑1 (MQ=255)                                                                                                                                  ccACCGCCAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCAc  <  1:1355027/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                     ACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTT‑CCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCAC  >  NC_000913/2842455‑2842722

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