Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
JC	NC_000913	2,966,777	Δ4 bp	1.8%	coding (1655‑1658/2247 nt)	ptsP ←	PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2966776	167 (1.500)	3 (0.030)	3/546	NT	1.8%	coding (1659/2247 nt)	ptsP	PEP‑protein phosphotransferase enzyme I; GAF domain containing protein
	?	NC_000913	2966781 =	169 (1.530)					coding (1654/2247 nt)	ptsP	PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

TCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2966478‑2966776 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACT  >  NC_000913/2966781‑2966990                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               TCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATAGG                                                                                                                                                                                                                  <  2:397444/301‑1                                                                                                     GGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGACTAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCCCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGC                                                                                                                 >  1:828937/1‑298                                                                                                                                                                                                                  ATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATAC   <  2:420629/299‑1                                                                                                                                                                                                                  ATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGAAACc  >  1:420622/1‑299                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               TCGCCAGAGCTCGTAACATTGCAGGATGAAGACTGTCATAAATGTTCGCCACCCGGGTATTGTTGCGATCAACGGCCAGAATGTATTGAGTCAGATCGTTGGTGCCAACAGAGATGAAATCGACCCGCTTTGCCAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2966478‑2966776 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑aGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACT  >  NC_000913/2966781‑2966990

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ 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.