New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_012971	550553 =	0 (0.000)	44 (0.440)	18/588	4.1	100%	pseudogene (4/364 nt)	ECD_RS02680	prophage tail fiber N‑terminal domain‑containing protein
	?	NC_012971	= 788130	0 (0.000)					pseudogene (6/1155 nt)	ECD_RS23105	prophage tail fiber N‑terminal domain‑containing protein
Rejected: Coverage evenness skew score above cutoff.

TACGGGCGCCAGAAATGGCGCCTTTTTTATTTGTGGAGTGAATATG‑‑‑GCAGTACAGATTTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTGGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTTGAGCATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCAGGTACGC  >  NC_012971/550507‑550805                                                  |                                                                                                                                                                                                                                                             atgacagcccgccggttcaggcgggcttttttGTGGAGTGGATATGGCAGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTAc    <  1:768578/271‑1 (MQ=25)        cccgccggttcaggcgggcttttttGTGGAGTGGATATGGCAGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGc  <  1:609414/273‑1 (MQ=25)        cccgccggttcaggcgggcttttttGTGGAGTGGATATGGCAGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCACTACCGTCTATGAAGGTTCCAGACCTGGTACGc  <  1:203931/273‑1 (MQ=25)        cccgccggttcaggcgggcttttttGTGGAGTGGATATGGCAGCAGTAAAAATCTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTAGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTCGAGTATGGCAAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCTGGTACGc  <  1:77083/273‑1 (MQ=25)                                                  |                                                                                                                                                                                                                                                             TACGGGCGCCAGAAATGGCGCCTTTTTTATTTGTGGAGTGAATATG‑‑‑GCAGTACAGATTTCAGGTGTGCTGAAAGATGGTGCGGGAAAACCAATACAGAACTGCACTATTCAACTGAAGGCAAAGCGTAACAGCACCACGGTACTGGTGAACACGGTGGCCTCTGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTTGAGCATGGCCAGTACAGCGTCACCCTGCTGGTTGAAGGTTTTCCGCCTTCACATGCCGGAACCATTACCGTCTATGAAGGTTCCAGACCAGGTACGC  >  NC_012971/550507‑550805

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