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

GCCAGCACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGTAAACCGGAGGTGACGTGGTA  >  NC_000913/526205‑526373                                                                       |                                                                                                   gccagcACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACG                                                                       <  1:1129855‑M2/100‑1 (MQ=255) gccagcACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACG                                                                       <  1:11505‑M2/100‑1 (MQ=255) gccagcACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACG                                                                       <  1:335207‑M2/100‑1 (MQ=255) gccagcACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACG                                                                       >  2:2596830‑M2/1‑100 (MQ=255)  ccagcACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACGG                                                                      >  1:2254563‑M2/1‑100 (MQ=255)  ccagcACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACGG                                                                      <  1:393662‑M2/100‑1 (MQ=255)   cagcACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACGGG                                                                     >  1:3027755‑M2/1‑100 (MQ=255)   cagcACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACGGG                                                                     >  2:525667‑M2/1‑100 (MQ=255)   cagcACCGCCCGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACGGG                                                                     >  1:2863104‑M2/1‑100 (MQ=255)    agcACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACGGGT                                                                    <  1:2984914‑M2/100‑1 (MQ=255)     gcACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACGGGTC                                                                   >  1:2225597‑M2/1‑100 (MQ=255)     gcACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACGGGTC                                                                   >  2:735562‑M2/1‑100 (MQ=255)                                                                 ccTCTACGACCCGCTGGGACGGCGAATGGCAAAACGGGTCgtagaggtagcggctctcgaccagtggctcgccatgctgtatccgcgtgtagaacaccag       <  1:2205146‑M2/100‑61 (MQ=255)                                                                    ctACGACCCGCTGGGACGGCGAATGGCAAAACGGGTCgtagaggtagcggctctcgaccagtggctcgccatgctgtatccgcgtgtagaacaccaggcg    <  1:1251267‑M2/100‑64 (MQ=255)                                                                      aCGACCCGCTGGGACGGCGAATGGCAAAACGGGTCgtagaggtagcggctctcgaccagtggctcgccatgctgtatccgcgtgtagaacaccaggcggt  <  2:807677‑M2/100‑66 (MQ=255)                                                                       |                                                                                                   GCCAGCACCGCCTGGTGTTCTACACGCGGATACAGCATGGCGAGCCACTGGTCGAGAGCCGCTACCTCTACGACCCGCTGGGACGGCGAATGGCAAAACGGGTCTGGCGGCGGGAGCGTGACCTGACGGGGTGGATGTCGCTGTCGCGTAAACCGGAGGTGACGTGGTA  >  NC_000913/526205‑526373

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