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.

GCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAG  >  NZ_CP009273/2837229‑2837395
                                                                                        |                                                                              
gCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTa                                                                               <  1:388799/90‑1 (MQ=255)
 cACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTaa                                                                              <  2:366687/90‑1 (MQ=255)
   cgtcatgcggatattttactgtttacgggcgcggtcacccgtgcaatgcgatcccctgcgctgcgtgcgtggcagtccgcgccggTAAAc                                                                            <  1:395859‑M2/5‑1 (MQ=255)
        ggAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTaa                                                                       <  1:274495/90‑1 (MQ=255)
             tACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAatat                                                                  >  2:300526/1‑90 (MQ=255)
                aaaTTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCg                                                               >  1:176042/1‑90 (MQ=255)
                       gtttaCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCCGTAAACAGTAAAATATCCGCATGAcg                                                        <  1:395859‑M1/85‑1 (MQ=255)
                        gggTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGAcgc                                                       <  1:259355/90‑1 (MQ=255)
                        gggTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGAcgc                                                       <  1:98979/90‑1 (MQ=255)
                         ttaccGGCGCGGTCTGCCAAGCACTCCGCGACGGGGGTCGCGTTGCACGGGTGACCGCGCCGGTAAACAGGAAAATGTCCGCATGGcgcg                                                      >  1:48447/4‑90 (MQ=255)
                                                  gcgatcccctgcgctgcgtgcgtggcagtccgcgccggTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGaa                             >  2:94389‑M2/39‑90 (MQ=255)
                                                  cagcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTACTGCCGaa                             >  1:337037/1‑90 (MQ=255)
                                                  cagcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGaa                             >  2:601913/1‑90 (MQ=255)
                                                   cgatcccctgcgctgcgtgcgtggcagtccgcgccggTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAt                            <  1:94389‑M2/53‑1 (MQ=255)
                                                   agcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAt                            <  1:430559/90‑1 (MQ=255)
                                                    gcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATc                           <  1:358822/90‑1 (MQ=255)
                                                        agGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTc                       <  1:337064/90‑1 (MQ=255)
                                                        agGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTc                       <  2:405664/90‑1 (MQ=255)
                                                                        cGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcggc       >  2:264725/1‑90 (MQ=255)
                                                                           gTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGaa    >  1:212122/1‑90 (MQ=255)
                                                                             gACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAg  >  1:216587/1‑90 (MQ=255)
                                                                             gACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAg  >  2:74799/1‑90 (MQ=255)
                                                                                        |                                                                              
GCACCGTAGGAGATACAAATTTTCGGGTCCGGCGCGGACTGCCACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAG  >  NZ_CP009273/2837229‑2837395

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