New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 4569286 =35 (1.130)4 (0.140) 4/152 1.9 11.3% intergenic (‑132/+30) mcrB/symE 5‑methylcytosine‑specific restriction endonuclease subunit McrB/endoribonuclease SymE
?NZ_CP009273 4569307 = 31 (1.100)intergenic (‑153/+9) mcrB/symE 5‑methylcytosine‑specific restriction endonuclease subunit McrB/endoribonuclease SymE
Rejected: Frequency below/above cutoff threshold.

GCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTT  >  NZ_CP009273/4569197‑4569370
                                                                                         |                                                                                    
tgtgcaaactgtcggcgcgtaaacaaaagcaggtgcaggcgtttattggagtgattgccggtaaacagaaagtcgcgtaacagcgtttac                                                                                      <  1:359314‑M1/1‑1 (MQ=255)
  tgcaaactgtcggcgcgtaaacaaaagcaggtgcaggcgtttattggagtgattgccggtaaacagaaagtcgcgtaacagcgtttaCTc                                                                                    <  1:787707‑M1/3‑1 (MQ=255)
  tgcaaactgtcggcgcgtaaacaaaagcaggtgcaggcgtttattggagtgattgccggtaaacagaaagtcgcgtaacagcgtttaCTc                                                                                    >  2:226214‑M1/88‑90 (MQ=255)
   gcaaactgtcggcgcgtaaacaaaagcaggtgcaggcgtttattggagtgattgccggtaaacagaaagtcgcgtaacagcgtttaCTcc                                                                                   <  2:148756‑M1/4‑1 (MQ=255)
    cGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTccc                                                                                  <  1:699964/90‑1 (MQ=255)
       actgtcggcgcgtaaacaaaagcaggtgcaggcgtttattggagtgattgccggtaaacagaaagtcgcgtaacagcgtttaCTCCCggt                                                                               <  2:8353‑M1/8‑1 (MQ=255)
             aCAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTaa                                                                         >  1:351539/1‑90 (MQ=255)
               aaTAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAAcc                                                                       >  2:591397/1‑90 (MQ=255)
               aaTAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAAcc                                                                       <  1:11783/90‑1 (MQ=255)
               aaTAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAAcc                                                                       >  1:700306/1‑90 (MQ=255)
                aTAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCg                                                                      >  1:536855/1‑90 (MQ=255)
                   gATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGa                                                                   <  1:546211/90‑1 (MQ=255)
                      ttttACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTa                                                                <  2:67131/90‑1 (MQ=255)
                        ttACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCg                                                                      <  2:270472/82‑1 (MQ=255)
                        ttACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCg                                                                      >  1:270472/1‑82 (MQ=255)
                          aCTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACg                                                            <  2:119342/90‑1 (MQ=255)
                           agcaggtgcaggcgtttattggagtgattgccggtaaacagaaagtcgcgtaacagcgtttaCTccc                                                                                  <  2:133272‑M1/5‑1 (MQ=255)
                           agcaggtgcaggcgtttattggagtgattgccggtaaacagaaagtcgcgtaacagcgtttaCTccc                                                                                  >  1:133272‑M1/63‑67 (MQ=255)
                                          ggTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCggtg                                                                              >  2:857885/1‑56 (MQ=255)
                                          ggTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCggtg                                                                              <  1:857885/56‑1 (MQ=255)
                                          ggTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTc                                            <  1:413268/90‑1 (MQ=255)
                                               atatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGttt                                       <  1:304494/90‑1 (MQ=255)
                                                tatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTa                                      >  1:519312/1‑90 (MQ=255)
                                                    tACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCgg                                  >  2:680356/1‑90 (MQ=255)
                                                    tACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCgg                                  >  1:98774/1‑90 (MQ=255)
                                                       tGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCaa                               <  2:524973/90‑1 (MQ=255)
                                                        gAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                              <  2:581323/90‑1 (MQ=255)
                                                         aGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATc                             >  1:499961/1‑90 (MQ=255)
                                                         aGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATc                             <  1:4798/90‑1 (MQ=255)
                                                         aGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATc                             >  2:414217/1‑90 (MQ=255)
                                                          gggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCa                            >  2:549092/1‑90 (MQ=255)
                                                           ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCAc                           <  2:700306/90‑1 (MQ=255)
                                                            gggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACt                          <  1:705014/90‑1 (MQ=255)
                                                              gAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTcc                        <  2:536855/90‑1 (MQ=255)
                                                                aaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCaa                      >  1:582022/1‑90 (MQ=255)
                                                                   gaaagtcgcgtaacagcgtttaCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATaa                   <  1:797164‑M1/68‑1 (MQ=255)
                                                                   gaaagtcgcgtaacagcgtttaCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATaa                   >  2:797164‑M1/23‑90 (MQ=255)
                                                                    ttCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATaaa                  >  2:803757/1‑90 (MQ=255)
                                                                          gcgtaacagcgtttaCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTg            <  1:529457‑M1/75‑1 (MQ=255)
                                                                            aGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCa          >  1:135422/1‑90 (MQ=255)
                                                                             gTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAc         <  1:725728/90‑1 (MQ=255)
                                                                              tCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAcc        <  1:648308/90‑1 (MQ=255)
                                                                                aGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTg      >  1:515486/1‑90 (MQ=255)
                                                                                 gTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGc     >  2:140785/1‑90 (MQ=255)
                                                                                  tGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCt    <  1:494689/90‑1 (MQ=255)
                                                                                    gtttaCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCttt  <  1:335216‑M1/85‑1 (MQ=255)
                                                                                         |                                                                                    
GCAACGTACTGAAACAATAGATTTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTT  >  NZ_CP009273/4569197‑4569370

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