New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3765534 =	NA (NA)	19 (0.220)	4/168	NT	100%	coding (3352/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

GGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCG  >  NC_000913/3765449‑3765627                                                                                                   |                                                                                 tggtcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagc                                                                                 >  2:2064190‑M2/99‑100 (MQ=255) tggtcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagc                                                                                 <  1:1233372‑M2/2‑1 (MQ=255)  ggtcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccaGCg                                                                                >  1:1476491‑M2/98‑100 (MQ=255)   gtcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccaGCgg                                                                               >  2:405919‑M2/97‑100 (MQ=255)    tcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccaGCGGc                                                                              <  2:1395408‑M2/5‑1 (MQ=255)    tcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccaGCGGc                                                                              <  1:821796‑M2/5‑1 (MQ=255)    tcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccaGCGGc                                                                              >  2:389456‑M2/96‑100 (MQ=255)                                                                     ctccacggtcaggccgcacgatgccagccaGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCt             >  1:1833492‑M2/31‑100 (MQ=255)                                                                     cacggtcaggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCt             >  1:1046074‑M2/31‑100 (MQ=255)                                                                      accgtcaggccgcacgatgccagccagcgGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTg            >  2:937818‑M2/30‑100 (MQ=255)                                                                      accgtcaggccgcacgatgccagccagcgGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTg            >  2:2071519‑M2/30‑100 (MQ=255)                                                                        cgtcaggccgcacgatgccagccagcgGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTa          >  1:1004128‑M2/28‑100 (MQ=255)                                                                           ggtcaggccgcacgatgccagccaGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCa       >  1:1289766‑M2/25‑100 (MQ=255)                                                                            gtcaggccgcacgatgccagccaGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCAc      >  1:2176152‑M2/24‑100 (MQ=255)                                                                            gtcaggccgcacgatgccagccaGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAATATCCACCTGTACCAc      >  2:2558124‑M2/24‑100 (MQ=255)                                                                             tcaggccgcacgatgccagccaGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACt     >  2:507011‑M2/23‑100 (MQ=255)                                                                             ggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt     >  2:1124429‑M2/23‑100 (MQ=255)                                                                             ggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt     >  1:2525574‑M2/23‑100 (MQ=255)                                                                             ggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt     >  1:251628‑M2/23‑100 (MQ=255)                                                                             ggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt     >  1:1405966‑M2/23‑100 (MQ=255)                                                                             ggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt     >  1:1033283‑M2/23‑100 (MQ=255)                                                                             ggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGCTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt     >  2:2470330‑M2/23‑100 (MQ=255)                                                                             ggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt     >  2:248925‑M2/23‑100 (MQ=255)                                                                             ggccgcacgatgccagccagcgGCGGCCTGACCGTGAAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt     >  2:1829720‑M2/23‑100 (MQ=255)                                                                              caggccgcacgatgccagccaGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTg    <  2:1833502‑M2/79‑1 (MQ=255)                                                                                cgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCg  >  2:1124115‑M2/20‑100 (MQ=255)                                                                                                   |                                                                                 GGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCG  >  NC_000913/3765449‑3765627

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