New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 6610358 (1.150)3 (0.060) 3/246 NT 5.0% coding (1161/1653 nt) betA‑I choline dehydrogenase
?NC_002947 = 66147 58 (1.210)coding (1205/1653 nt) betA‑I choline dehydrogenase

GGCCGCTGGCCATGGCGCCGAGCCAGTTAGGCGCGTTCGTGCGTTCAAGCCCCGAACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAAC‑GCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/65952‑66103
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcggtgCTGCAGATGTCGATGCGCCCGCGGCTGGCCGGGCGCAGGTTGCACACCGAGG  <  NC_002947/66147‑66096
                                                                                                                                                                                                             
ggactCTGGCCATGGCGCCTAGCCAGTTAGGCGCGTTCGTTTGTTCAAGCCCCGAACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAAC‑GCTTCGGTGAGCCGTAGCATCAGTTCCCGGCCTTTACTGCCTCG                                                         <  1:179888/144‑1
  CCGCTGGCCATGGCGCCGAGCCAGTTAGGCGCGTTCGTGCGTTCAAGCCCCGAACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAAC‑GCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGT                                                       >  2:252951/1‑149
  CCGCTGGCCATGGCGCCGAGCCAGTTAGGCGCGTTCGTGCGTTCAAGCCCCGAACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAAC‑GCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGT                                                       <  6:279779/149‑1
  CCGCTGGCCATGGCGCCGAGCCAGTTAGGCGCGTTCGTGCGTTCAAGCCCCGAACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAAC‑GCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGT                                                       >  8:235999/1‑149
                                                       ACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAAC‑GCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGCTGCAGATGTCGATGCGCCCGCGGCTGGCCGGGCGCAGGTTGCACACCGAGG  <  6:220211/149‑1
                                                             CACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAAC‑GCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGCTGCAGATGTCGATGCGCCCGCGGCTGGCCGGGCGCAGGTTGCACACCGAG   <  7:79595/142‑1
                                                             CACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAAC‑GCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTGCTGCAGATGTCGATGCGCCCGCGGCTGGCCGGGCGCAGGTTGCACACCGAG   >  8:79595/1‑142
                                                                                               TCACTGGAACGGCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGT                                                       <  1:232520/57‑1
                                                                                               TCACTGGAAC‑GCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGT                                                       >  2:232520/1‑56
                                                                                                                                                                                                             
GGCCGCTGGCCATGGCGCCGAGCCAGTTAGGCGCGTTCGTGCGTTCAAGCCCCGAACAGGCCACCGCCAACCTTCAGTATCACGTGCAGCCGCTGTCACTGGAAC‑GCTTCGGTGAGCCGTTGCATCAGTTCCCGGCCTTTACTGCCTCGGTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/65952‑66103
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcggtgCTGCAGATGTCGATGCGCCCGCGGCTGGCCGGGCGCAGGTTGCACACCGAGG  <  NC_002947/66147‑66096

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