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

ACCGTCCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATC  >  NC_000913/3619998‑3620262                                                                                                                       |                                                                                                                                                   aCCGTCCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACG                                                                                                                     <  2:1187082‑M2/150‑1 (MQ=255) aCCGTCCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACG                                                                                                                     >  2:1230796‑M2/1‑150 (MQ=255)  ccGTCCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGC                                                                                                                    <  1:1366561‑M2/150‑1 (MQ=255)   cGTCCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCG                                                                                                                   >  2:920903‑M2/1‑150 (MQ=255)    gTCCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGG                                                                                                                  >  1:777272‑M2/1‑150 (MQ=255)     tCCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGC                                                                                                                 >  2:1397247‑M2/1‑150 (MQ=255)      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)                                                                                                                       |                                                                                                                                                   ACCGTCCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATC  >  NC_000913/3619998‑3620262

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.