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%	115.9 / NA	31	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 (17/14);  total (17/14)

GAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCAC  >  NC_000913/2842466‑2842722                                                                                                                              |                                                                                                                                    gAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAgc                                                                                                                        >  1:542312/1‑139 (MQ=255)             cccGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATttcctt                                                                                                            >  2:96174/1‑139 (MQ=255)               cGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGTGCGAAGCCGCCATttccttcc                                                                                                          <  2:245264/139‑1 (MQ=255)                 aTAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCa                                                                                                        >  2:561734/1‑139 (MQ=255)                 aTAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCa                                                                                                        >  1:40035/1‑139 (MQ=255)                             caTCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCg                                                                                            <  1:366391/139‑1 (MQ=255)                                    cAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATaa                                                                                     <  2:40035/139‑1 (MQ=255)                                      aTTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAAc                                                                                   <  1:96174/139‑1 (MQ=255)                                       tttAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACg                                                                                  >  2:401586/1‑139 (MQ=255)                                         tAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGt                                                                                <  1:722479/139‑1 (MQ=255)                                                tGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTAcc                                                                         >  1:577115/1‑139 (MQ=255)                                                 ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                                        <  2:542312/139‑1 (MQ=255)                                                 ggCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                                        >  2:756035/1‑139 (MQ=255)                                                         agCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTa                                                                >  2:455240/1‑139 (MQ=255)                                                                   ccccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTAcc                                                                         <  2:145503/120‑1 (MQ=255)                                                                   ccccATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTAcc                                                                         >  1:145503/1‑120 (MQ=255)                                                                      cATCAAGCACCTCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCa                                                   <  2:189544/139‑1 (MQ=255)                                                                       aTCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCAt                                                  <  2:577115/139‑1 (MQ=255)                                                                        tCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                                        <  2:37859/116‑1 (MQ=255)                                                                        tCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCg                                                                        >  1:37859/1‑116 (MQ=255)                                                                                 tCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAAATGCGCGAAGCCGCCATTTCTTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTaa                                        >  2:469742/1‑139 (MQ=255)                                                                                  cccTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAgccgcc                                                                                                                    >  2:278900/1‑62 (MQ=255)                                                                                  cccTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAgccgcc                                                                                                                    <  1:278900/62‑1 (MQ=255)                                                                                            aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat                             >  1:90555/1‑139 (MQ=255)                                                                                            aTTCCTTAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGat                             >  2:521668/1‑139 (MQ=255)                                                                                            aTTCCTTAATG‑AAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATttccttcc                                                                                                          <  1:543552/61‑1 (MQ=255)                                                                                            aTTCCTTAATG‑AAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATttccttcc                                                                                                          >  2:543552/1‑61 (MQ=255)                                                                                                 ttAATGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCg                        >  2:165603/1‑139 (MQ=255)                                                                                                    aTGAAAAACTACTCTTCTTCCACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCt                     <  1:521668/139‑1 (MQ=255)                                                                                                                       ccACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCAc  >  2:632708/1‑139 (MQ=255)                                                                                                                       ccACCGCCAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCAc  <  1:238270/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                    GAACCATGAAGGCCCGATAAGCGCAGCGCATCGGGCAATTTAGCATTTGGCATCAGAGCTGGTTCACCCCATCAAGAACATCCCTGTCCTGATTCCTTAATGAAAAACTACTCTTCTTCCACCGCTAACTGCGCGAAGCCGCCATTTCCTTCCCAGCCTTCCAGTCGTTGATAAACGGTTTCTACCGCATCTTTAATCGGCTGGGTCATCGGGTAGTAAAAGCCGACGATATCCGGCTGAATGCCGAGGAAAATCAC  >  NC_000913/2842466‑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: ‑