New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 2837257 =27 (1.240)3 (0.150) 3/152 1.6 12.5% coding (313/768 nt) hycG formate hydrogenlyase subunit HycG
?NZ_CP009273 2837317 = 18 (0.920)coding (253/768 nt) hycG formate hydrogenlyase subunit HycG
Rejected: Frequency below/above cutoff threshold.

CGTCATGCGGATATTTTACTGTTTACCGGCGCGGTCACCCGTGCAATGCGATCCCCTGCGCTGCGTGCGTGGCAGTCCGCGCCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/2837341‑2837257
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccgcgccggTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAT  >  NZ_CP009273/2837317‑2837369
                                                                                                                                          
CGTCATGCGGATATTTTACTGTTTACGGGCGCGGTCACCCGTGCAATGCGATCCCCTGCGCTGCGTGCGTGGCAGTCCGCGCCGGTAAAC                                                  <  1:395859/90‑1
                                               GCGATCCCCTGCGCTGCGTGCGTGGCAGTCCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAA   >  2:94389/1‑90
                                                CGATCCCCTGCGCTGCGTGCGTGGCAGTCCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAT  <  1:94389/90‑1
                                                                                                                                          
CGTCATGCGGATATTTTACTGTTTACCGGCGCGGTCACCCGTGCAATGCGATCCCCTGCGCTGCGTGCGTGGCAGTCCGCGCCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/2837341‑2837257
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccgcgccggTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAT  >  NZ_CP009273/2837317‑2837369

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