New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3622015 =	0 (0.000)	3 (0.020)	3/512	NT	2.6%	coding (2824/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	3765561 =	223 (1.500)					coding (3379/4134 nt)	rhsA	Rhs protein with putative toxin 55 domain; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor
Rejected: Frequency below/above cutoff threshold.

AGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATAC  >  NC_000913/3621728‑3622303                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                 ggatcctgagtgatataccgcccctgcagcgggtcataatagcggtggcggttgtaatacaggccggactcctcatcatactgctgccccggcaggcggataagctgctgcagctgatgcgggttctcttcattcagcaggttgccccattcatcgtattctgcgcaccattctgttgtcccttccttgctgataagggccagcggcaggccgcgatggtcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctcCa                                                                                                                                                                                                                                                                                                 >  1:802247‑M1/288‑289 (MQ=255)   atcctgagtgatataccgcccctgcagcgggtcataatagcggtggcggttgtaatacaggccggactcctcatcatactgctgccccggcaggcggataagctgctgcagctgatgcgggttctcttcattcagcaggttgccccattcatcgtattctgcgcaccattctgttgtcccttccttgctgataagggccagcggcaggccgcgatggtcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctcCAc                                                                                                                                                                                                                                                                                                >  1:445656‑M1/286‑288 (MQ=255)                                                                                                                                                       tcatcgtattctgcgcaccattctgttgtcccttccttgctgataagggccagcggcaggccgcgatggtcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctcCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCCGGCCGAAAgtcggta                                                                                                                                           <  2:133301‑M1/152‑4 (MQ=255)                                                                                                                                                                                                                   ccgcgatggtcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctcCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGAc                                                                                 >  1:133301‑M1/78‑287 (MQ=255)                                                                                                                                                                                                                                            ggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctcCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACCGGCTGGATGTCGCTGTCACgg                                                      <  1:1250655‑M1/237‑1 (MQ=255)                                                                                                                                                                                                                                                                                              gTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGAt    >  1:137678/1‑289 (MQ=31)                                                                                                                                                                                                                                                                                                CACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAACGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGGCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACGGGCCGCCCACgatca  <  2:1139688‑M1/289‑3 (MQ=255)                                                                                                                                                                                                                                                                                                CACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGAACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGGAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACgatcc  <  1:261761‑M1/289‑3 (MQ=255)                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                 AGCACCTGAACAGCCTGCTGTCTGACCGCGATTACACCTGGAACGACAACGGCGAACTCATCCGCATCAGCAGCCCGCGCCAGACCCGGAGTTACAGCTACAGCACCACCGGCAGGCTGACCGGCGTTCACACCACCGCAGCGAATCTGGATATCCGCATCCCGTATGCCACAGACCCGGCAGGTAACCGCCTGCCCGACCCGGAGCTGCACCCGGACAGCACCCTCAGCATGTGGCCGGATAACCGTATCGCCCGTGACGCGCACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATAC  >  NC_000913/3621728‑3622303

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