New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3620862	NA (NA)	3 (0.020)	3/428	NT	2.6%	coding (1671/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
	?	NC_000913	= 3621013	114 (0.910)					coding (1822/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Rejected: Frequency below/above cutoff threshold.

CGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGA  >  NC_000913/3620773‑3621098                                                                                          |                                                                                                                                                                                                                                             cGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTG                                                                                          >  2:651801‑M2/1‑238 (MQ=255)                                                                            gaccatgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAc              <  2:525480/239‑1 (MQ=31)                                                                             accatgaccGCAGCGGGCAGATGACGGCGGTGGACCGCGAGTAAGGGCTGAGGCAGTACCGCGCATACGACAGCGGTGGACAGTTAATTGCCGAGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAc              <  2:564748/238‑1 (MQ=255)                                                                              ccatgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAcg             <  1:1224242/238‑1 (MQ=31)                                                                               tgggacctgtccGGGCGGAGGACGGCGTGGCCCCGCGGGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGCCAGTTAATTGCCGGGAAAGACACGCAGGGCCATGAAACGCGGTATGAATCCAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAcgcg           <  2:224567/228‑1 (MQ=11)                                                                               catgcccGCTGCGGGCAGATGGCGGAGGTGCACGGGGAGGAAGGGCGGAGTCAGTACGGGGCATACGACAGCCGTGGAAAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTAGGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAcgcg           <  2:408945/239‑1 (MQ=11)                                                                               catgaccGCTTCGGGCAGATGACGGCGGTGCACCGGGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAcgcg           <  1:1086201/239‑1 (MQ=14)                                                                               catgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCGGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAcgcg           <  2:976938/239‑1 (MQ=14)                                                                               catgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAcgcg           <  2:1337502/239‑1 (MQ=31)                                                                               catgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAcgc            >  1:866922/1‑238 (MQ=31)                                                                               catgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAcgc            >  2:706211/1‑238 (MQ=31)                                                                                atgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAcgcg           >  1:65561/1‑238 (MQ=31)                                                                                atgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCAGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACgcgc          >  2:1037354/1‑239 (MQ=14)                                                                                  gaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAg        >  1:421172/1‑239 (MQ=31)                                                                                  gaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAg        >  1:847697/1‑239 (MQ=31)                                                                                  gaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAg        >  1:598976/1‑239 (MQ=31)                                                                                   accGGTTCGGGAAGAGGCCGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACGGCGCATACGACAGCCGTGGCAAGTTAAGTCCCGTGAAAGACACGCGGGGCCATGAAACGCGGTATGAATACAACAGCGCGGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCa         <  2:696034/237‑1 (MQ=11)                                                                                   accGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCCGCAGAAACGGGACACAGTACGAAGCGGGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCCGt       >  2:1217062/1‑239 (MQ=9)                                                                                   accGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGGAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGt       >  1:701027/1‑239 (MQ=14)                                                                                    ccGTTTGGGGCAGAGGACGGGGGTGCACCGCGAGGAAGGGCTGAGTCAGTACGGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTa      <  2:721348/239‑1 (MQ=9)                                                                                    ccGCTTCGGGCAGAGGGCGGCGGTGCACCGCGAGGAGGGGCTGGGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGCCACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTa      <  2:455472/239‑1 (MQ=9)                                                                                       cTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGGCGCGCAGTATgg   >  2:193504/1‑239 (MQ=14)                                                                                       cTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTa      >  2:591732/1‑236 (MQ=31)                                                                                       cTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTa      <  1:591732/236‑1 (MQ=31)                                                                                       cTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATgg   >  1:748103/1‑239 (MQ=31)                                                                                        ttCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTAAGGa  >  2:246955/1‑239 (MQ=14)                                                                                          |                                                                                                                                                                                                                                             CGGAAAACCATGACGTGGAGCCGTTACGGTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGA  >  NC_000913/3620773‑3621098

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