New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3620116	NA (NA)	4 (0.040)	4/262	NT	100%	coding (925/4236 nt)	rhsB	rhs element protein RhsB
	?	NC_000913	= 3620152	0 (0.000)					coding (961/4236 nt)	rhsB	rhs element protein RhsB

CCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAG  >  NC_000913/3620003‑3620301                                                                                                                                                      |                                                                                                                                                      ccGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCG                                                                                                                                                       >  1:512362‑M2/1‑150 (MQ=255) ccGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTATGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGACGCTGGTGCGATATGGCTGGACGCCACGCGGCG                                                                                                                                                       >  2:347410‑M2/1‑150 (MQ=255)                                                                                                       aTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacagacggatgccgttgtcccggccatattcggtgtaacccggcagggtat                                                 <  1:1588775‑M2/150‑103 (MQ=255)                                                                                                            ccTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacagacggatgccgttgtcccggccatattcggtgtaacccggcagggtatcagga                                            <  1:840592‑M2/150‑108 (MQ=255)                                                                                                               gccgcgccgcTGGTGCGCTATGGCTGGACGCCACGCGGCGcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacagacggatgccgttgtcccggccatattcggtgtaacccggcagggtatcaggaaaa                                         <  2:1331788‑M2/150‑111 (MQ=255)                                                                                                                          ggTGCGCTATGGCTGGACGCCACGCGGCGcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacagacggatgccgttgtcccggccatattcggtgtaacccggcagggtatcaggaaaagcggacggttc                              >  2:840592‑M2/1‑29 (MQ=255)                                                                                                                                                cgcggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacagacggatgccgttgtcccggccatattcggtgtaacccggcagggtatcaggaaaagcggacggttccgtcccgccggaaatggcctgc        >  1:1331788‑M2/1‑7 (MQ=255)                                                                                                                                                 gcggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacagacggatgccgttgtcccggccatattcggtgtaacccggcagggtatcaggaaaagcggacggttccgtcccgccggaaatggcctgct       <  1:962304‑M2/150‑145 (MQ=255)                                                                                                                                                 gcggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacagacggatgccgttgtcccggccatattcggtgtaacccggcagggtatcaggaaaagcggacggttccgtcccgccggaaatggcctgct       <  2:1423209‑M2/150‑145 (MQ=255)                                                                                                                                                   ggcgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacagacggatgccgttgtcccggccatattcggtgtaacccggcagggtatcaggaaaagcggacggttccgtcccgccggaaatggcctgctgc     <  1:1333855‑M2/150‑147 (MQ=255)                                                                                                                                                    aagcgtcatgtaaattatcagggtattccgggtagtgcgtcaacaacaccgcagacagacggatgcagttgtcccggccatattaggtgtaacccggcagggtatcaggaaaagcggacggttccgtcccgccggaaatggcctgctgcc    <  2:805481‑M2/148‑148 (MQ=255)                                                                                                                                                     cgcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacagacggatgccgttgtcccggccatattcggtgtaacccggcagggtatcaggaaaagcggacggttccgtcccgccggaaatggcctgctgccg   >  2:1491334‑M2/1‑2 (MQ=255)                                                                                                                                                      gcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacagacggatgccgttgtcccggccatattcggtgtaacccggcagggtatcaggaaaagcggacggttccgtcccgccggaaatggcctgctgccgg  >  2:1093414‑M2/1‑1 (MQ=255)                                                                                                                                                      gcggcaggtaaattctccgggtattccgggtcgtgcgtcagccacacggcagacagacggatgccgttgtcccggccatattcggtgtaacccggcagggtatcaggaaaagcggacggttccgtcccgccggaaatggcctgctgccgg  >  2:956925‑M2/1‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAG  >  NC_000913/3620003‑3620301

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