New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 284547215 (0.690)3 (0.150) 3/154 1.6 12.9% coding (250/273 nt) hypC hydrogenase 3 maturation protein HypC
?NZ_CP009273 = 2845516 27 (1.370)coding (22/1122 nt) hypD hydrogenase formation protein HypD
Rejected: Frequency below/above cutoff threshold.

CGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCGGA  >  NZ_CP009273/2845429‑2845603
                                                                                       |                                                                                       
cGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCgc                                                                                       >  1:410880‑M2/1‑88 (MQ=255)
   cTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCgcgcg                                                                                    <  1:554070‑M2/90‑6 (MQ=255)
         aaaCATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGaa                                                                              <  2:170923/90‑1 (MQ=255)
                 ttGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGAt                                                                      >  2:248321/1‑90 (MQ=255)
                  tGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATg                                                                     <  2:346298/90‑1 (MQ=255)
                       ttGAGCCGGATGCCGGCGCGCTGTTGGATGGCGAGGAAAAATACTGCTTTTTGTTGATGAATATCgcgcgccgacatccggctcaacgtc                                                                >  1:522252‑M2/1‑65 (MQ=255)
                                 ttccGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCAt                                                      >  2:185250/4‑90 (MQ=255)
                                     ggCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTgc                                                  >  2:157041/1‑90 (MQ=255)
                                      gcgcgcTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCgcgcgccgacatccggctcaacgtcaaacatgttttgtaa                                                 <  1:450638‑M2/90‑41 (MQ=255)
                                        gcgcTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGcgcg                                               >  2:294346/1‑90 (MQ=255)
                                        gcgcTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGcgcg                                               >  2:163963/1‑90 (MQ=255)
                                         cgcTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGa                                              >  2:182901/1‑90 (MQ=255)
                                          gcTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGaa                                             <  1:143158/90‑1 (MQ=255)
                                          gcTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGaa                                             >  2:201048/1‑90 (MQ=255)
                                               tgtATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCgcgcgccgacatccggctcaacgtcaaacatgttttgtaaggcgtcgag                                        >  2:258632‑M2/1‑41 (MQ=255)
                                                     gCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTcaca                                  >  2:428043/1‑90 (MQ=255)
                                                       gAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACAtc                                <  1:519370/90‑1 (MQ=255)
                                                       gAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACAtc                                <  2:559406/90‑1 (MQ=255)
                                                       gAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACAtc                                <  2:554070/90‑1 (MQ=255)
                                                       gAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACAtc                                <  1:160949/90‑1 (MQ=255)
                                                       gAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACAtc                                <  1:163963/90‑1 (MQ=255)
                                                       gAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACAtc                                <  1:355259/90‑1 (MQ=255)
                                                       gAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACAtc                                >  2:171338/1‑90 (MQ=255)
                                                         ggAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACAtctc                              >  2:146825/1‑90 (MQ=255)
                                                                  atGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACATCTCTCTTACAcc                     >  1:113388/1‑90 (MQ=255)
                                                                      gTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACATCTCTCTTACAccgccg                 >  2:479276/1‑90 (MQ=255)
                                                                                   atatCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCg    >  2:167421/1‑90 (MQ=255)
                                                                                    tatCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCgg   <  1:157041/90‑1 (MQ=255)
                                                                                    tatCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCgg   <  1:355519/90‑1 (MQ=255)
                                                                                    tatCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCgg   <  1:146825/90‑1 (MQ=255)
                                                                                    tatCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCgg   <  1:437104/90‑1 (MQ=255)
                                                                                     atCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCGGa  >  2:553328/1‑90 (MQ=255)
                                                                                       |                                                                                       
CGCCTTACAAAACATGTTTGACGTTGAGCCGGATGTCGGCGCGCTGTTGTATGGCGAGGAAAAATAATGCGTTTTGTTGATGAATATCGCGCGCCGGAACAGGTGATGCAGTTAATTGAGCATCTGCGCGAACGTGCTTCACATCTCTCTTACACCGCCGAACGCCCTCTGCGGA  >  NZ_CP009273/2845429‑2845603

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