New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1390544 =	41 (0.600)	4 (0.060)	4/202	6.7	8.5%	coding (55/729 nt)	mpaA	murein peptide amidase A
	?	NC_000913	1390553 =	45 (0.670)					coding (46/729 nt)	mpaA	murein peptide amidase A
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGCCTGTAATACTTCCTGAACCAGAGACGATCATGACCGTAACCCGCCCACGCGCCGAACGCGGCGCATTTCCGCCCGGAACAGAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1390630‑1390544 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cCGGGCGGAAATGCGCCGCGTTCGGCGCG  >  NC_000913/1390553‑1390580                                                                                                                     TGCCTGTAATACTTCCTGAACCAGAGACGATCATGACCGTAACCCGCCCACGCGCCGAACGCGGCGCATTTCCGCCCGGAACAGAACCGGGCGGAAATGCGCCGCGTTCGGCGCG  >  1:952352/1‑115 TGCCTGTAATACTTCCTGAACCAGAGACGATCATGACCGTAACCCGCCCACGCGCCGAACGCGGCGCATTTCCGCCCGGAACAGAACCGGGCGGAAATGCGCCGCGTTCGGCGCG  <  2:952352/115‑1          TACTTCCTGAACCAGAGACGATCATGACCGTAACCCGCCCACGCGCCGAACGCGGCGCATTTCCGCCCGGAACAGAAC‑‑GGCGGAAATGCGCCGCGTT         <  1:498485/97‑1          TACTTCCTGAACCAGAGACGATCATGACCGTAACCCGCCCACGCGCCGAACGCGGCGCATTTCCGCCCGGAACAGAAC‑‑GGCGGAAATGCGCCGCGTT         >  2:498485/1‑97                                                                                                                     TGCCTGTAATACTTCCTGAACCAGAGACGATCATGACCGTAACCCGCCCACGCGCCGAACGCGGCGCATTTCCGCCCGGAACAGAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1390630‑1390544 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cCGGGCGGAAATGCGCCGCGTTCGGCGCG  >  NC_000913/1390553‑1390580

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