New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3620220NA (NA)19 (0.220) 10/166 NT 100% coding (1029/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 = 3620255 0 (0.000)coding (1064/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTA  >  NC_000913/3620147‑3620319
                                                                         |                                                                                                   
nCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCacac                                                                           >  1:1650613‑M2/2‑100 (MQ=255)
GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCgtcactc                                                                           >  1:1565862‑M2/1‑98 (MQ=255)
GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCacac                                                                           <  1:1839209‑M2/100‑1 (MQ=255)
  GGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACgg                                                                         >  1:1466337‑M2/1‑100 (MQ=255)
   GCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACggg                                                                        <  1:615283‑M2/100‑1 (MQ=255)
   GCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACggg                                                                        <  2:341867‑M2/100‑1 (MQ=255)
   GCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAAGACCGGGGCCGGATGGTGGCGCACCGTCACACggg                                                                        <  1:1883932‑M2/100‑1 (MQ=255)
    cGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGC                                                                       <  1:17229‑M2/100‑1 (MQ=255)
    cGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGC                                                                       <  2:1206472‑M2/100‑1 (MQ=255)
    cGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGC                                                                       >  2:206709‑M2/1‑100 (MQ=255)
    cGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGC                                                                       >  2:2094652‑M2/1‑100 (MQ=255)
    CGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACggg                                                                        >  1:2030433‑M2/1‑99 (MQ=255)
     gcactGGCGGTGGTCTATGACCGTCGTGGCAAACCGGTGCGCAGCTTTACCTACGATGATAAATACCCGGGCCGGATGGTGGCGCACCGTCACACGGGCC                                                                      >  1:2083089‑M2/3‑100 (MQ=255)
     gAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATGCCGGGGCCGGATGGTGGCGCACCGTCACACGGGCC                                                                      >  1:1868940‑M2/1‑100 (MQ=255)
     gAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCC                                                                      >  2:1945682‑M2/1‑100 (MQ=255)
       aCTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCG                                                                     >  1:618494‑M2/1‑99 (MQ=255)
       aCTGGCGGTGGTATATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGG                                                                    >  1:2511713‑M2/1‑100 (MQ=255)
        cTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                   <  2:187248‑M2/100‑1 (MQ=255)
        cTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                   >  1:1805844‑M2/1‑100 (MQ=255)
          ggCGGTGGTGTATGAGCGTAGTGGCAAGCAGGTGCGGAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                   >  2:63142‑M2/1‑98 (MQ=255)
          ggCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                  <  1:174399‑M2/99‑1 (MQ=255)
                  tgtATGACCGTAGTGGCAACCAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCgcaactcc                                                          <  1:1495118‑M2/99‑9 (MQ=255)
                             aGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGGTGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaag                                              >  1:1870268‑M2/1‑80 (MQ=255)
                              gTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagt                                             <  1:537683‑M2/100‑22 (MQ=255)
                                  cAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGAGCCGGCCccggtatttctcatcgtaagtaaag                                         <  1:1591702‑M2/100‑26 (MQ=255)
                                                                 aaaTACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcataca          <  1:220162‑M2/100‑57 (MQ=255)
                                                                     aCCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccac      <  2:1613455‑M2/100‑61 (MQ=255)
                                                                      ccGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccacc     <  1:2350942‑M2/100‑62 (MQ=255)
                                                                      ccGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccacc     <  1:1084542‑M2/100‑62 (MQ=255)
                                                                      ccGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccac      <  2:720403‑M2/99‑61 (MQ=255)
                                                                      ccGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcaacgtaagtaaagctgcgcacctgtttgccactacggtcatacaccacc     <  1:355073‑M2/100‑62 (MQ=255)
                                                                       nggggcCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccg    <  2:1095506‑M2/99‑63 (MQ=255)
                                                                        ggggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgc   <  1:1988918‑M2/100‑64 (MQ=255)
                                                                        ggggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgc   <  1:1521886‑M2/100‑64 (MQ=255)
                                                                        ggggCCGGATGGTGGCGCACAGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgc   <  1:1584364‑M2/100‑64 (MQ=255)
                                                                         gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgggcacctgtttgccactacggtcatacaccaccgcc  <  2:260616‑M2/100‑65 (MQ=255)
                                                                         gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgcc  <  2:969637‑M2/100‑65 (MQ=255)
                                                                         gggCCGGATGGTGGCGCACCGTCACACGGCCCGGCCccggtatttatcatcgtaagtaacgctgcgcacctgtttgccactacgctcatacaccaccgc   <  2:1043930‑M2/99‑64 (MQ=255)
                                                                         |                                                                                                   
GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTA  >  NC_000913/3620147‑3620319

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