New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 326549 =35 (1.010)7 (0.220) 7/276 1.9 15.8% coding (699/1671 nt) betA choline dehydrogenase, a flavoprotein
?NC_000913 326652 = 42 (1.300)coding (596/1671 nt) betA choline dehydrogenase, a flavoprotein
Rejected: Frequency below cutoff threshold.

GCAGGAAGGTTTTGGTCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/326708‑326549
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acgcgcggtgCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCAC  >  NC_000913/326652‑326752
                                                                                                                                                                                                                                                                     
GCAGGAAGGTTTTGGTCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAA                                                                                                                >  2:252324/1‑151
   GGAAGGTTTTGGTCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACG                                                                                                              >  1:169593/1‑150
        GTTTTGGTCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGg                                                                                                        >  1:237138/1‑150
          TTTGGTCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGC                                                                                                      <  2:518116/151‑1
               TCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCG                                                                                                 <  1:275156/151‑1
                        TCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCGCGACGGCCC                                                                                        <  2:110997/151‑1
                             CCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCGCGACGGCCCTGCGG                                                                                   <  1:33561/151‑1
                                             CCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCC                                                                   <  2:964/151‑1
                                                                                   CCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGAT                              <  2:511045/150‑1
                                                                                                               TACTCACGCCATTACCGATCACATCATTTTTGACGGCAAACGCGCGGTGCTGGCGCGACGGCCCAGCGGCGTGACGGTGCGCTCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCAC  <  1:254090/150‑1
                                                                                                                                                                                                                                                                     
GCAGGAAGGTTTTGGTCCGATGGATCGCACCGTCACGCCGCAGGGCCGTCGCGCCAGCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/326708‑326549
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑acgcgcggtgCTGGCGCGACGGCCCTGCGGCGTGACGGTGCGATCCATCGGACCAAAACCTTCCTGCTGATAACCGTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCAC  >  NC_000913/326652‑326752

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