New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3622466 =NA (NA)16 (0.180) 7/170 NT 100% coding (3275/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 3622497 = 0 (0.000)coding (3306/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

GCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCG  >  NC_000913/3622398‑3622596
                                                                                                   |                                                                                                   
tcaggccgcacgatgccagccagcggcggctttcctcactcgcccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgg                                                                                                     <  2:1067979‑M2/1‑1 (MQ=255)
tcaggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgg                                                                                                     >  1:1192789‑M2/100‑100 (MQ=255)
tcaggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgg                                                                                                     <  1:1680324‑M2/1‑1 (MQ=255)
 caggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcggg                                                                                                    <  1:790608‑M2/2‑1 (MQ=255)
 caggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcggg                                                                                                    <  2:2079999‑M2/2‑1 (MQ=255)
 caggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcggg                                                                                                    <  1:1648716‑M2/2‑1 (MQ=255)
  aggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGa                                                                                                   <  2:739468‑M2/3‑1 (MQ=255)
  aggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGa                                                                                                   <  1:1936804‑M2/3‑1 (MQ=255)
  aggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGa                                                                                                   <  2:2031417‑M2/3‑1 (MQ=255)
  aggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGa                                                                                                   <  2:1821079‑M2/3‑1 (MQ=255)
  aggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGa                                                                                                   >  1:4226‑M2/98‑100 (MQ=255)
  aggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGa                                                                                                   <  1:1251798‑M2/3‑1 (MQ=255)
   ggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgcgcatctgcaccagcaccggcgGGaa                                                                                                  >  1:512439‑M2/97‑100 (MQ=255)
   ggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcggg                                                                                                    <  1:1349889‑M2/2‑1 (MQ=255)
   ggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGaa                                                                                                  <  2:2396391‑M2/4‑1 (MQ=255)
   ggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGaa                                                                                                  >  1:2198163‑M2/97‑100 (MQ=255)
   ggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGaa                                                                                                  >  2:1674043‑M2/97‑100 (MQ=255)
   ggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGaa                                                                                                  >  2:1491884‑M2/97‑100 (MQ=255)
   ggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGaa                                                                                                  >  2:113139‑M2/97‑100 (MQ=255)
   ggccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGa                                                                                                   >  1:760465‑M2/97‑99 (MQ=255)
   ggccgcacgaggccagccagcggcatctatcctcactcacccagtcagccaggatttcacttaccagccggtcgagcatctgcaccagcaccggcgGGaa                                                                                                  <  2:734960‑M2/4‑1 (MQ=255)
    gccgcacgatgccagccagcggcggctttcctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccagcaccggcgGGaa                                                                                                  <  1:2235454‑M2/4‑1 (MQ=255)
     ccgaacgatgccaggcagcggcggctgtcctcactcacgcggtcagccaggatgtcactttccagccggtcgagcatctgcaccagcaccggcggg                                                                                                    <  2:782671‑M2/2‑1 (MQ=255)
                               ttcctcactcacccggacagccaggatttcactttccagccggtcgagcatctgcaccagccacggcgGGAAAGGGAAATCCTGGCTGACCGGgtgagtg                                                                      >  1:814714‑M2/69‑100 (MQ=255)
                                  cctcactcacccggtcagccaggatttcactttccagccggtcgagcatctgcaccgtgccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAgg                                                                   >  1:142238‑M2/66‑100 (MQ=255)
                                                       gatttcactttccagccggtcgagcatctgcaccagcaccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCAt                                              >  1:83110‑M2/45‑100 (MQ=255)
                                                            cactttccagccggtcgagcatctgcaccagcaccggcgGGAAAGTGAAATCATGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGc                                         >  1:2533694‑M2/40‑100 (MQ=255)
                                                              ctttccagccggtcgagcatctgcaccagcaccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCgg                                       >  1:867273‑M2/38‑100 (MQ=255)
                                                              ctttccagccggtcgagcatctgcaccagcaccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCgg                                       >  1:989722‑M2/38‑100 (MQ=255)
                                                              ctttccagccggtcgagcatctgcaccagcaccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCgg                                       >  1:724499‑M2/38‑100 (MQ=255)
                                                              ctttccagccggtcgagcatctgcaccagcaccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCgg                                       >  2:1586072‑M2/38‑100 (MQ=255)
                                                               ctttccagccggtcgagcatctgcaccgtgccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGc                                      >  1:2027409‑M2/37‑100 (MQ=255)
                                                               ctttccagccggtcgagcatctgcaccgtgccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGc                                      >  2:2548876‑M2/37‑100 (MQ=255)
                                                                 tccagccggtcgagcatctgcaccagcaccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCt                                    >  1:1546141‑M2/35‑100 (MQ=255)
                                                                  tccagccggtcgagcatctgcaccgtgccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTg                                   >  1:2510255‑M2/34‑100 (MQ=255)
                                                                  ccagccggtcgagcatctgcaccagcaccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTg                                   >  1:100275‑M2/34‑100 (MQ=255)
                                                                  ccagccggtcgagcatctgcaccagcaccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTg                                   >  1:2089379‑M2/34‑100 (MQ=255)
                                                                   ccagccggtcgagcatctgcaccgtgccggcgGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGa                                  >  1:234987‑M2/33‑100 (MQ=255)
                                                                     gccggtcgagcatctgcaccagcaccggctGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGAcc                                >  1:209680‑M2/31‑100 (MQ=255)
                                                                                                   GGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCg  <  1:2018531‑M2/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCG  >  NC_000913/3622398‑3622596

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