New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3765531 =	NA (NA)	13 (0.160)	3/168	NT	100%	coding (3349/4134 nt)	rhsA	Rhs protein with putative toxin 55 domain; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor
	?	NC_000913	3765547 =	0 (0.000)					coding (3365/4134 nt)	rhsA	Rhs protein with putative toxin 55 domain; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor

GGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGG  >  NC_000913/3765449‑3765637                                                                                                   |                                                                                           tcgcagtggtacaggtggatttttctcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggc                                                                                           <  1:1106679‑M2/2‑1 (MQ=255)  cgcagtgttacaggtggatttttcgcgccggcgtgtacaccggctgcatctggttttgcatctgctccacggtcaggtctcacgatgccagccagcggag                                                                                          <  2:1100639‑M2/3‑3 (MQ=255)  cgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcg                                                                                          <  1:2170717‑M2/3‑1 (MQ=255)  cgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcg                                                                                          >  2:534575‑M2/98‑100 (MQ=255)  cgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcg                                                                                          <  1:183094‑M2/3‑1 (MQ=255)  cgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcg                                                                                          >  1:1777973‑M2/98‑100 (MQ=255)  cgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcg                                                                                          >  1:1745821‑M2/98‑100 (MQ=255)   gcagtggtacaggtggatttttcgcgctggcgtgtacaccgggtccatatggttttgcatatggtcaacggtcaggccgcacgatgccagccagcggc                                                                                           >  2:71846‑M2/97‑98 (MQ=255)    cagtggtacaggtggctttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcggc                                                                                        <  1:2362455‑M2/5‑1 (MQ=255)     agtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcggc                                                                                        <  1:779302‑M2/5‑1 (MQ=255)                                                                       cggtcagaccgcacgatgctagcgagggGCGGCCTGACCGTGGAGCAGAGGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGt                     >  2:2049414‑M2/29‑100 (MQ=255)                                                                         gtcaggccgcacgatgccagccagcgGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTAc                   >  1:1367123‑M2/27‑100 (MQ=255)                                                                         gtcaggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTAc                   >  1:262816‑M2/27‑100 (MQ=255)                                                                         gtcaggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTCCACGCCGGCGCGAAAaatccacataaac                   >  1:1592110‑M2/27‑94 (MQ=255)                                                                          cggtcaggccgcacgatgccagccaGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTAcc                  >  2:2187622‑M2/26‑100 (MQ=255)                                                                           caggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCa                 >  1:646792‑M2/25‑100 (MQ=255)                                                                            aggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCAc                >  1:1655713‑M2/24‑100 (MQ=255)                                                                            aggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCAc                >  2:2329646‑M2/24‑100 (MQ=255)                                                                             tcaggccgcacggtgccagccaGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt               >  1:1179073‑M2/23‑100 (MQ=255)                                                                             ggccgcacgatgccagccagcgGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACt               >  1:828534‑M2/23‑100 (MQ=255)                                                                             ggccgcacgatgccagccagcgGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt               >  1:414414‑M2/23‑100 (MQ=255)                                                                              gccgcacgatgccagccagcgGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCAc                >  1:2159132‑M2/22‑98 (MQ=255)                                                                              gccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTg              >  1:363471‑M2/22‑100 (MQ=255)                                                                                          tgccagccaGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCgg  <  2:1179062‑M2/91‑1 (MQ=255)                                                                                                   |                                                                                           GGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGG  >  NC_000913/3765449‑3765637

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.