New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	326549 =	44 (0.810)	8 (0.160)	8/276	2.7	16.0%	coding (699/1671 nt)	betA	choline dehydrogenase, a flavoprotein
	?	NC_000913	326652 =	43 (0.840)					coding (596/1671 nt)	betA	choline dehydrogenase, a flavoprotein
Rejected: Frequency below cutoff threshold.

GGTTTTGGTCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/326701‑326549 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acgcgcggtgCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGG  >  NC_000913/326652‑326799                                                                                                                                                                                                                                                                                                               GGTTTTGGTCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGG                                                                                                                                                        >  2:317280/1‑151  GTTTTGGTCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGT                                                                                                                                                       >  2:293937/1‑151    TTTGGTCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGC                                                                                                                                                     <  2:182940/151‑1         TCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCG                                                                                                                                                <  2:434014/151‑1                          TCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCGCGACGGCCCTGCGGCG                                                                                                                                <  1:481137/150‑1                               CCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGG                                                                                                                          <  1:167747/151‑1                                                    ACCGCGCGTGGCTATCGCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAAC                                                                                                     <  2:59995/151‑1                                                                            GCCAAATCGCGTCCTAACCTGACCATTCGTACTCCCGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCGCGAAGGCCCTGCGGCGTGGAGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAG                                                                               <  2:421899/149‑1                                                                            GCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCGCGACGGTCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGA                                                                              <  2:197982/150‑1                                                                                                         TACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCAC                                                 <  2:113751/150‑1                                                                                                                                                   CGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGT      >  2:863681/1‑151                                                                                                                                                     CGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTT    >  2:283596/1‑151                                                                                                                                                       GTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGG  >  2:44082/1‑151                                                                                                                                                       GTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGG  >  2:59024/1‑151                                                                                                                                                        TGCTGGCGCGCCGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGG  <  1:27115/150‑1                                                                                                                                                                                                                                                                                                               GGTTTTGGTCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/326701‑326549 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acgcgcggtgCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGG  >  NC_000913/326652‑326799

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