New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3650664 =	61 (1.020)	4 (0.070)	4/230	NT	6.8%	coding (202/609 nt)	PP_3216	cysteine dioxygenase
	?	NC_002947	3650697 =	52 (0.910)					coding (169/609 nt)	PP_3216	cysteine dioxygenase

CCCGACCCGCAGCGCTATCAGCAATACCTGCTGCACTGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/3650703‑3650664 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgcgGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGC  >  NC_002947/3650697‑3650834                                                                                                                                                                                    CCCGACCCGCAGCGCTATCAGCAATACCTGCTGCACTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGC                              >  2:908810/1‑150       CCGCAGCGCTATCAGCAATACCTGCTGCACTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTT                                      <  1:415479/136‑1       CCGCAGCGCTATCAGCAATACCTGCTGCACTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTT                                      >  2:415479/1‑136                             TGCTGCACTGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGC  <  1:908810/150‑1                                     TGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTC             >  1:182049/1‑131                                     TGCGGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTC             <  2:182049/131‑1                                                                                                                                                                                    CCCGACCCGCAGCGCTATCAGCAATACCTGCTGCACTGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/3650703‑3650664 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgcgGGTCGGGCTGGGCGAGGTCGTTCGGCAGCCAGTCGTCATGGGCGACCAGATCGGCCAGCAGGCCTTGGCCTGCGTCGAGAAGGGTGGCTTCGTCAGGCTGTTGGTCAAGCAGCGTGGCGAGGTGGTCGATGAAGTGGC  >  NC_002947/3650697‑3650834

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