New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3620350NA (NA)13 (0.150) 8/166 NT 100% coding (1159/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 = 3620425 0 (0.000)coding (1234/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

GCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGC  >  NC_000913/3620321‑3620448
                             |                                                                                                  
gCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                              >  2:2410645‑M2/1‑100 (MQ=255)
gCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                              >  1:2315310‑M2/1‑100 (MQ=255)
gCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGT                              <  1:2546743‑M2/100‑1 (MQ=255)
  tACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                            <  1:1359125‑M2/100‑1 (MQ=255)
  tACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                            >  2:198368‑M2/1‑100 (MQ=255)
  tACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGG                            <  1:315241‑M2/100‑1 (MQ=255)
   acacGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                           >  2:1466345‑M2/1‑100 (MQ=255)
   acacGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGT                           >  2:1385249‑M2/1‑100 (MQ=255)
    cacGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTG                          >  1:1320576‑M2/1‑100 (MQ=255)
    cacGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTG                          >  1:187247‑M2/1‑100 (MQ=255)
    cacGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTG                          >  1:1064507‑M2/1‑100 (MQ=255)
        tATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGCCGAAGCCGGGCTGAAGCGGGTGGTG                          <  1:567876‑M2/96‑1 (MQ=255)
                         ncgCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactga     <  1:1897001‑M2/99‑21 (MQ=255)
                         ccGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGTGCGAAGCCGGGCTGAAGCGGGTGGTGAggcggtctttctcatagtga     <  1:2015752‑M2/100‑21 (MQ=255)
                         ccGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactga     <  1:638705‑M2/100‑21 (MQ=255)
                           gCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgata   <  2:342580‑M2/100‑23 (MQ=255)
                            catcaccatcacCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAtgcggtctttctcatactgatac  <  2:2091674‑M2/100‑24 (MQ=255)
                             |                                                                                                  
GCTACACGTATCAGTATGAGAAAGACCGCATCACCATCACCGACAGCCTGGACCGCCGTGAAGTGCTGCACACGCAGGGCGAAGCCGGGCTGAAGCGGGTGGTGAAAAAGGAACACGCGGACGGCAGC  >  NC_000913/3620321‑3620448

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.