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

ACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAG  >  NC_000913/3620204‑3620397                                                    |                                                                                                                                               acGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccag                                                                              >  1:276811‑M2/1‑52 (MQ=255) acGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccag                                                                              <  2:276811‑M2/118‑67 (MQ=255)              ccGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcg                                        <  1:809781‑M2/143‑105 (MQ=255)                                          ccccccccgccccggtatctctcgtcgtaagtaaagctgcgcccctgtttgccactccggtcatacaccactgccagttcgccgcgtggcgtccagccatagcgcacctgcggcgcggcaggtaaattctccgggtattccgg            >  2:704690‑M2/6‑11 (MQ=255)                                               ccggccccggtatttatcatagtaagtaaagctgcgcacctggttgacactacggtcatacaccaccgccagttcgccgcgtggcgtccaggcatagcgaaccagcggcgcgggaggtacatcccccgggtattccggcacgt       <  2:211801‑M2/143‑138 (MQ=255)                                                cggccccggtatttatcatcgtaagtaacgctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatcgcgcaccagcggcgcggcaggtaaattctccgggtattccgggtcgtg      >  2:551052‑M2/1‑5 (MQ=255)                                                 ggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggcaggtaaattctccgggtattccgggtcgtgc     >  2:644140‑M2/1‑4 (MQ=255)                                                  gccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggcaggtaaattctcgggggattccgggtcgtgcg    <  1:74984‑M2/143‑141 (MQ=255)                                                    cccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggcaggtaaattctccgggtattccgcgtcgtgcctc  <  1:156990‑M2/143‑143 (MQ=255)                                                    |                                                                                                                                               ACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAG  >  NC_000913/3620204‑3620397

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.