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

GCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCA  >  NC_000913/3620113‑3620396                                                                                                                                               |                                                                                                                                              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)                                                                                                   tAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagctcgccgcgtggcgtccagccagagcgcaccc                                             <  2:75526‑M2/143‑99 (MQ=255)                                                                                                            gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcgg                                         <  1:364281‑M2/138‑103 (MQ=255)                                                                                                            gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcgg                                         >  2:364281‑M2/1‑36 (MQ=255)                                                                                                                                            tggcccggtatttcccatcgtaagtacagcggcgcacctgtgtgccactacggtggttcgctgccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggcaggtaaattctccgggtattccgggtcgtgc    >  2:697954‑M2/4‑4 (MQ=255)                                                                                                                                            ggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatcgcgcaccagcggcgcggccggtacattctcccggtcttccgggtcctac    <  2:167531‑M2/143‑140 (MQ=255)                                                                                                                                              ccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccaggtcgccgcgtggcgtccagccatagcgcactagcggcgcgtcaggtaaattctccgggtattccgggtcgtgcgt  >  1:554215‑M2/1‑2 (MQ=255)                                                                                                                                               |                                                                                                                                              GCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCA  >  NC_000913/3620113‑3620396

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