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%	133.7 / NA	42	*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 (25/17);  total (25/17)

TATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCG  >  NC_000913/2842445‑2842700                                                                                                                                 |                                                                                                                                tATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                       <  1:287158/139‑1 (MQ=255)       gCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCg                                                                                                                 >  1:432451/1‑139 (MQ=255)            cATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAAc                                                                                                            >  1:628039/1‑139 (MQ=255)             aTGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCttcttc                                                                                                                      <  2:523078/128‑1 (MQ=255)             aTGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCttcttc                                                                                                                      >  1:523078/1‑128 (MQ=255)                  gCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                       <  2:584915/122‑1 (MQ=255)                  gCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                       >  1:584915/1‑122 (MQ=255)                  gCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGcgcg                                                                                                      >  1:433335/1‑139 (MQ=255)                          cATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAgccgcc                                                                                              <  2:432451/139‑1 (MQ=255)                                  cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTcc                                                                                                                     >  2:328563/1‑108 (MQ=255)                                  cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTcc                                                                                                                     <  1:328563/108‑1 (MQ=255)                                    cGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAgccgcc                                                                                              >  2:133841/1‑129 (MQ=255)                                    cGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAgccgcc                                                                                              <  1:133841/129‑1 (MQ=255)                                            gcagcgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTc                                                                            <  2:34912/139‑1 (MQ=255)                                            gcagcgcaTCGGACAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTc                                                                            <  1:769823/139‑1 (MQ=255)                                                cgcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGt                                                                        <  1:208656/139‑1 (MQ=255)                                                 gcaTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTc                                                                       >  1:389526/1‑139 (MQ=255)                                                   aTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGt                                                                     >  1:237279/1‑139 (MQ=255)                                                          aaTTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATaaa                                                              <  2:433335/139‑1 (MQ=255)                                                           aTTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAAc                                                             <  2:389526/139‑1 (MQ=255)                                                                     tGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTAcc                                                   <  1:645247/139‑1 (MQ=255)                                                                      ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                       >  1:252672/1‑70 (MQ=255)                                                                      ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTcttctt                                                                                                                       <  2:252672/70‑1 (MQ=255)                                                                      ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                  >  2:112212/1‑139 (MQ=255)                                                                                           cATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                             <  1:732004/139‑1 (MQ=255)                                                                                             tCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATc                           >  2:596281/1‑139 (MQ=255)                                                                                               aaGAACATCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCgg                         <  1:68018/139‑1 (MQ=255)                                                                                                      tCCCTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTaa                  >  1:55837/1‑139 (MQ=255)                                                                                                       cccTGTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTaaa                 >  2:270012/1‑139 (MQ=255)                                                                                                           gTCCTGATTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGcc             >  1:83026/1‑139 (MQ=255)                                                                                                                 aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTACCTGCGCGAAGCCGCCAtt                                                                                           <  1:110658/55‑1 (MQ=255)                                                                                                                 aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCAtt                                                                                           >  2:110658/1‑55 (MQ=255)                                                                                                                 aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat       >  2:744490/1‑139 (MQ=255)                                                                                                                 aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat       >  2:43273/1‑139 (MQ=255)                                                                                                                 aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat       >  2:520655/1‑139 (MQ=255)                                                                                                                 aTTCCTTAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTAATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat       >  2:472533/1‑139 (MQ=255)                                                                                                                      ttAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAgccg                                                                                                <  2:80607/45‑1 (MQ=255)                                                                                                                      ttAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAgccg                                                                                                >  1:80607/1‑45 (MQ=255)                                                                                                                      ttAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCg  >  2:238690/1‑139 (MQ=255)                                                                                                                      ttAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCg  >  1:72231/1‑139 (MQ=255)                                                                                                                      ttAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCg  >  1:69065/1‑139 (MQ=255)                                                                                                                      ttAATGAAAAATTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCg  >  1:490034/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                TATCCGGCCTACATGTCGCACGAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCG  >  NC_000913/2842445‑2842700

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