New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3620569 =	NA (NA)	12 (0.150)	7/166	NT	100%	coding (1378/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	3620582 =	0 (0.000)					coding (1391/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

TCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCA  >  NC_000913/3620483‑3620681                                                                                                    |                                                                                                    caattcatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcc                                                                                                     <  1:1900253‑M2/1‑1 (MQ=255)    ttcatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGc                                                                                                  <  2:561323‑M2/4‑1 (MQ=255)     tcatcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGCg                                                                                                 <  2:948478‑M2/5‑1 (MQ=255)      catcatattcccggcgcaattccagcccgtcaggcccggtggctgccgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGCGt                                                                                                <  2:1264457‑M2/6‑1 (MQ=255)      catcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatcCGGCGt                                                                                                >  1:464805‑M2/95‑100 (MQ=255)      catcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggtgatagtaaaacgccgatgccctgccatcCGGCGt                                                                                                >  1:880993‑M2/95‑100 (MQ=255)       atcatattcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtcaaacgccgatgccctgccatcCg                                                                                                    >  2:1486580‑M2/94‑95 (MQ=255)              tcccggcgcaattccagcccgtcaggcccggtggctgacgttaactggttgtggtggttatagtaaaacgccgatgccctgccatccggcgtcaaactgg                                                                                        >  1:326644‑M2/87‑92 (MQ=255)                                                                         atagtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGa                             >  1:587027‑M2/28‑100 (MQ=255)                                                                         atagtaaaacgccgatgccctgccatcCGGCGGTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCggg                              >  1:2264277‑M2/28‑99 (MQ=255)                                                                          tagtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCggg                              >  2:2209868‑M2/27‑98 (MQ=255)                                                                            gtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAata                          >  1:1470383‑M2/25‑100 (MQ=255)                                                                            gtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAat                           >  2:783721‑M2/25‑99 (MQ=255)                                                                            gtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAata                          >  1:1557950‑M2/25‑100 (MQ=255)                                                                            gtaaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAat                           >  2:331823‑M2/25‑99 (MQ=255)                                                                             taaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAatat                         >  2:1087091‑M2/24‑100 (MQ=255)                                                                              aaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAata                          >  2:207951‑M2/23‑98 (MQ=255)                                                                                aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatga                       >  1:2306940‑M2/21‑99 (MQ=255)                                                                                aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCACCGGGAATatgat                      >  1:671095‑M2/21‑100 (MQ=255)                                                                                 acgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatgatg                     >  2:119202‑M2/20‑100 (MQ=255)                                                                                                    CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTc   <  1:1439240‑M2/99‑1 (MQ=255)                                                                                                    CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTc   <  2:1124358‑M2/99‑1 (MQ=255)                                                                                                    CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCa  >  2:1128137‑M2/1‑100 (MQ=255)                                                                                                    CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCa  <  1:1458415‑M2/100‑1 (MQ=255)                                                                                                    CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCa  >  1:1012164‑M2/1‑100 (MQ=255)                                                                                                    |                                                                                                    TCAGGGCACAGACGGATGCCGCAGGCAGGACAACAGAGTACAGCCCGGATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCA  >  NC_000913/3620483‑3620681

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