New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3765531 =	NA (NA)	12 (0.150)	5/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

GTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCG  >  NC_000913/3765450‑3765627                                                                                                  |                                                                                 cgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcg                                                                                >  2:603212‑M2/98‑100 (MQ=255)  gcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcgg                                                                               <  1:1356913‑M2/4‑1 (MQ=255)  gcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcgg                                                                               <  2:232716‑M2/4‑1 (MQ=255)   cagttgtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcgg                                                                               <  2:536952‑M2/4‑1 (MQ=255)   cagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcggc                                                                              <  1:1442579‑M2/5‑1 (MQ=255)   cagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcggc                                                                              >  1:2347956‑M2/96‑100 (MQ=255)   cagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcggc                                                                              >  2:1917499‑M2/96‑100 (MQ=255)           caggtgcatttttcgctccggcgtctacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagccgcggc                                                                              <  2:2095088‑M2/5‑1 (MQ=255)                                                               tgctccacggtcaggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATc                  >  2:231469‑M2/36‑100 (MQ=255)                                                                      cggtcaggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGGTGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGt           >  1:482770‑M2/29‑100 (MQ=255)                                                                      ccacggtcaggccgcacgatgccagccaGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGt           >  1:914009‑M2/29‑100 (MQ=255)                                                                       ggtcaggccgcacgatgccagccagcgGCGGCTTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTa          >  1:1055046‑M2/28‑100 (MQ=255)                                                                          caggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCa       >  1:2036905‑M2/25‑100 (MQ=255)                                                                            ggccgcacgatgccagccagcgGCGGCGTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAAACCACCTGTACCAc      >  2:991484‑M2/23‑99 (MQ=255)                                                                            ggccgcacgatgccagccagcgGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACt     >  2:774398‑M2/23‑100 (MQ=255)                                                                            ggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt     >  1:1383009‑M2/23‑100 (MQ=255)                                                                             gccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTg    >  2:624763‑M2/22‑100 (MQ=255)                                                                              ccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGc   >  1:2042506‑M2/21‑100 (MQ=255)                                                                              ccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGc   >  1:1177400‑M2/21‑100 (MQ=255)                                                                               ggccgcacgatgccagccaGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCg  >  1:2256412‑M2/20‑100 (MQ=255)                                                                                                  |                                                                                 GTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCG  >  NC_000913/3765450‑3765627

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