New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 65250	77 (1.480)	4 (0.080)	4/258	NT	5.5%	coding (308/1653 nt)	betA‑I	choline dehydrogenase
	?	NC_002947	= 65282	64 (1.270)					coding (340/1653 nt)	betA‑I	choline dehydrogenase

CCCGTTGGTTACCTCTACTGCATCGGCAACCCGCGCACAGACTGGTGCTTCAAGACCGAGGCACAACCTGGCCTCGGCGGTCGCGCTTTGGGTTATCCAAGGGGCAAGGTGCTAGGTGGCTGTTCGTCGATCAACGGCATGATCTACATGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/65099‑65250 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcgCGGCCCAATGGTCGTAGTCGGCGGCTTGGCCCCGCATGTAGATCATGCCGTTGATCGACGAACAGCCACCTAGCACCTTGCCCCTTGGATAACCCAAAGCGCGACCGCCGAGGCCAGGTTGTGCCTCGGTCTTGAAGCACCAGTC  <  NC_002947/65282‑65138                                                                                                                                                                                                                                                                                                           CCCGTTGGTTACCTCTACTGCATCGGCAACCCGCGCACAGACTGGTGCTTCAAGACCGAGGCACAACCTGGCCTCGGCGGTCGCGCTTTGGGTTATCCAAGGGGCAAGGTGCTAGGTGGCTGTTCGTCGATCAACGGCATGATCTACAT                                                                                                                                                      <  5:99556/149‑1  CCGTTGGTTACCTCTACTGCATCGGCAACCCGCGCACAGACTGGTGCTTCAAGACCGAGGCACAACCTGGCCTCGGCGGTCGCGCTTTGGGTTATCCAAGGGGCAAGGTGCTAGGTGGCTGTTCGTCGATCAACGGCATGATCTACATt                                                                                                                                                     >  3:226299/1‑148                                   GCACAGACTGGTGCTTCAAGACCGAGGCACAACCTGGCCTCGGCGGTCGCGCTTTGGGTTATCCAAGGGGCAAGGTGCTAGGTGGCTGTTCGTCGATCAACGGCATGATCTACATGCGCGGCCCAATGGTCGTAGTCGGCGGCTTGGCC                                                                                                                    >  2:322526/1‑149                                        GACTGGTGCTTCAAGACCGAGGCACAACCTGGCCTCGGCGGTCGCGCTTTGGGTTATCCAAGGGGCAAGGTGCTAGGTGGCTGTTCGTCGATCAACGGCATGATCTACATGCGCGGCCCAATGGTCGTAGTCGGCGGCTTGGCCCCGCA                                                                                                               <  1:322526/149‑1                                                                   CCTGGCCTCGGCGGTCGCGCTTTGGGTTATCCAAGGGGCAAGGTGCTAGGTGGCTGTTCGTCGATCAACGGCATGATCTACATGCGCGGCCCAATGGTCGTAGTCGGCGGCTTGGCCCCGCATGT                                                                                                            <  3:251264/125‑1                                                                   CCTGGCCTCGGCGGTCGCGCTTTGGGTTATCCAAGGGGCAAGGTGCTAGGTGGCTGTTCGTCGATCAACGGCATGATCTACATGCGCGGCCCAATGGTCGTAGTCGGCGGCTTGGCCCCGCATGT                                                                                                            >  4:251264/1‑125                                                                                                                                                      GCGCGGCCCAATGGTCGTAGTCGGCGGCTTGGCCCCGCATGTAGATCATGCCGTTGATCGACGAACAGCCACCTAGCACCTTGCCCCTTGGATAACCCAAAGCGCGACCGCCGAGGCCAGGTTGTGCCTCGGTCTTGAAGCACCAGTC  <  1:116862/148‑1                                                                                                                                                      GCGCGGCCCAATGGTCGTAGTCGGCGGCTTGGCCCCGCATGTAGATCATGCCGTTGATCGACGAACAGCCACCTAGCACCTTGCCCCTTGGATAACCCAAAGCGCGACCGCCGAGGCCAGGTTGTGCCTCGGTCTTGAAGCACCAGT   <  8:253640/147‑1                                                                                                                                                      GCGCGGCCCAATGGTCGTAGTCGGCGGCTTGGCCCCGCATGTAGATCATGCCGTTGATCGACGAACAGCCACCTAGCACCTTGCCCCTTGGATAACCCAAAGCGCGA                                           <  3:146422/107‑1                                                                                                                                                      GCGCGGCCCAATGGTCGTAGTCGGCGGCTTGGCCCCGCATGTAGATCATGCCGTTGATCGACGAACAGCCACCTAGCACCTTGCCCCTTGGATAACCCAAAGCGCGA                                           >  4:146422/1‑107                                                                                                                                                                                                                                                                                                           CCCGTTGGTTACCTCTACTGCATCGGCAACCCGCGCACAGACTGGTGCTTCAAGACCGAGGCACAACCTGGCCTCGGCGGTCGCGCTTTGGGTTATCCAAGGGGCAAGGTGCTAGGTGGCTGTTCGTCGATCAACGGCATGATCTACATGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/65099‑65250 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcgCGGCCCAATGGTCGTAGTCGGCGGCTTGGCCCCGCATGTAGATCATGCCGTTGATCGACGAACAGCCACCTAGCACCTTGCCCCTTGGATAACCCAAAGCGCGACCGCCGAGGCCAGGTTGTGCCTCGGTCTTGAAGCACCAGTC  <  NC_002947/65282‑65138

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