New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3620569 =	NA (NA)	39 (0.280)	10/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‑GCCTGACGGGCTGGAATTGCGCCGGGA  >  NC_000913/3620531‑3620654                                       |                                                                                       aTGTGGTGCCGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCC                          >  1:2220597‑M2/1‑100 (MQ=255) aTGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCC                          >  1:1837470‑M2/1‑100 (MQ=255) aTGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCC                          >  1:1837723‑M2/1‑100 (MQ=255) aTGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCC                          >  1:2019041‑M2/1‑100 (MQ=255)  tgtgGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCT                         <  2:914192‑M2/100‑1 (MQ=255)   gtggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTG                        <  1:1578399‑M2/100‑1 (MQ=255)   gtggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTG                        >  2:3872521‑M2/1‑100 (MQ=255)    tggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGA                       <  1:1367505‑M2/100‑1 (MQ=255)    tggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGA                       <  2:2088717‑M2/100‑1 (MQ=255)     ggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGAC                      <  2:964365‑M2/100‑1 (MQ=255)     ggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGAC                      <  2:1115716‑M2/100‑1 (MQ=255)     ggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGAC                      <  1:4027671‑M2/100‑1 (MQ=255)     ggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGAC                      >  1:277996‑M2/1‑100 (MQ=255)      gtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGACG                     <  2:2701149‑M2/100‑1 (MQ=255)      gtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGACG                     >  1:709482‑M2/1‑100 (MQ=255)           gggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGCCCGGT‑GGCTGACGttaac                >  2:1610370‑M2/1‑95 (MQ=255)           gggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaac                >  1:1261210‑M2/1‑95 (MQ=255)           gggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaac                >  1:2633745‑M2/1‑95 (MQ=255)           gggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCAGCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGCCCGGT‑GGCTGACGttaac                >  1:1912140‑M2/1‑95 (MQ=255)            ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaact               >  1:3350471‑M2/1‑94 (MQ=255)            ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaact               >  1:790495‑M2/1‑94 (MQ=255)            ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaact               >  1:543917‑M2/1‑94 (MQ=255)            ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaact               >  1:3049687‑M2/1‑94 (MQ=255)            ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaact               >  1:2078136‑M2/1‑94 (MQ=255)            ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTGTTACTATAACCACCACAACCAGTTAACGTCAGGCCCGGT‑GGCTGACGttaact               >  2:4037950‑M2/1‑94 (MQ=255)             gCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGC‑CCGGTGGCTGACGttaactg              >  1:1327573‑M2/1‑93 (MQ=255)                 catcaCGCGCATAACCACGCAGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGCCCGGT‑GGCTGACGttaactggttg          >  1:3294327‑M2/1‑89 (MQ=255)                         gcATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGCCCGGT‑GGCTGACGttaactggttgtggtggtt  >  2:2145512‑M2/1‑81 (MQ=255)                          cATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCCGGTACCGTCAGCCACCGG‑GCCTCACGttcactggttgtggtggt   <  1:604037‑M2/98‑19 (MQ=255)                                       |                                                                                       ATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGG‑GCCTGACGGGCTGGAATTGCGCCGGGA  >  NC_000913/3620531‑3620654

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