New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	526259 =	NA (NA)	9 (0.080)	6/168	NT	100%	coding (2999/4281 nt)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	526275 =	0 (0.000)					coding (3015/4281 nt)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

CGAGCGGACCCACCACTACCACTACGACAGCCAGCACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGTAAACCGGAGGTGACGTGGTAT  >  NC_000913/526176‑526374                                                                                                    |                                                                                                    gtttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccggcgggtcgtagaggtagcggctctc                                                                                                     <  2:1871657‑M2/1‑1 (MQ=255) gtttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctc                                                                                                     >  1:1121431‑M2/100‑100 (MQ=255) gtttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctc                                                                                                     >  2:1276139‑M2/100‑100 (MQ=255) gtttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctc                                                                                                     <  2:358268‑M2/1‑1 (MQ=255) gtttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctc                                                                                                     <  2:2751718‑M2/1‑1 (MQ=255) gtttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctc                                                                                                     >  2:1034849‑M2/100‑100 (MQ=255) gtttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctc                                                                                                     <  1:2108462‑M2/1‑1 (MQ=255)  tttacgctacagcgacatccacgccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagctggtcgtagaggtagcggctctCg                                                                                                    <  1:704751‑M2/2‑1 (MQ=255)  tttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCg                                                                                                    <  1:132378‑M2/2‑1 (MQ=255)   ttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCGa                                                                                                   <  2:1654886‑M2/3‑1 (MQ=255)   ttacgcgacagcgacagccgccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCGa                                                                                                   <  1:328873‑M2/3‑1 (MQ=255)    tacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCGAc                                                                                                  <  1:1070575‑M2/4‑1 (MQ=255)    tacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgccccagcgggtcgtagaggtagcggctctCGa                                                                                                   <  1:1161532‑M2/3‑1 (MQ=255)     acgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCGAcc                                                                                                 <  1:2108575‑M2/5‑1 (MQ=255)     acgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCGAcc                                                                                                 <  2:1818700‑M2/5‑1 (MQ=255)     acgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCGAc                                                                                                  <  1:3085698‑M2/4‑1 (MQ=255)      cgcgacagctagatccaccccgtcaggtcgcgctcccgccgccagacccgttttgctattcgccgtcccagcgggtcgtagaggtagcggctctCGAc                                                                                                  <  2:2244328‑M2/4‑1 (MQ=255)                                                                     cgtcccagcgggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGAtgtc                                 >  1:2707018‑M2/32‑100 (MQ=255)                                                                       tcccagcgggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGAtgtcgc                               >  2:2080627‑M2/30‑100 (MQ=255)                                                                             cgggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGAtgtcgctgtcgc                         >  1:509219‑M2/24‑100 (MQ=255)                                                                             cgggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGAtgtcgctgtcgc                         >  1:1163779‑M2/24‑100 (MQ=255)                                                                             cgggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGACTGGCGGCGGGAGCGTGAACTGACGGGGTGGAtgtcgctgtcg                          >  1:1172633‑M2/24‑99 (MQ=255)                                                                              gggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTcgcg                        >  2:2869516‑M2/23‑100 (MQ=255)                                                                              gggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGGGGAtgtcgctgtcg                          >  1:2087390‑M2/23‑98 (MQ=255)                                                                               ggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGt                       >  2:75688‑M2/22‑100 (MQ=255)                                                                                         gtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGTAAACCGGAgg             <  2:1163783‑M2/89‑1 (MQ=255)                                                                                                    CGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGTAAACCGGAGGTGACGTGGTAt  <  1:2455716‑M2/100‑1 (MQ=255)                                                                                                    CGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGTAAACCGGAGGTGACGGGGTAt  >  1:3075597‑M2/1‑100 (MQ=255)                                                                                                    |                                                                                                    CGAGCGGACCCACCACTACCACTACGACAGCCAGCACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGTAAACCGGAGGTGACGTGGTAT  >  NC_000913/526176‑526374

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