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%	181.4 / NA	57	*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/26);  total (31/26)

TACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGA  >  NC_000913/2842454‑2842705                                                                                                                        |                                                                                                                                     tACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTa                                                                                                                   <  1:642822/139‑1 (MQ=255)   cATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAAc                                                                                                                 >  2:184602/1‑139 (MQ=255)    aTGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCttcttc                                                                                                                           <  2:307609/128‑1 (MQ=255)    aTGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCttcttc                                                                                                                           >  1:307609/1‑128 (MQ=255)            cGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAg                                                                                                        >  1:21923/1‑139 (MQ=255)             gAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAgc                                                                                                       <  2:769977/139‑1 (MQ=255)                         cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                            >  2:775041/1‑106 (MQ=255)                         cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                            <  1:775041/106‑1 (MQ=255)                             aTAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCa                                                                                       <  2:205639/139‑1 (MQ=255)                             aTAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCa                                                                                       <  2:893601/139‑1 (MQ=255)                             aTAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCa                                                                                       <  1:93847/139‑1 (MQ=255)                                  cgcagcgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                            >  2:708085/1‑97 (MQ=255)                                  cgcagcgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                            <  1:708085/97‑1 (MQ=255)                                  cgcagcgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCAtt                                                                                                >  2:1398324/1‑125 (MQ=255)                                  cgcagcgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCAtt                                                                                                <  1:1398324/125‑1 (MQ=255)                                       cgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGt                                                                             <  1:184602/139‑1 (MQ=255)                                                     tAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGt                                                               <  1:785283/139‑1 (MQ=255)                                                          tttGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATaa                                                                    <  2:1382107/129‑1 (MQ=255)                                                          tttGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATaa                                                                    >  1:1382107/1‑129 (MQ=255)                                                                  cAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATaa                                                                    <  2:686113/121‑1 (MQ=255)                                                                  cAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATaa                                                                    >  1:686113/1‑121 (MQ=255)                                                                               ccccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCttcttc                                                                                                                           >  1:818198/1‑53 (MQ=255)                                                                               ccccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCttcttc                                                                                                                           <  2:818198/53‑1 (MQ=255)                                                                               ccccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTggg                                     >  2:567011/1‑139 (MQ=255)                                                                                 ccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAgccgcc                                                                                                   <  2:1259769/75‑1 (MQ=255)                                                                                 ccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAgccgcc                                                                                                   >  1:1259769/1‑75 (MQ=255)                                                                                  cATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                  <  2:1415612/139‑1 (MQ=255)                                                                                    tCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATc                                >  2:714700/1‑139 (MQ=255)                                                                                    tCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATc                                >  1:1223901/1‑139 (MQ=255)                                                                                    tCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATc                                >  2:452137/1‑139 (MQ=255)                                                                                          aCATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCTCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGgtag                          >  2:637066/1‑139 (MQ=255)                                                                                            aTCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                  <  2:620665/129‑1 (MQ=255)                                                                                            aTCCCTGTCCTGATTCCTTAATGAAAAATTACCCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                  >  1:620665/1‑129 (MQ=255)                                                                                              cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATttccttcc                                                                                         <  2:425355/72‑1 (MQ=255)                                                                                              cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATttccttcc                                                                                         >  1:425355/1‑72 (MQ=255)                                                                                              cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTaaa                      >  1:876470/1‑139 (MQ=255)                                                                                                   tCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCcg                 >  2:859336/1‑139 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCATCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTc                                                                                 <  1:2020/70‑1 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCAtt                                                                                                <  2:578630/55‑1 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCAtt                                                                                                >  1:578630/1‑55 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTc                                                                                 >  2:2020/1‑70 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATaa                                                                    <  1:1104838/83‑1 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATaa                                                                    >  2:1104838/1‑83 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                       <  2:1449345/96‑1 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                       >  1:1449345/1‑96 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCt                                        >  2:182992/1‑111 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCt                                        <  1:182992/111‑1 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat            >  2:974232/1‑139 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat            >  2:1258810/1‑139 (MQ=255)                                                                                                        aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat            >  1:1127176/1‑139 (MQ=255)                                                                                                             ttAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATttccttcc                                                                                         <  2:1381824/57‑1 (MQ=255)                                                                                                             ttAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATttccttcc                                                                                         >  1:1381824/1‑57 (MQ=255)                                                                                                               aaTGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGa                                                                                                          <  2:182829/38‑1 (MQ=255)                                                                                                               aaTGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGa                                                                                                          >  1:182829/1‑38 (MQ=255)                                                                                                                 tGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTg   <  2:1223901/139‑1 (MQ=255)                                                                                                                  gAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGa  >  2:796433/1‑139 (MQ=255)                                                                                                                  gAAAAATAACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGa  >  1:964067/1‑139 (MQ=255)                                                                                                                        |                                                                                                                                     TACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGA  >  NC_000913/2842454‑2842705

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: ‑