New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3620220	NA (NA)	5 (0.060)	5/264	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

AGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCC  >  NC_000913/3620102‑3620368                                                                                                                       |                                                                                                                                                     agaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCC                                                                                                                       >  2:1176015‑M2/1‑150 (MQ=255)  gaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGAGGCGCACCGTCACACGGGCCG                                                                                                                      <  1:746147‑M2/150‑1 (MQ=255)   aaTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGG                                                                                                                     >  1:1483982‑M2/1‑150 (MQ=255)   aaTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGG                                                                                                                     >  1:347332‑M2/1‑150 (MQ=255)    aTTTACCTGCCGTGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                                                                    >  2:200600‑M2/1‑150 (MQ=255)    aTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                                                                    <  1:1279250‑M2/150‑1 (MQ=255)     tttACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                                                                   <  2:104400‑M2/150‑1 (MQ=255)     tttACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                                                                   >  2:290161‑M2/1‑150 (MQ=255)                                                                           aGTGGCAAACAGGTGCGCAGCTTTACTTAAGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCCCACGGGCCGGCCacggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcg                                             <  2:47632‑M2/150‑71 (MQ=255)                                                                                                tttACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACCCGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatag                        <  2:682440‑M2/150‑92 (MQ=255)                                                                                                               aaaTACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcg         <  2:1098935‑M2/150‑107 (MQ=255)                                                                                                                   aCCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggcag     >  1:1098935‑M2/1‑40 (MQ=255)                                                                                                                      ggggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggcaggta  <  2:1282178‑M2/150‑114 (MQ=255)                                                                                                                       |                                                                                                                                                     AGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCC  >  NC_000913/3620102‑3620368

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