New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	734797 =	48 (0.550)	9 (0.110)	3/168	NT	28.4%	pseudogene (665/1521 nt)	rhsO	pseudogene, Rhs family
	?	NC_000913	3622561 =	0 (0.000)					coding (3370/4236 nt)	rhsB	Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

CCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCG  >  NC_000913/3622464‑3622659                                                                                                  |                                                                                                   tccagccggtcgagcatctgcaccgtgccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTg                                                                                                  >  1:2254133‑M2/34‑100 (MQ=255) ccagccggtcgagcatctgcaccagcaccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTg                                                                                                  >  1:377141‑M2/34‑100 (MQ=255)  ccagccggtcgagcatctgcaccgtgccggctGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGa                                                                                                 >  1:218095‑M2/33‑100 (MQ=255)                                  nGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCg                                                                 >  2:1545099‑M2/2‑100 (MQ=255)                                  GGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCGGCCACTGACTGGCATCGTGCGGCGTGACCTTGGAGCAGATGCAAAACCAGATGGACCCg                                                                 <  2:748487‑M2/100‑1 (MQ=255)                                  GGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTTACCGTGGAGCAGATGCAAAACCAGATGGACCCg                                                                 <  2:1932592‑M2/100‑1 (MQ=255)                                  GGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGAccc                                                                  >  2:424625‑M2/1‑99 (MQ=255)                                  GGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCg                                                                 >  1:1296212‑M2/1‑100 (MQ=255)                                  GGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCg                                                                 >  2:457824‑M2/1‑100 (MQ=255)                                   GAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCgg                                                                <  1:618226‑M2/100‑1 (MQ=255)                                   GAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCgg                                                                <  1:580520‑M2/100‑1 (MQ=255)                                   GAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGATGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCgg                                                                >  1:999769‑M2/1‑100 (MQ=255)                                    AAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGgt                                                               <  1:562392‑M2/100‑1 (MQ=255)                                    AAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGgt                                                               <  1:2084775‑M2/100‑1 (MQ=255)                                     AAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGgtg                                                              <  2:2027335‑M2/100‑1 (MQ=255)                                     AAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGgtg                                                              >  1:1010423‑M2/1‑100 (MQ=255)                                     AAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGgtg                                                              >  1:1064314‑M2/1‑100 (MQ=255)                                                                 ccgtcaggccgcacgatgccagccagcggcggcCTGACGGTGGAGCAGATGCAAAACCAGATGGCCCCGGTGTACACGCCGGCGCGAAAAATCCACCTGt                                  >  2:107457‑M2/34‑100 (MQ=255)                                                                      gtcaggccgcacgatgccagccagcggcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCAc                             >  1:1248988‑M2/29‑100 (MQ=255)                                                                      aggccgcacgatgccagccagcggcggcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCAc                             >  2:911671‑M2/29‑100 (MQ=255)                                                                       ggccgcacgatgccagccagcggcggcCTGACCGTGGGGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt                            >  1:2119789‑M2/28‑100 (MQ=255)                                                                       ggccgcacgatgccagccagcggcggcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt                            >  1:635367‑M2/28‑100 (MQ=255)                                                                       ggccgcacgatgccagccagcggcggcCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt                            >  2:2213734‑M2/28‑100 (MQ=255)                                                                       ggccgcacgatgccagccagcggcggcCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCGGTACCACt                            >  2:2213757‑M2/28‑100 (MQ=255)                                                                           gccgcacgatgccagccagcggcCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGa                        >  1:574896‑M2/24‑100 (MQ=255)                                                                           gcacgatgccagccagcggcggcCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGa                        >  2:2287527‑M2/24‑100 (MQ=255)                                                                                             gcggcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCt      <  1:783407‑M2/95‑1 (MQ=255)                                                                                              cggcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCt      >  2:608665‑M2/5‑99 (MQ=255)                                                                                              cggcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTg     <  2:1063028‑M2/96‑1 (MQ=255)                                                                                               ggcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTgg    <  1:1352384‑M2/97‑1 (MQ=255)                                                                                               ggcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTgg    <  1:665446‑M2/97‑1 (MQ=255)                                                                                                gcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgc   <  1:304610‑M2/98‑1 (MQ=255)                                                                                                gcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgc   <  1:2242054‑M2/98‑1 (MQ=255)                                                                                                gcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgc   <  2:1200850‑M2/98‑1 (MQ=255)                                                                                                gcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgc   >  1:697817‑M2/3‑100 (MQ=255)                                                                                                gcCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACAAGAAAGAGAGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgc   <  2:1844860‑M2/98‑1 (MQ=255)                                                                                                 cCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGgcg  >  1:2255233‑M2/2‑100 (MQ=255)                                                                                                  |                                                                                                   CCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCG  >  NC_000913/3622464‑3622659

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