New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 32226938 (2.480)3 (0.210) 3/238 1.3 8.9% coding (690/1671 nt) betA choline dehydrogenase
?NZ_CP009273 = 322352 26 (1.840)coding (607/1671 nt) betA choline dehydrogenase
Rejected: Frequency below/above cutoff threshold.

CCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/322124‑322269
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caccgcgcgtGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGGGCGTCGAATGGCTGGAAGGCGACAGCACCATCCCAACCCG  <  NZ_CP009273/322352‑322219
                                                                                                                                                                                                                                                                                        
CCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCAC                                                                                                                                               >  2:26805/1‑139
                                                                                        CGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTCTGACCGATCACATCATTTTTGACGG                                                        >  1:258001/1‑138
                                                                                                            TGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTT                                                               <  1:152029/111‑1
                                                                                                            TGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTT                                                               >  2:152029/1‑111
                                                                                                                                            GCGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGGGCGTCGAATGGCTGGAAGGCGACAGCACCATCCCAACCC   >  1:145134/1‑139
                                                                                                                                             CGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACTCGCGGTGGGCGTCGAATGGCTGGAAGGCGACAGCACCATCCCAACCCG  >  2:84589/1‑139
                                                                                                                                             CGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGGGCGTCGAATGGCTGGAAGGCGACAGCACCATCCCAACCCG  >  2:138837/1‑139
                                                                                                                                             CGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGGGCGTCGAATGGCTGGAAGGCGACAGCACCATCCCAACCCG  >  2:85192/1‑139
                                                                                                                                             CGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGGGCGTCGAATGGCTGGAA                         <  1:141233/116‑1
                                                                                                                                             CGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGGGCGTCGAATGGCTGGAA                         >  2:141233/1‑116
                                                                                                                                             CGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACG                                                                                    >  1:66043/1‑57
                                                                                                                                             CGCGTGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACG                                                                                    <  2:66043/57‑1
                                                                                                                                                 TGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGCCCGAT                                                                         >  1:240162/1‑64
                                                                                                                                                 TGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGAT                                                                         <  2:240162/64‑1
                                                                                                                                                                                                                                                                                        
CCAGCAGTTCAGCGTTGCCGACGCCGGAGCGTTGCAGGATCTGCGGTGAGGCAATCGCGCCTGCACATAACAGCACTTCTTTGTTGGCCGTTGCGCGGGTTGGGATGGTGCTGTCGCCTTCCAGCCATTCGACGCCCACCGCGCGT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/322124‑322269
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caccgcgcgtGGCTATCTCGATCAGGCCAAATCGCGTCCTAACCTGACCATTCGTACTCACGCTATGACCGATCACATCATTTTTGACGGCAAACGCGCGGTGGGCGTCGAATGGCTGGAAGGCGACAGCACCATCCCAACCCG  <  NZ_CP009273/322352‑322219

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