New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 66118	68 (0.930)	6 (0.090)	6/232	NT	8.4%	coding (1176/1653 nt)	betA‑I	choline dehydrogenase
	?	NC_002947	= 66126	67 (0.980)					coding (1184/1653 nt)	betA‑I	choline dehydrogenase

ATGGCGCCGAGCCAGTTAGGCGCGTTCGTGCGTTCAAGCCCCGAACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAACGCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/65963‑66118 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcgcccgCGGCTGGCCGGGCGCAGG  <  NC_002947/66126‑66109                                                                                                                                                                                ATGGCGCCGAGCCAGTTAGGCGCGTTCGTGCGTTCAAGCCCCGAACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAACGCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTG                          <  2:174835/150‑1     CGCCGAGCCAGTTAGGCGCGTTCGTGCGTTCAAGCCCCGAACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAACGCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCC                      <  2:293237/150‑1                      GCGTTCGTGCGTTCAAGCCCCGAACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAACGCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGCGGCTGGCCGG         <  1:1835893/146‑1                      GCGTTCGTGCGTTCAAGCCCCGAACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAACGCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGCGGCTGGCCGG         >  2:1835893/1‑146                                                                                CGCTGTCACTGGAACGCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGC                       >  1:82359/1‑74                                                                                CGCTGTCACTGGAACGCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGC                       <  2:82359/74‑1                                                                                          GGAACGCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGCGGCTGGCCGGGCGCAG   >  1:255336/1‑84                                                                                          GGAACGCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGCGGCTGGCCGGGCGCAG   <  2:255336/84‑1                                                                                                   CGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGCGGCTGGCCGGGCGCAGG  >  1:1894629/1‑76                                                                                                   CGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCGCGGCTGGCCGGGCGCAGG  <  2:1894629/76‑1                                                                                                                                                                                ATGGCGCCGAGCCAGTTAGGCGCGTTCGTGCGTTCAAGCCCCGAACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAACGCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGTGCAACCTGCGCCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/65963‑66118 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcgcccgCGGCTGGCCGGGCGCAGG  <  NC_002947/66126‑66109

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