New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	3038086 =	55 (1.050)	4 (0.080)	4/234	NT	7.4%	coding (315/1143 nt)	gbd	4‑hydroxybutyrate dehydrogenase
	?	NC_002947	3038112 =	48 (0.960)					coding (341/1143 nt)	gbd	4‑hydroxybutyrate dehydrogenase

GTATCAGGATGGCACCGCGCCCGACCTCGCTCCCCGTGCCGGCCGTGGTCGGCACCGCAATCACCGGCGCGGTTGCCGGGGTGATCCGGGCCAGGCCGCCCTCGATCACCGCAAAGCTCTTCAACGGGCCGT‑C‑GTGCCGGGCACACACGGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/3038237‑3038086 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acggccCGTTGAAGAGCTTTGCGGTGATCGAGGGCGGCCTGGCCCGGATC  >  NC_002947/3038112‑3038155                                                                                                                                                                                                        GTATCAGGATGGCACCGCGCCCGACCTCGCTCCCCGTGCCGGCCGTGGTCGGCACCGCAATCACCGGCGCGGTTGCCGGGGTGATCCGGGCCAGGCCGCCCTCGATCACCGCAAAGCTCTTCAACGGGCCGT‑C‑GTGCCGGGCACACACGG                                                >  1:809369/1‑150                                            CGTGGTCGGCACCGCAATCACCGGCGCGGTTGCCGGGGTGATCCGGGCCAGGCCGCCCTCGATCACCGCAAAGCTCTTCAACGGGCCGTGT‑GTGCCCGGCAC‑GACGGCCCGTTGAA                                       <  1:836751/116‑1                                            CGTGGTCGGCACCGCAATCACCGGCGCGGTTGCCGGGGTGATCCGGGCCAGGCCGCCCTCGATCACCGCAAAGCTCTTCAACGGGCCGT‑GTGTGCCCGGCAC‑GACGGCCCGTTGAA                                       >  2:836751/1‑116                                                                                                                                GCCGT‑C‑GTGCCGGGCACACACGGCCCGTTGAAGAGCTTTGCGGTGATCGAGGGCGGCCTGGCCCGGATC  <  1:207973/69‑1                                                                                                                                GCCGT‑C‑GTGCCGGGCACACACGGCCCGTTGAAGAGCTTTGCGGTGATCGAGGGCGGCCTGGCCCGGATC  >  2:207973/1‑69                                                                                                                                                        GCCCGTTGAAGAGCTTTGCGGTGATCGAGGGCGGCCTGGCCCGGAT   <  1:1395014/46‑1                                                                                                                                                        GCCCGTTGAAGAGCTTTGCGGTGATCGAGGGCGGCCTGGCCCGGAT   >  2:1395014/1‑46                                                                                                                                                                                                        GTATCAGGATGGCACCGCGCCCGACCTCGCTCCCCGTGCCGGCCGTGGTCGGCACCGCAATCACCGGCGCGGTTGCCGGGGTGATCCGGGCCAGGCCGCCCTCGATCACCGCAAAGCTCTTCAACGGGCCGT‑C‑GTGCCGGGCACACACGGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/3038237‑3038086 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acggccCGTTGAAGAGCTTTGCGGTGATCGAGGGCGGCCTGGCCCGGATC  >  NC_002947/3038112‑3038155

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