New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	526259 =	NA (NA)	7 (0.050)	4/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                                                                                                    |                                                                                                    gtttacgcgacagcgacatcccccccgtcaggtcacgctcccgccgccggacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctc                                                                                                     <  1:2743539‑M2/1‑1 (MQ=255) gtttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctc                                                                                                     <  1:733484‑M2/1‑1 (MQ=255) gtttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctc                                                                                                     <  2:842483‑M2/1‑1 (MQ=255) gtttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctc                                                                                                     <  1:2340208‑M2/1‑1 (MQ=255)  tttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCg                                                                                                    <  2:2209830‑M2/2‑1 (MQ=255)  tttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCg                                                                                                    <  1:1742347‑M2/2‑1 (MQ=255)  tttacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCg                                                                                                    <  1:588350‑M2/2‑1 (MQ=255)    tacgcgacagcgacatccacctcgtcagttcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCg                                                                                                    <  1:889747‑M2/2‑1 (MQ=255)    tacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCGAc                                                                                                  <  2:3381461‑M2/4‑1 (MQ=255)    tacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCGAc                                                                                                  <  2:316271‑M2/4‑1 (MQ=255)    tacgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCGAc                                                                                                  <  2:3423475‑M2/4‑1 (MQ=255)     acgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCGAcc                                                                                                 <  2:2456888‑M2/5‑1 (MQ=255)     acgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCGAcc                                                                                                 <  1:1882490‑M2/5‑1 (MQ=255)     acgcgacagcgacatccaccccgtcaggtcacgctcccgccgccagacccgttttgccattcgccgtcccagcgggccgtagaggtagcggctctCGAcc                                                                                                 >  1:3672785‑M2/96‑100 (MQ=255)     acgcgacagcgacatccaccccgtcagggcacgcgcccgccgcctgacccgttttgccattcgccgtcccagcgggtcgtagaggtagcggctctCGAcc                                                                                                 <  2:1801836‑M2/5‑1 (MQ=255)                                                                          cagcgggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGAtgtcgctgt                            >  1:411030‑M2/27‑100 (MQ=255)                                                                             cgggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGAtgtcgctgtcgc                         >  2:2094136‑M2/24‑100 (MQ=255)                                                                             cgggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGAtgtcgctgtcgc                         >  1:2655432‑M2/24‑100 (MQ=255)                                                                              gggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTcgcg                        >  1:211350‑M2/23‑100 (MQ=255)                                                                              gggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTcgcg                        >  2:1933994‑M2/23‑100 (MQ=255)                                                                               ggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGt                       >  1:1053932‑M2/22‑100 (MQ=255)                                                                               ggtcgtagaggtagcggctctCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGt                       >  2:1642810‑M2/22‑100 (MQ=255)                                                                                                    CGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGTAAACCGGAGGTGACGTGGTa   >  1:1589065‑M2/1‑99 (MQ=255)                                                                                                    CGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGTAAACCGGAGGTGACGTGGTAt  <  2:1105030‑M2/100‑1 (MQ=255)                                                                                                    CGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGTAAACCGGAGGTGACGTGGTAt  >  1:668668‑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 < 22 ≤ 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.