New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 27822328 (1.300)3 (0.150) 3/148 1.7 9.8% coding (235/909 nt) yagE 2‑keto‑3‑deoxygluconate aldolase
?NZ_CP009273 = 278229 30 (1.540)coding (241/909 nt) yagE 2‑keto‑3‑deoxygluconate aldolase
Rejected: Frequency below/above cutoff threshold.

GCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATC  >  NZ_CP009273/278134‑278312
                                                                                         |                                                                                         
gCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGC                                                                                           >  2:372970‑M1/1‑90 (MQ=255)
gCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGC                                                                                           >  1:105468‑M1/1‑90 (MQ=255)
gCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGC                                                                                           <  1:595950‑M1/90‑1 (MQ=255)
     ggCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATc                                                                                      <  2:295039/90‑1 (MQ=255)
     ggCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATc                                                                                      >  2:412773/1‑90 (MQ=255)
          gTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcac                                                                                 <  1:118392‑M2/90‑5 (MQ=255)
             ctcCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcaccgg                                                                              <  2:219651‑M2/90‑8 (MQ=255)
             ctcCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcaccgg                                                                              <  2:506068‑M2/90‑8 (MQ=255)
                    cTCGGCGCCGAAGAGCGTAAGGCCATTGCCCGCTTTGCTATCGATCATGTCGGTCGCCGCGTGCCGGTGCTGATCGgcaccggcacgcga                                                                       >  2:295379‑M2/1‑76 (MQ=255)
                      cGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCaa                                                                     >  1:401137/1‑90 (MQ=255)
                       ggCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAAc                                                                    <  1:231547/90‑1 (MQ=255)
                                  gCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAAc                                                         >  2:624434/1‑90 (MQ=255)
                                      aaaGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATc                                                     <  1:7414/90‑1 (MQ=255)
                                       aaGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCg                                                    >  1:477126/1‑90 (MQ=255)
                                          ccATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCAcc                                                                       >  1:161058/1‑68 (MQ=255)
                                          ccATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCAcc                                                                       <  2:161058/68‑1 (MQ=255)
                                            aTTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTc                                               >  2:310338/1‑90 (MQ=255)
                                            aTTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTc                                               >  2:482639/1‑90 (MQ=255)
                                             ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcaccggcacgcgacgatcgacat                                                             >  2:307218‑M2/1‑51 (MQ=255)
                                             ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcaccggcacgcgacgatcgacat                                                             <  1:307218‑M2/75‑25 (MQ=255)
                                                  cGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagccag                                         >  1:26688/1‑90 (MQ=255)
                                                   gCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagccagc                                        >  1:644145/1‑90 (MQ=255)
                                                         cTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACgcgc                                  <  1:30609/90‑1 (MQ=255)
                                                           aTCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCgcag                                <  1:412773/90‑1 (MQ=255)
                                                              gATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGcagcag                             >  2:550694/1‑90 (MQ=255)
                                                                  aTGTCGATCGTCGCGTGCCGGTGCcgatcagcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaa                                        <  1:307218‑M1/75‑52 (MQ=255)
                                                                  aTGTCGATCGTCGCGTGCCGGTGCcgatcagcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaa                                        >  2:307218‑M1/1‑24 (MQ=255)
                                                                  aTGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACg                                                                   >  2:201616/1‑48 (MQ=255)
                                                                  aTGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACg                                                                   <  1:201616/48‑1 (MQ=255)
                                                                  aTGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATc                                                     >  2:315235/1‑62 (MQ=255)
                                                                  aTGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATc                                                     <  1:315235/62‑1 (MQ=255)
                                                                       gATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCAt                                                      <  2:356398/56‑1 (MQ=255)
                                                                       gATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCAt                                                      >  1:356398/1‑56 (MQ=255)
                                                                          cgtcgCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGAc                 <  2:105468/90‑1 (MQ=255)
                                                                          cgtcgCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGAc                 <  1:482639/90‑1 (MQ=255)
                                                                            tcgCGTGCCGGTGCcgatcagcaccggcacgcggcgaccgacatgatcgatagcaaagcgggcaatggccttacgctcttcggcgccgag               >  2:295379‑M1/1‑14 (MQ=255)
                                                                                   ccggtgccgatcagcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggag        <  2:219651‑M1/90‑84 (MQ=255)
                                                                                   ccggtgccgatcagcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggag        <  2:506068‑M1/90‑84 (MQ=255)
                                                                                     ggTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtggt      >  1:8087/1‑90 (MQ=255)
                                                                                      gtgccgatcagcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaac     <  1:118392‑M1/90‑87 (MQ=255)
                                                                                      gTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtggtg     >  1:108077/1‑90 (MQ=255)
                                                                                         cTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATc  <  2:293979/90‑1 (MQ=255)
                                                                                         cTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATc  <  2:401137/90‑1 (MQ=255)
                                                                                         |                                                                                         
GCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATC  >  NZ_CP009273/278134‑278312

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