New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3621534	NA (NA)	7 (0.070)	4/162	NT	100%	coding (2343/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	= 3621592	0 (0.000)					coding (2401/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

GCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCATGAAACTCGGCGACACACCGCTGGTGGA  >  NC_000913/3621486‑3621632                                                 |                                                                                                   gCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGGAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGC                                                 <  1:1555724‑M2/100‑1 (MQ=255)  ctctGGCAGCATGAGACCAGACATGCGTACAACTCGCAGGGGCTGGCGAACCGCTGTATACCGGACGGCCTGCCCGCCGTGGAATGGCTGACCTACGGCA                                                <  1:1646010‑M2/100‑1 (MQ=255)   tctGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGGATACCGAACAGCCTGCCCGCCCTGGAATGGCTGACCTACGGCA                                                <  1:2020433‑M2/99‑1 (MQ=255)   tctGGCAGCATCAGACCAGACATGCGTACAACGCGCAGGGGCGGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAG                                               >  1:251695‑M2/1‑100 (MQ=255)    ctGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGC                                              >  1:1017967‑M2/1‑100 (MQ=255)    ctGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGC                                              >  1:2278952‑M2/1‑100 (MQ=255)    ctGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGC                                              <  1:875677‑M2/100‑1 (MQ=255)    ctGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGC                                              >  1:408725‑M2/1‑100 (MQ=255)     tGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCG                                             >  2:855711‑M2/1‑100 (MQ=255)     tGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCG                                             >  2:1084096‑M2/1‑100 (MQ=255)      ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                            <  1:666520‑M2/100‑1 (MQ=255)      ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                            <  2:1482788‑M2/100‑1 (MQ=255)      ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                            <  2:1934992‑M2/100‑1 (MQ=255)      ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                            >  2:2113970‑M2/1‑100 (MQ=255)      ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                            >  2:2600242‑M2/1‑100 (MQ=255)       gcagcaTGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGT                                           <  2:564441‑M2/100‑1 (MQ=255)       gcagcaTGAGACCAGACATGCGTACAACGCACAGGGGCTGGCGAACGGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGggg                                            <  1:970769‑M2/99‑4 (MQ=255)        cagcaTGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTT                                          <  1:1961120‑M2/100‑1 (MQ=255)        cagcaTGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTT                                          >  2:2593298‑M2/1‑100 (MQ=255)                                               gCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTcgccagcccctgcgcgttgtacgcatgtctggtctcatg   <  1:358100‑M2/100‑40 (MQ=255)                                                cGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTcgccagcccctgcgcgttgtacgcatgtctggtctcatgc  <  1:878856‑M2/100‑41 (MQ=255)                                                cGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTcgccagcccctgcgcgttgtacgcatgtctggtctcatgc  <  2:377119‑M2/100‑41 (MQ=255)                                                cGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTcgccagcccctgcgcgttgtacgcatgtctggtctcatgc  <  1:1136958‑M2/100‑41 (MQ=255)                                                 |                                                                                                   GCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCATGAAACTCGGCGACACACCGCTGGTGGA  >  NC_000913/3621486‑3621632

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