New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	523984 =	23 (0.810)	3 (0.120)	3/116	2.3	12.2%	coding (724/4281 nt)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	524076 =	23 (0.930)					coding (816/4281 nt)	rhsD	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Rejected: Frequency below/above cutoff threshold.

GCTGACGTACCGGCGTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACCTC  >  NC_000913/523917‑524045                                                                    |                                                              gcggatacccctgtcggggccgtattcggtaccgggcagtgtgtcggggaacgctgaggctgagagag                                                               <  2:652505‑M1/1‑1 (MQ=255) gcggatacccctgtcggggccgtattcggtaccgggcagtgtgtcggggaacgctgaggctgagagag                                                               <  2:466205‑M1/1‑1 (MQ=255)  cTGACGTACCGGCTTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCcgg                                                              <  2:1039577/68‑1 (MQ=255)   tGACGTACCGGCGTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCcggg                                                             <  2:139580/68‑1 (MQ=255)   ggatacccctgtcggggccgtattcggtaccgggcagtgtgtcggggaacgctgaggctgagagaggg                                                             >  1:451421‑M1/66‑68 (MQ=255)   ggatacccctgtcggggccgtattcggtaccgggcagtgtgtcggggaacgctgaggctgagagaggg                                                             >  1:466201‑M1/66‑68 (MQ=255)     aCGTACCGGCGTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCcgggcg                                                           >  2:979982/1‑68 (MQ=255)         aCCGGCGTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAg                                                       >  2:946443/1‑68 (MQ=255)           cGGCGTGAGGCCGCCGGTGACCTGGCCGGGGAAATCGCCGGCGTGACGGACGGTGCCGGGCGGGAGtt                                                     >  1:390269/1‑68 (MQ=255)              cGTGAGGCCGCCGGTGACCTGGCCGGGGGAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCg                                                  >  2:663296/1‑68 (MQ=255)               gTGAGGCCGCCGGTGACCTGGCCGGGGATATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGt                                                 <  2:413777/68‑1 (MQ=255)               gTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGt                                                 <  1:128063/68‑1 (MQ=255)                tGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTc                                                >  2:44466/1‑68 (MQ=255)                    gccgccGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGt                                            <  2:915181/68‑1 (MQ=255)                            tGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCa                                    >  2:976286/1‑68 (MQ=255)                               accgggcagtgtgtcggggaacgctgaggctgagagaGGGCGGGAGTTCCGTCTGGTGCTGACCACGc                                 <  1:197167‑M1/31‑1 (MQ=255)                                        ggATATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGc                        <  2:693389/68‑1 (MQ=255)                                         gagaTCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCAGACCACGCAGGCGCAGCg                       >  1:626008/4‑68 (MQ=255)                                         gAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCg                       >  2:792824/1‑68 (MQ=255)                                          tgtcggggaacgctgaggctgagagaGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGt                      >  2:1033064‑M1/27‑68 (MQ=255)                                               aCCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTgg                                             <  2:33323/40‑1 (MQ=255)                                               aCCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTgg                                             >  1:33323/1‑40 (MQ=255)                                               aCCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGa                 <  1:976278/68‑1 (MQ=255)                                                    cGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAgg            <  2:530625/68‑1 (MQ=255)                                                    cGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAgg            <  2:927498/68‑1 (MQ=255)                                                          ggacggTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCa      <  1:460268/68‑1 (MQ=255)                                                          ggacggTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCa      >  1:311569/1‑68 (MQ=255)                                                           gacggTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCAc     >  2:381436/1‑68 (MQ=255)                                                            acggTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCAcc    >  2:916571/1‑68 (MQ=255)                                                              gagagaGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACctc  >  2:197168‑M1/7‑68 (MQ=255)                                                                    |                                                              GCTGACGTACCGGCGTGAGGCCGCCGGTGACCTGGCCGGGGAAATCACCGGCGTGACGGACGGTGCCGGGCGGGAGTTCCGTCTGGTGCTGACCACGCAGGCGCAGCGTGCGGAAGAGGCCCGCACCTC  >  NC_000913/523917‑524045

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