New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 386412	35 (1.320)	3 (0.110)	3/276	NT	7.9%	coding (101/1296 nt)	phoR	sensory histidine kinase in two‑component regulatory system with PhoB
	?	REL606	= 666130	NA (NA)					noncoding (1443/1443 nt)	IS150	repeat region

GGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCTTCCCGGCTTTCATCCTGGGTGCATTTTTTGGTTACCTGCCCTGGTTt‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/386268‑386412 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGTACTGACCCCAAAAAGTTGGACAGTTAAACACGAGGCATATAGGTCTGATTCCGATATTCAATTGGAGTCAGACCTTTTAATTTCAGGCTAATTCTTCTGCTGTTGTAGTATTCAATAT  <  REL606/666130‑666010                                                                                                                                                                                                                                                                            GGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCTTCCCGGCTTTCATCCTGGGTGCATTTTTTGGTTACCTGCCCTGGTTTGTAC                                                                                                                      <  2:68714/150‑1                                      ATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCTTCCCGGCTTTCATCCTGGGTGCATTTTTTGGTTACCTGCCCTGGTTTGTACTGACCCCAAAAAGTTGGACAGTTAAACACGAGGCATA                                                                                 >  2:143758/1‑150                                          TACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCTTCCCGGCTTTCATCCTGGGTGCATTTTTTGGTTACCTGCCCTGGTTTGTACTGACCCCAAAAAGTTGGACAGTTAAACACGAGGCATATAGG                                                                             <  1:285080/150‑1                                                                                                                     GGTTCATTTTTTGGTTACCTGCCCTGGTTTGTACTGACCCCAAAAAGTTTGACAGTTAGACACGAGGCATATAGGTCTGATTCCGATATTCAATTGGAGTCAGACCTTTTAATTTCAGGCTAATTCTTCTGCTGTTGTAGTATTCAATAT  <  1:143347/150‑1                                                                                                                                                                                                                                                                            GGACGCAGAGCAGGGCTTATGATTTCTTAACTGGAGTATCTTACGTGCTGGAACGGCTGTCGTGGAAAAGGCTGGTGCTGGAGCTGCTACTTTGCTGCTTCCCGGCTTTCATCCTGGGTGCATTTTTTGGTTACCTGCCCTGGTTt‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/386268‑386412 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGTACTGACCCCAAAAAGTTGGACAGTTAAACACGAGGCATATAGGTCTGATTCCGATATTCAATTGGAGTCAGACCTTTTAATTTCAGGCTAATTCTTCTGCTGTTGTAGTATTCAATAT  <  REL606/666130‑666010

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