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.

TTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACA  >  NZ_CP009273/4569219‑4569395
                                                                                        |                                                                                        
ttttACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTa                                                                                         <  2:67131/90‑1 (MQ=255)
   aaagcaggtgcaggcgtttattggagtgattgccggtaaacagaaagtcgcgtaacagcgtttactcccggttaaccaccgggagTAAAc                                                                                      <  1:335216‑M2/5‑1 (MQ=255)
    aCTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACg                                                                                     <  2:119342/90‑1 (MQ=255)
             caggcgtttattggagtgattgccggtaaacagaaagtcgcgtaacagcgtttactcccggttaaccaccgggagTAAACGCTGTTAcgc                                                                            <  1:529457‑M2/15‑1 (MQ=255)
                    ttattggagtgattgccggtaaacagaaagtcgcgtaacagcgtttactcccggttaaccaccgggagTAAACGCTGTTACGCGACTTTc                                                                     >  2:797164‑M2/69‑90 (MQ=255)
                    ttattggagtgattgccggtaaacagaaagtcgcgtaacagcgtttactcccggttaaccaccgggagTAAACGCTGTTACGCGACTTTc                                                                     <  1:797164‑M2/22‑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                                                           >  1:98774/1‑90 (MQ=255)
                              tACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCgg                                                           >  2:680356/1‑90 (MQ=255)
                                 tGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCaa                                                        <  2:524973/90‑1 (MQ=255)
                                  gAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAAt                                                       <  2:581323/90‑1 (MQ=255)
                                   aGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATc                                                      >  2:414217/1‑90 (MQ=255)
                                   aGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATc                                                      >  1:499961/1‑90 (MQ=255)
                                   aGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATc                                                      <  1:4798/90‑1 (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)
                                                                           ggtTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCg              <  2:723575/90‑1 (MQ=255)
                                                                            gtTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGa             <  2:268846/90‑1 (MQ=255)
                                                                            gtTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGa             >  2:537706/1‑90 (MQ=255)
                                                                            gtTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGa             <  1:442733/90‑1 (MQ=255)
                                                                             ttAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGAc            <  1:469206/90‑1 (MQ=255)
                                                                              tAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACa           <  2:833765/90‑1 (MQ=255)
                                                                              tAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACa           <  2:321510/90‑1 (MQ=255)
                                                                                accgggagTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGt         <  2:8353‑M2/82‑1 (MQ=255)
                                                                                   gggagTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCt                             <  2:133272‑M2/62‑1 (MQ=255)
                                                                                   gggagTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCt                             >  1:133272‑M2/6‑67 (MQ=255)
                                                                                    ggagTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGc     <  2:148756‑M2/86‑1 (MQ=255)
                                                                                     gagTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGca    <  1:787707‑M2/87‑1 (MQ=255)
                                                                                     gagTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGca    >  2:226214‑M2/4‑90 (MQ=255)
                                                                                       gTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcaca  <  1:359314‑M2/89‑1 (MQ=255)
                                                                                        |                                                                                        
TTTTACTGCAAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACA  >  NZ_CP009273/4569219‑4569395

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