New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3620350	NA (NA)	3 (0.040)	3/264	NT	100%	coding (1159/4236 nt)	rhsB	rhs element protein RhsB
	?	NC_000913	= 3620425	0 (0.000)					coding (1234/4236 nt)	rhsB	rhs element protein RhsB

TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAGACGGA  >  NC_000913/3620272‑3620498                                                                               |                                                                                                                                                     tACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTG                                                                               <  1:802040‑M2/150‑1 (MQ=255)  aCGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                                                                              >  1:1057250‑M2/1‑150 (MQ=255)  aCGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                                                                              >  1:404311‑M2/1‑150 (MQ=255)  aCGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                                                                              <  2:1073102‑M2/150‑1 (MQ=255)   cGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                             <  1:115489‑M2/150‑1 (MQ=255)   cGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                             <  1:1220800‑M2/150‑1 (MQ=255)   cGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                             <  1:277847‑M2/150‑1 (MQ=255)   cGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                             <  1:682068‑M2/150‑1 (MQ=255)   cGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                                                                             >  2:1087740‑M2/1‑150 (MQ=255)     aCAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGA                                                                           >  2:189387‑M2/1‑150 (MQ=255)                                                                        aGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagctgttctgtcacccgccc        <  2:168361‑M2/150‑68 (MQ=255)                                                                           ccGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagctgttctgtcacccgcccgtc     <  2:1143672‑M2/150‑71 (MQ=255)                                                                              catcaccatcacCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatacgtgtagcttaagcctgccgggtttagctgttctgtcacccgcccgtcgct  <  2:151145‑M2/150‑74 (MQ=255)                                                                               |                                                                                                                                                     TACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGCGTCACGCAGAGTCAGTTTGACGCCGTGGGCAGGCTCAGGGCACAGACGGA  >  NC_000913/3620272‑3620498

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