New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3622174	NA (NA)	5 (0.040)	3/168	NT	100%	coding (2983/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	= 3622205	0 (0.000)					coding (3014/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

AGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATAC  >  NC_000913/3622106‑3622303                                                                                                    |                                                                                                   aGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGT                                                                                                    >  1:1235124‑M2/1‑100 (MQ=255) aGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGT                                                                                                    <  2:2671041‑M2/100‑1 (MQ=255)                                                          gCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgca                                           <  2:3155709‑M2/100‑58 (MQ=255)                                                               cTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctctgcatattgtgtccgcgtgtagtgcaccagc                                      <  1:2266872‑M2/100‑63 (MQ=255)                                                                 ttACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccg                                    <  1:2521609‑M2/100‑65 (MQ=255)                                                                      aCCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcggctttcgaccagcggctctgcatattgtgtccgcgtgtagtgcaccagccggtgct                               <  1:1948817‑M2/100‑70 (MQ=255)                                                                      aCCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgct                               <  1:948841‑M2/100‑70 (MQ=255)                                                                                                cGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagt      <  1:2557427‑M2/99‑95 (MQ=255)                                                                                                cGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGa     >  2:3140969/1‑100 (MQ=31)                                                                                                 gggTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtcc    >  1:150290‑M2/1‑4 (MQ=255)                                                                                                  ggTcgtaaagatagcgactttcgaccagcggatcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccg   <  1:1603725‑M2/100‑98 (MQ=255)                                                                                                   gTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgc  >  1:2645916‑M2/1‑2 (MQ=255)                                                                                                   gTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgc  <  1:3074471‑M2/100‑99 (MQ=255)                                                                                                   gTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgc  >  2:1342866‑M2/1‑2 (MQ=255)                                                                                                   gTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccg   <  1:1116007‑M2/99‑98 (MQ=255)                                                                                                    |                                                                                                   AGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATAC  >  NC_000913/3622106‑3622303

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