Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	3,382,371	(GGT)2→1	3.6%	coding (174‑176/1068 nt)	degS →	serine endoprotease, periplasmic

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3382370	108 (1.010)	4 (0.040)	4/424	NT	3.6%	coding (173/1068 nt)	degS	serine endoprotease, periplasmic
	?	NC_000913	3382374 =	107 (1.020)					coding (177/1068 nt)	degS	serine endoprotease, periplasmic

TCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGC‑GG‑T‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3382143‑3382370 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gg‑tTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGT  >  NC_000913/3382374‑3382608                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   TCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGCCGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCCGGCTGTTCGCCGCGCCGCGCCAGC‑Gctggcttacgctt                                                                                                                                                                                                                                   >  1:25133/1‑226                                                                                                                                                     CTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGC‑GG‑TTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCAT                                                                                                                                          <  1:628288/179‑1                                                                                                                                                     CTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGC‑GG‑TTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCAT                                                                                                                                          >  2:628288/1‑179                                                                                                                                                               CGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGCGGG‑TTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCTTGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGA                                                                     >  2:66770/1‑239                                                                                                                                                                 cgATTTGACTGTACCGAAAGGTCGCCTGCCTGCTATAGTCTTGCGTATCGCCGCGCCGCGCCAGC‑GGTTTAACGTTTACAACCGGAGTTTGAACACCATCTCTCACAACCAGCTTGAGATTCGCACCCTGGGATCCGGTGTAATCATTGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAG                                                                   <  1:66770/237‑1                                                                                                                                                                                                                                    G‑TTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGT  >  1:423148/1‑237                                                                                                                                                                                                                                    G‑TTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGT  <  2:423148/237‑1                                                                                                                                                                                                                                    G‑TTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACT        >  1:517361/1‑231                                                                                                                                                                                                                                    G‑TTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACT        <  2:517361/231‑1                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   TCCGGTTAACTCGTGGTATGCTGCTGCCGTTCCCTTTTTTAATGACGCCTCCATCATGTTTGTGAAGCTCTTACGTTCCGTTGCGATTGGATTAATTGTCGGCGCTATTCTGCTGGTTGCCATGCCTTCGCTGCGCAGCCTTAACCCGCTTTCCACTCCGCAATTTGACAGTACCGATGAGACGCCTGCCAGCTATAATCTGGCGGTTCGCCGCGCCGCGCCAGC‑GG‑T‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/3382143‑3382370 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gg‑tTAACGTTTACAACCGTGGTTTGAACACCAACTCTCACAACCAGCTTGAGATCCGCACCCTGGGATCCGGTGTAATCATGGATCAACGCGGTTATATCATCACCAATAAACACGTCATCAACGACGCCGATCAGATCATCGTCGCCTTACAGGATGGACGTGTATTTGAAGCATTGCTGGTGGGATCTGACTCTCTAACCGATCTGGCGGTACTTAAAATTAATGCCACTGGCGGT  >  NC_000913/3382374‑3382608

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.