New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4866957	59 (1.200)	4 (0.090)	4/244	NT	6.4%	coding (154/1455 nt)	xdhA	xanthine dehydrogenase subunit XdhA
	?	NC_002947	= 4866973	60 (1.280)					coding (170/1455 nt)	xdhA	xanthine dehydrogenase subunit XdhA

GATCCAGTTTCTCGTCAACCAGGAGCTGCGTAGTGAGCACGCCCTGGACCCTAACATGACGGTGCTGCAGTACCTGCGCGAGCACCTGGGCAAAACCGGCACCAAGGAGGGCTGCGCCAGTGGAGACTGCGGCGCTTGCACCGTGGTGGTCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/4866806‑4866957 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggtcgTCCTGGATCAGTTCGCCGACCACCACGGT  <  NC_002947/4866973‑4866945                                                                                                                                                                                       GATCCAGTTTCTCGTCAACCAGGAGCTGCGTAGTGAGCACGCCCTGGACCCTAACATGACGGTGCTGCAGTACCTGCGCGAGCACCTGGGCAAAACCGGCACCAAGGAGGGCTGCGCCAGTGGAGACTGCGGCGCTTGCACCGTGGTGG                                  >  6:228849/1‑149    CCAGTTTCTCGTCAACCAGGAGCTGCGTAGTGAGCACGCCCTGGACCCTAACATGACGGTGCTGCAGTACCTGCGCGAGCACCTGGGCAAAACCGGCACCAAGGAGGGCTGCGCCAGTGGAGACTGCGGCGCTTGCACCGTGGTGGTCG                               >  3:105782/1‑149     CAGTTTCTCGTCAACCAGGAGCTGCGTAGTGAGCACGCCCTGGACCCTAACATGACGGTGCTGCAGTACCTGCGCGAGCACCTGGGCAAAACCGGCACCAAGGAGGGCTGCGCCAGTGGAGACTGCGGCGCTTGCACCGTGGTGGTCG                               <  2:84304/148‑1                                                 CCCTAACATGACGGTGCTGCAGTACCTGCGCGAGCACCTGGGCAAAACCGGCACCAAGGAGGGCTGCGCCAGTGGAGACTGCGGCGCTTGCACCGTGGTGGTCGTCCTGGATCAGTTCGCCGACCACCACGG   >  5:166043/1‑132                                                 CCCTAACATGACGGTGCTGCAGTACCTGCGCGAGCACCTGGGCAAAACCGGCACCAAGGAGGGCTGCGCCAGTGGAGACTGCGGCGCTTGCACCGTGGTGGTCGTCCTGGATCAGTTCGCCGACCACCACGG   <  6:166043/132‑1                                                                                                 CGGCACCAAGGAGGGCTGCGCCAGTGGAGACTGCGGCGCTTACACCGTGGTGGTCGTCCTGGATCAGTTCGCCGACCACCACGGT  >  3:132110/1‑85                                                                                                 CGGCACCAAGGAGGGCTGCGCCAGTGGAGACTGCGGCGCTTACACCGTGGTGGTCGTCCTGGATCAGTTCGCCGACCACCACGGT  <  4:132110/85‑1                                                                                                                                                                                       GATCCAGTTTCTCGTCAACCAGGAGCTGCGTAGTGAGCACGCCCTGGACCCTAACATGACGGTGCTGCAGTACCTGCGCGAGCACCTGGGCAAAACCGGCACCAAGGAGGGCTGCGCCAGTGGAGACTGCGGCGCTTGCACCGTGGTGGTCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/4866806‑4866957 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggtcgTCCTGGATCAGTTCGCCGACCACCACGGT  <  NC_002947/4866973‑4866945

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