New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3621314	NA (NA)	6 (0.040)	5/168	NT	100%	coding (2123/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	= 3621344	0 (0.000)					coding (2153/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

AACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTGGCAGTATGACGAACGTGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGGCGGTGCA  >  NC_000913/3621242‑3621413                                                                         |                                                                                                    aaCTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCG                                                                          <  2:908848‑M2/100‑1 (MQ=255) aaCTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCG                                                                          <  2:4235694‑M2/100‑1 (MQ=255) aaCTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCG                                                                          >  2:2997541‑M2/1‑100 (MQ=255) aaCTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCG                                                                          <  2:1414456‑M2/100‑1 (MQ=255) aaCTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCG                                                                          >  2:1253848‑M2/1‑100 (MQ=255) aaCTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCG                                                                          <  2:1236246‑M2/100‑1 (MQ=255)   cctatCCGCAGCTAGGATGAGTGTCTGGTCACCCACTGGCACTATGACGAAGCAGCCCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGC                                                                           >  1:3712711‑M2/3‑97 (MQ=255)   cTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACGGCATCACGCCCCGCACCGTGAAGGGTGAAACCGCAGAGCGGT                                                                        <  1:3087592‑M2/100‑1 (MQ=255)   cTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGT                                                                        <  2:1945649‑M2/100‑1 (MQ=255)   cTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGT                                                                        <  2:4013156‑M2/100‑1 (MQ=255)   cTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGT                                                                        >  2:1992890‑M2/1‑100 (MQ=255)   cTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGT                                                                        <  2:1215157‑M2/100‑1 (MQ=255)   cTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGT                                                                        >  1:1535260‑M2/1‑100 (MQ=255)   cTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGG                                                                         <  2:242265‑M2/99‑1 (MQ=255)   cTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGGGAAGGGTGAAACCGCAGAGCGG                                                                         <  2:4402675‑M2/99‑1 (MQ=255)   cTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCCCACGCACCGCACCGTGAAGGGTGAAACCGCACAGCGct                                                                        <  1:1019188‑M2/100‑3 (MQ=255)    ttATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTG                                                                       >  2:2435500‑M2/1‑100 (MQ=255)    ttATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTG                                                                       >  2:2585210‑M2/1‑100 (MQ=255)    ttATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTG                                                                       >  2:2953227‑M2/1‑100 (MQ=255)    ttATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTG                                                                       <  2:3081992‑M2/100‑1 (MQ=255)    ttATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTG                                                                       >  2:3511795‑M2/1‑100 (MQ=255)    ttATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCGCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTG                                                                       >  1:296738‑M2/1‑100 (MQ=255)     tATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAACGGTGAAACCGCAGAGCGGTG                                                                       >  2:4139527‑M2/1‑99 (MQ=255)                        ggTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTGcgtgaggcggtctgcttcgt                                                   >  1:1297288‑M2/1‑80 (MQ=255)                                                              ccTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTGcggtgcgtgaggcggtctgcttcgtcatagtgccagtgggtgaccagaccctcatcct             <  1:1177588‑M2/100‑59 (MQ=255)                                                                    gcaccgcaccGTGAAGGGTGAAACCGCAGAGCGGTGcggtgcgtgaggcggtctgcttcgtcatagtgccagtgggtgaccagaccctcatcctcgctgc       <  2:1796500‑M2/100‑65 (MQ=255)                                                                       ccgcaccgTGAAGGGTGAAACCGCAGAGCGGTGctgtgcgtgaggcggtctgcttcgtcatagtgagagtgggtgacaagaccgtcgtactcgctgcg      <  1:117169‑M2/98‑66 (MQ=255)                                                                         gcaccgTGAAGGGTGAAACCGCAGAGCGGTGcggtgcgtgaggcggtctgcttcgtcatagtgccagtgggtgaccagaccctcatcctcgctgcggata  <  2:2646685‑M2/100‑70 (MQ=255)                                                                         gcaccgTGAAGGGTGAAACCGCAGAGCGGTGcggtgcgtgaggcggtctgcttcgtcatagtgccagtgggtgaccagaccctcatcctcgctgcggata  <  1:481074‑M2/100‑70 (MQ=255)                                                                         |                                                                                                    AACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTGGCAGTATGACGAACGTGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGGCGGTGCA  >  NC_000913/3621242‑3621413

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.