New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 2416060 =37 (0.790)3 (0.070) 3/146 3.1 7.9% coding (874/894 nt) yfcH TIGR01777 family oxidoreductase
?NZ_CP009273 2416107 = 38 (0.940)intergenic (+27/+21) yfcH/rpnB TIGR01777 family oxidoreductase/recombination‑promoting nuclease RpnB
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAGTACTGGTATTAGGTGGACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTGCGCGATAAGA  >  NZ_CP009273/2415971‑2416149
                                                                                         |                                                                                         
cAGTACTGGTATTAGGTGGACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGgcgc                                                                                           >  2:559785/1‑90 (MQ=255)
       ggTATTAGGTGGACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGa                                                                                    >  2:789689/1‑90 (MQ=255)
       ggTATTAGGTGGACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGa                                                                                    >  2:1065783/1‑90 (MQ=255)
       ggTATTAGGTGGACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGa                                                                                    >  2:700673/1‑90 (MQ=255)
          aTTAGGTGGACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGt                                                                                 <  2:775567/90‑1 (MQ=255)
            tAGGTGGACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCg                                                                               >  2:321310/1‑90 (MQ=255)
                 ggACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGc                                                                          >  1:1157357/1‑90 (MQ=255)
                 ggACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGc                                                                          >  1:1153238/1‑90 (MQ=255)
                 ggACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGc                                                                          >  2:440065/1‑90 (MQ=255)
                 ggACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGc                                                                          >  2:427691/1‑90 (MQ=255)
                 ggACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGc                                                                          >  2:1115047/1‑90 (MQ=255)
                      aCGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGAtg                                                                     >  2:236785/1‑90 (MQ=255)
                      aCGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGAtg                                                                     <  2:1176735/90‑1 (MQ=255)
                       cgcgcgCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGAtgt                                                                    <  2:179784/90‑1 (MQ=255)
                                         cTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATcc                                                  >  2:696827/1‑90 (MQ=255)
                                         cTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATcc                                                  >  2:1132145/1‑90 (MQ=255)
                                         cTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATcc                                                  >  1:1098741/1‑90 (MQ=255)
                                         cTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATcc                                                  >  1:18366/1‑90 (MQ=255)
                                                     ggTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTc                                      <  1:321310/90‑1 (MQ=255)
                                                         ttGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCgg                                  <  1:516321/90‑1 (MQ=255)
                                                           gCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGgtg                                >  2:659847/1‑90 (MQ=255)
                                                               ttCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTAc                            >  1:410394/1‑90 (MQ=255)
                                                                 cGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAg                          <  1:832817/90‑1 (MQ=255)
                                                                 cGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAg                          <  1:1065783/90‑1 (MQ=255)
                                                                       tACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAg                    <  1:1136978/90‑1 (MQ=255)
                                                                                gaagaGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTgc           <  1:306376/90‑1 (MQ=255)
                                                                                  agagGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTgcgc         <  1:659847/90‑1 (MQ=255)
                                                                                    agGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTGCGCGa       <  1:310163/90‑1 (MQ=255)
                                                                                    agGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTGCGCGa       >  2:5743/1‑90 (MQ=255)
                                                                                     ggCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTgcgc         <  2:810388/87‑1 (MQ=255)
                                                                                     ggCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTgcgc         <  2:780294/87‑1 (MQ=255)
                                                                                     ggCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTgcgc         >  1:780294/1‑87 (MQ=255)
                                                                                     ggCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTgcgc         >  1:810388/1‑87 (MQ=255)
                                                                                     ggCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTGCGCGAt      >  1:21493/1‑90 (MQ=255)
                                                                                     ggCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTGCGCGAt      >  2:343952/1‑90 (MQ=255)
                                                                                     ggCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTGCGCGAt      >  2:1267502/1‑90 (MQ=255)
                                                                                     ggCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTGCGCGAt      >  2:782084/1‑90 (MQ=255)
                                                                                         cTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTGCGCGATAAGa  <  2:14815/90‑1 (MQ=255)
                                                                                         |                                                                                         
CAGTACTGGTATTAGGTGGACAACGCGCGCTGCCTAAACGGCTGGAAGAAGCGGGTTTTGCGTTTCGCTGGTACGATTTAGAAGAGGCGCTGGCGGATGTCGTTCGCTGATGTGGTTTACAGCAAACATCCGCCAGTTAACTCCCGGTGTTACAGGATTAGTGGCTTTGCGCGATAAGA  >  NZ_CP009273/2415971‑2416149

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