New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3620066 =	NA (NA)	3 (0.150)	3/104	2.0	100%	coding (875/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	3620490 =	0 (0.000)					coding (1299/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

GAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCA  >  NC_000913/3620424‑3620557                                                                   |                                                                    agcgcaccagcggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagac                                                                    >  1:211512‑M2/67‑68 (MQ=255)    gcaccagcggctcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacagt                                                                 <  1:510275‑M2/5‑2 (MQ=255)    gcaccagcggcgcggcaggtaagttctccgggtattccgggtcgtgcgtcagccacacggcagacaga                                                                 >  1:346086‑M2/64‑68 (MQ=255)    gcaccagcggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacaga                                                                 <  2:9202‑M2/5‑1 (MQ=255)    gcaccagcggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacaga                                                                 <  2:707940‑M2/5‑1 (MQ=255)        cagcggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagACAGACGGa                                                             <  2:951836‑M2/9‑1 (MQ=255)            ggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagACAGACGGATGcc                                                         <  1:841836‑M2/13‑1 (MQ=255)            ggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagACAGACGGATGcc                                                         <  1:916350‑M2/13‑1 (MQ=255)            ggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagACAGACGGATGcc                                                         >  1:779370‑M2/56‑68 (MQ=255)            ggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagACAGACGGATGcc                                                         <  2:834283‑M2/13‑1 (MQ=255)            ggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagACAGACGGATGcc                                                         <  1:684081‑M2/13‑1 (MQ=255)                         aattctccgggtattccgggtcgtgcgtcagccacacggcagACAGACGGATGCCGCAGGCAGGacaa                                            >  1:954209‑M2/43‑68 (MQ=255)                                  ggtattccgggtcgtgcgtcagccacacggcagACAGACg                                                               <  2:635083‑M2/7‑1 (MQ=255)                                  ggtattccgggtcgtgcgtcagccacacggcagACAGACg                                                               >  1:635083‑M2/34‑40 (MQ=255)                                    tattccggggcgtgcgtcagccacacggcagACAGACGGATGCCGCAGGCAGGACAACAGAGTACAGc                                 <  1:836001‑M2/37‑1 (MQ=255)                                     attccgtgtcgtgcgtcagccacacggcagACAGACGGATGCCTCAGGCAGGACAACAGAGTACAGcc                                <  1:864397‑M2/38‑1 (MQ=255)                                                                   ACAGACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCa  <  2:665965‑M2/68‑1 (MQ=255)                                                                   |                                                                    GAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCA  >  NC_000913/3620424‑3620557

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