New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3620868	NA (NA)	5 (0.120)	5/276	2.8	13.1%	coding (1677/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	= 3621008	36 (0.830)					coding (1817/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Rejected: Frequency below cutoff threshold.

GACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATT  >  NC_000913/3620848‑3621009                     |                                                                                                                                              ttctgctgccgtccggggca‑atgacggcggtgcACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGAC            >  2:177019‑M2/16‑151 (MQ=255) gACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTACTTGCCGTGAAACACACGCAGGGCCATGAAACGCGGCATGAATCCAACATCGCCGGTGACCTGA             <  1:424357‑M2/151‑1 (MQ=255) gACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGGAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGA             <  2:214011‑M2/151‑1 (MQ=255) gACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGA             <  1:284498‑M2/151‑1 (MQ=255) gACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGA             <  1:404507‑M2/151‑1 (MQ=255) gACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGA             <  2:623149‑M2/151‑1 (MQ=255)  aCCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGA             <  1:575412‑M2/150‑1 (MQ=255)  aCCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGA             <  2:461227‑M2/150‑1 (MQ=255)   ccATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGAC            <  1:322954‑M2/150‑1 (MQ=255)            cTTCGGGCAGAAGACGGCGGTGCACCGCGAGAAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGCCACGCAGGGCCATGAACCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATc  >  1:404806‑M2/1‑150 (MQ=255)                     |                                                                                                                                              GACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATT  >  NC_000913/3620848‑3621009

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