New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3620220	NA (NA)	3 (0.080)	3/252	NT	100%	coding (1029/4236 nt)	rhsB	rhs element protein RhsB
	?	NC_000913	= 3620255	0 (0.000)					coding (1064/4236 nt)	rhsB	rhs element protein RhsB

TACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCAT  >  NC_000913/3620108‑3620357                                                                                                                 |                                                                                                                                          tACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGCACTGGCGGTGGTGTATGACCGTAGTGGCGAACAGGTGCGCAGCTTTACTTACGATGATAAAGACCGGGGCCGGATGGTGGCGCACCGTCCCGCGGGC                                                                                                             <  2:236190‑M2/143‑1 (MQ=255) tACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGC                                                                                                             >  1:708341‑M2/1‑143 (MQ=255)     tGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAGCTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGCGCCTGATGGTGGCGCACCGTCACACGGGCCGGC                                                                                                         >  2:748620‑M2/1‑143 (MQ=255)     tGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                                                         >  2:686976‑M2/1‑143 (MQ=255)      gccgcggcgcTGGGGCGATATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACTATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                                                        >  2:100980‑M2/1‑143 (MQ=255)      gccgcgccgcTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTGTGACCGTAGTGGCAAACAGGTGCGCAGCTTTCCTTACGATGATAGGTACCGGGGCCGGATGGCGGCGCACCGTCGCACGGGCCGGCC                                                                                                        <  2:202172‑M2/143‑1 (MQ=255)      gccgcgccgcTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                                                        <  1:476623‑M2/143‑1 (MQ=255)      gccgcgccgcGGGGGCGCTATGGCTGGACGCCACGCGGCGAACGGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                                                        >  1:557696‑M2/1‑143 (MQ=255)         gcgccgcTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGGCCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGAAGAAAAATACCGGGGCCGGAAGGCGGCGCACCGACACAGTGGCCGGacaga                                                                                                     <  2:128952‑M2/143‑6 (MQ=255)                                                                                                        tAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagctcgccgcgtggcgtccagccagagcgcaccc      <  2:75526‑M2/143‑99 (MQ=255)                                                                                                                 gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcgg  <  1:364281‑M2/138‑103 (MQ=255)                                                                                                                 gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcgg  >  2:364281‑M2/1‑36 (MQ=255)                                                                                                                 |                                                                                                                                          TACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCAT  >  NC_000913/3620108‑3620357

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