New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3620220NA (NA)15 (0.120) 6/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

GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTT  >  NC_000913/3620147‑3620318
                                                                         |                                                                                                  
GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCacac                                                                          >  2:1627102‑M2/1‑100 (MQ=255)
GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCacac                                                                          <  1:1422475‑M2/100‑1 (MQ=255)
GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCacac                                                                          <  2:1777014‑M2/100‑1 (MQ=255)
 CGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACg                                                                         >  2:1694747‑M2/1‑100 (MQ=255)
 CGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACGTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTcaca                                                                           <  1:2198669‑M2/98‑1 (MQ=255)
  GGCGAACTGGCGGTGGTGTATGACTGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGtcacacag                                                                        >  2:1910783‑M2/1‑98 (MQ=255)
  GGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACgg                                                                        >  1:3289264‑M2/1‑100 (MQ=255)
  GGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACgg                                                                        <  2:2069191‑M2/100‑1 (MQ=255)
   GCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACggg                                                                       <  1:2331289‑M2/100‑1 (MQ=255)
   GCGAACTGGCGGTGGTGTATAACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACggg                                                                       >  1:2261104‑M2/1‑100 (MQ=255)
    cGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGC                                                                      >  1:1255288‑M2/1‑100 (MQ=255)
    cGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGC                                                                      >  2:2421694‑M2/1‑100 (MQ=255)
     gAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCC                                                                     >  2:2070225‑M2/1‑100 (MQ=255)
     gAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCC                                                                     >  1:3060541‑M2/1‑100 (MQ=255)
     gAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCC                                                                     >  1:3258799‑M2/1‑100 (MQ=255)
     gAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCC                                                                     >  2:1017804‑M2/1‑100 (MQ=255)
      aaCTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCG                                                                    <  2:1732855‑M2/100‑1 (MQ=255)
      aaCTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCG                                                                    >  2:2094818‑M2/1‑100 (MQ=255)
      aaCTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCC                                                                     >  2:1222133‑M2/1‑99 (MQ=255)
       aCTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGG                                                                   <  1:751856‑M2/100‑1 (MQ=255)
        cTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                  <  2:2997224‑M2/100‑1 (MQ=255)
        cTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                  >  2:1893246‑M2/1‑100 (MQ=255)
        cTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                  >  2:1790828‑M2/1‑100 (MQ=255)
        cTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                  <  1:301042‑M2/100‑1 (MQ=255)
        cTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                  <  1:3303537‑M2/100‑1 (MQ=255)
        cTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                  <  1:3670649‑M2/100‑1 (MQ=255)
        cTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC                                                                  <  2:1442046‑M2/100‑1 (MQ=255)
         tGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                 >  2:3554735‑M2/1‑100 (MQ=255)
         tGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                 >  2:2305571‑M2/1‑100 (MQ=255)
         tGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                 >  1:1694740‑M2/1‑100 (MQ=255)
          ggCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGGAGCTTTGCTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                 <  1:2583229‑M2/99‑1 (MQ=255)
                                                                  aaTACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGACccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacac        <  1:2019525‑M2/100‑58 (MQ=255)
                                                                     aCCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcctacaccac     <  1:2209382‑M2/100‑61 (MQ=255)
                                                                     aCCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccac     <  1:2403779‑M2/100‑61 (MQ=255)
                                                                     aCCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccac     <  2:2579434‑M2/100‑61 (MQ=255)
                                                                      ccGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccacc    <  2:395493‑M2/100‑62 (MQ=255)
                                                                      ccGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccacc    <  2:996491‑M2/100‑62 (MQ=255)
                                                                      ccGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccacc    <  2:3125091‑M2/100‑62 (MQ=255)
                                                                      ccGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccacc    <  2:2304718‑M2/100‑62 (MQ=255)
                                                                        ggggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgttagtaaagctgcgcacctgtttgccactacggtcatacaccaccgc  <  1:1069182‑M2/100‑64 (MQ=255)
                                                                        ggggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgc  <  1:2509843‑M2/100‑64 (MQ=255)
                                                                         |                                                                                                  
GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTT  >  NC_000913/3620147‑3620318

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