New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3620750 =	NA (NA)	4 (0.130)	3/128	2.8	100%	coding (1559/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	3621165 =	0 (0.000)					coding (1974/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

AAGCTCAGCAACTGACCGTAACGGCTCCACGTCATGGTTTTCCGGCTGCCGGTGGCATCTTCCGTTGcg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3620816‑3620750 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCT  >  NC_000913/3621165‑3621233                                                                                                                                          AAGCTCAGCAACTGACCGTAACGGCTCCACGTCATGGTTTTCCGGCTGCCGGTGGCATCTTCCGTTGC                                                                      >  1:444055/1‑68 AAGCTCAGCAACTGACCGTAACGGCTCCACGTCATGGTTTTCCGGCTGCCGGTGGCATCTTCCGTTGC                                                                      <  1:634462/68‑1 AAGCTCAGCAACTGACCGTAACGGCTCCACGTCATGGTTTTCCGGCTGCCGGTGGCATCTTCCGTTGC                                                                      >  1:953605/1‑68 AAGCTCAGCAACTGACCGTAACGGCTCCACGGCATGGTTTTCCGGCTGCCGGTGGCATCTTCCGTTGC                                                                      <  2:944323/68‑1                       GGCTCCACGTCATGGTTTTCCGGCTGCCGGTGGCATCTTCCGTTGCGATGTACTTGACCGGCTGATAC                                                <  1:318199/68‑1                       GGCTCCACGTCATGGTTTTCCGGCTGCCGGTGGCATCTTCCGTTGCGATGTACTTGACCGGCTGATAC                                                >  1:991217/1‑68                                    GGTTTTCCGGCTGCCGGTGGCATCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTT                                   >  1:532354/1‑68                                                         acCTTCCGTTGCGATGTACTTGACCGGCTGATACAGGAAGCCGGCTTTGACGGCCGCACACAGCGTTA              >  2:263588/3‑68                                                                     GATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCT  >  1:130763/1‑68                                                                                                                                          AAGCTCAGCAACTGACCGTAACGGCTCCACGTCATGGTTTTCCGGCTGCCGGTGGCATCTTCCGTTGcg‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3620816‑3620750 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCT  >  NC_000913/3621165‑3621233

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