New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3620569 =	NA (NA)	35 (0.240)	11/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

ATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGACGGGCTGGAATTGC  >  NC_000913/3620531‑3620647                                       |                                                                                aTGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGTATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCC                   >  1:1748687‑M2/1‑100 (MQ=255) aTGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCC                   >  1:1516303‑M2/1‑100 (MQ=255) aTGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCC                   >  2:1138120‑M2/1‑100 (MQ=255) aTGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCC                   >  2:602672‑M2/1‑100 (MQ=255) aTGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCC                   >  1:2762069‑M2/1‑100 (MQ=255)   gtggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTG                 >  1:607064‑M2/1‑100 (MQ=255)   gtggtgACGGGCCTCATCACGCGCATAACCACGACGGATGGCAGGGCATCGGCGTTTTACTATAACCACAACAACCAGTTAACGTCAGCCACCGG‑GCCT                  <  1:2573955‑M2/99‑1 (MQ=255)    tggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGA                <  2:1174536‑M2/100‑1 (MQ=255)     ggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGAC               <  1:2594093‑M2/100‑1 (MQ=255)     ggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGAC               <  2:345626‑M2/100‑1 (MQ=255)     ggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGAC               <  2:3012350‑M2/100‑1 (MQ=255)     ggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGAC               >  1:3086935‑M2/1‑100 (MQ=255)     ggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGAC               <  2:3555439‑M2/100‑1 (MQ=255)     ggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGAC               >  2:76038‑M2/1‑100 (MQ=255)     ggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCC                   <  1:2597433‑M2/96‑1 (MQ=255)      gtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGACG              >  1:1613521‑M2/1‑100 (MQ=255)      gtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGACG              >  1:1677165‑M2/1‑100 (MQ=255)      gtgACGGGCCTCATCACGCGCATAACCACGCCGGATGCCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGACG              <  2:1647580‑M2/100‑1 (MQ=255)        gACCGGCCTCATCACGCGCATAACCACGCCTGATGGCATGGCATCGGCGTTTTACTATCACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGACG              >  1:3588164‑M2/1‑98 (MQ=255)           gggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCAGAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaac         >  2:3011056‑M2/1‑95 (MQ=255)           gggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaac         >  2:1059505‑M2/1‑95 (MQ=255)           gggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaac         >  1:3023326‑M2/1‑95 (MQ=255)            ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGCCCGGT‑GGCTGACGttaact        >  2:2202957‑M2/1‑94 (MQ=255)            ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaact        >  1:1576312‑M2/1‑94 (MQ=255)            ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaact        >  2:328319‑M2/1‑94 (MQ=255)            ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaact        >  2:3121873‑M2/1‑94 (MQ=255)            ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaact        >  1:3457300‑M2/1‑94 (MQ=255)             gCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaactg       >  2:3305827‑M2/1‑93 (MQ=255)             gCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaactg       >  1:2047678‑M2/1‑93 (MQ=255)              ccTCATCACGCGCATACCCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaact        >  2:1745059‑M2/1‑92 (MQ=255)                  atcaCGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaactggttgt  >  2:946121‑M2/1‑88 (MQ=255)                                       |                                                                                ATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGACGGGCTGGAATTGC  >  NC_000913/3620531‑3620647

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.