New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4356370 =	172 (2.180)	4 (0.050)	4/250	NT	2.3%	coding (42/1458 nt)	dtpC	dipeptide and tripeptide permease
	?	NC_000913	4356389 =	174 (2.230)					coding (23/1458 nt)	dtpC	dipeptide and tripeptide permease
Rejected: Frequency below/above cutoff threshold.

TGATGGCCGCTATACCGTTAAGGTATTGAAAGAAGAAAGCAAAAAATAATTAGCTCGTACAAGGGAAGTGGCTTGCCACTTCCCTTTTTTGCTCATAAGGAGAACACATGAAAACACCCTCACAGCCGCGCGCGATATACTATATCGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4356518‑4356370 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGA  >  NC_000913/4356389‑4356426                                                                                                                                                                                             TGATGGCCGCTATACCGTTAAGGTATTGAAAGAAGAAAGCAAAAAATAATTAGCTCGTACAAGGGAAGTGGCTTGCCACTTCCCTTTTTTGCTCATAAGGAGAACACATGAAAACACCCTCACAGCCGCGCGCGATATACTATATCGTG                                        >  5:50321/1‑149   ATGGCCGCTATACCGTTAAGGTATTGAAAGAAGAAAGCAAAAAATAATTAGCTCGTACAAGGGAAGTGGCTTGCCACTTCCCTTTTTTGCTCATAAGGAGAACACATGAAAACACCCTCACAGCCGCGCGCGATATACTATATCGTG                                        <  7:46385/147‑1        CGCTATACCGTTAAGGTATTGAAAGAAGAAAGCAAAAAATAATTAGCTCGTACAAGGGAAGTGGCTTGCCACTTCCCTTTTTTGCTCATAAGGAGAACACATGAAAACACCCTCACAGCCGCGCACGATATAGTATATCGCGCGCGGCT                                 >  5:333043/1‑149                                       GCAAAAAATAATTAGCTCGTACAAGGGAAGTGGCTTGCCACTTCCCTTTTTTGCTCATAAGGAGAACACATGAAAACACCCTCACAGCCGCGCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGA  <  6:333043/149‑1                                                                         TTGCCACTTCCCTTTTTTGCTCATAAGGAGAACACATGAAAACACCCTCACAGCCGCGCGCGATATACTATATCGTG                                        >  3:347493/1‑77                                                                         TTGCCACTTCCCTTTTTTGCTCATAAGGAGAACACATGAAAACACCCTCACAGCCGCGCGCGATATACTATATCGTG                                        <  4:347493/77‑1                                                                                                                        TCACAGCCGCGCGCGATATACTATATCGTGCGCGGTTGTGAGGGTGTTTTCATGT               >  7:258837/1‑55                                                                                                                        TCACAGCCGCGCGCGATATACTATATCGTGCGCGGTTGTGAGGGTGTTTTCATGT               <  8:258837/55‑1                                                                                                                                                                                             TGATGGCCGCTATACCGTTAAGGTATTGAAAGAAGAAAGCAAAAAATAATTAGCTCGTACAAGGGAAGTGGCTTGCCACTTCCCTTTTTTGCTCATAAGGAGAACACATGAAAACACCCTCACAGCCGCGCGCGATATACTATATCGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4356518‑4356370 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGA  >  NC_000913/4356389‑4356426

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