New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AP009048	525265 =	NA (NA)	3 (0.050)	3/100	NT	100%	coding (2781/4281 nt)	rhsD	rhsD element protein
	?	AP009048	525295 =	0 (0.000)					coding (2811/4281 nt)	rhsD	rhsD element protein

GTGGCGGTAGACATAGTGCGCATCCTCCGCGATGCGGTTATCCGGCCACACTGTGAGTGTACTGTCCGGGt‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/525334‑525265 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TAACCGCATCGCGGAGGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGG  >  AP009048/525295‑525365                                                                                                                                               GTGGCGGTAGACATAGTGCGCATCCTCCGCGATGCGGTTATCCGGCCACACTGTGAGTGTACTGTCCGGGT                                                                        >  1:3110774/1‑71 GTGGCGGTAGACATAGTGCGCATCCTCCGCGATGCGGTTATCCGGCCACACTGTGAGTGTACTGTCCGGGT                                                                        <  1:4648658/71‑1 GTGGCGGTAGACATAGTGCGCATCCTCCGCGATGCGGTTATCCGGCCACACTGTGAGTGTACTGTCCGGGT                                                                        <  1:4944700/71‑1                                            GGCCACACTGTGAGTGTACTGTCCGGGTAACCGCATCGCGGAGGATGCGCACTATGTCTACCGCCACGAT                              >  1:5351511/1‑70                                            GGCCACACTGTGAGTGTACTGTCCGGGTAACCGCATCGCGGAGGATG                                                     <  1:2045296/47‑1                                                      TGAGTGTACTGTCCGGATAACCGCATCGCGGAGGATGCGCACTATGTCTACCGCCACGATGAATACGGCAG                   >  1:3932813/1‑71                                                                       TAACCGCATCGCGGAGGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGG  <  1:4467459/71‑1                                                                       TAACCGCACCGCGGAGGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGG  <  1:343836/71‑1                                                                                                                                               GTGGCGGTAGACATAGTGCGCATCCTCCGCGATGCGGTTATCCGGCCACACTGTGAGTGTACTGTCCGGGt‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/525334‑525265 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TAACCGCATCGCGGAGGATGCGCACTATGTCTACCGCCACGATGAATACGGCAGGCTGACGGAGAAGACGG  >  AP009048/525295‑525365

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