New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3620220NA (NA)11 (0.110) 7/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
                                                                         |                                                                                                   
GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCacac                                                                           <  2:2614539‑M2/100‑1 (MQ=255)
 CGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACg                                                                          <  2:1062491‑M2/100‑1 (MQ=255)
   GCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACacgtg                                                                        >  1:2735520‑M2/1‑98 (MQ=255)
   GCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACggg                                                                        <  1:433786‑M2/100‑1 (MQ=255)
   GCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACggg                                                                        <  1:450381‑M2/100‑1 (MQ=255)
   GCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACggg                                                                        <  2:1089265‑M2/100‑1 (MQ=255)
   GCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTGACTTACGATGATAAATACCGGGGCCGGATTGTGGCGCACCGTCACACggg                                                                        <  1:1345831‑M2/100‑1 (MQ=255)
    cGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGC                                                                       <  2:2702629‑M2/100‑1 (MQ=255)
    cGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGC                                                                       <  1:47581‑M2/100‑1 (MQ=255)
     gAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCC                                                                      >  1:743374‑M2/1‑100 (MQ=255)
      aaCTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCG                                                                     >  2:455527‑M2/1‑100 (MQ=255)
       aCTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGG                                                                    <  1:363990‑M2/100‑1 (MQ=255)
       aCTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGG                                                                    >  1:479648‑M2/1‑100 (MQ=255)
       aCTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGG                                                                    >  2:229015‑M2/1‑100 (MQ=255)
       aCTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCGTTACTTACGATGATAAATACCGGGGCCCGATCGTGGCGCACCGTCAAACGGGCCG                                                                     <  1:1724345‑M2/99‑1 (MQ=255)
         tagcgGTGGTGTATGACCGTAGTGCCAAACAGGCGCGCAGCTTGACTTACGATGATAAATACAGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGG                                                                    >  2:829021‑M2/3‑98 (MQ=255)
         tGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTCCTTACGATGATAAATACCGGGGTCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                  >  2:2770213‑M2/1‑100 (MQ=255)
         tGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                  >  2:1452352‑M2/1‑100 (MQ=255)
         tGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                  >  1:744724‑M2/1‑100 (MQ=255)
         tGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC                                                                  >  2:280676‑M2/1‑100 (MQ=255)
                                                               atAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcata            <  2:2082182‑M2/100‑55 (MQ=255)
                                                                    tACCGGGGCCGGATGGTGGCGCACCGTCACTCGGGCCGGCCccggtgcttatcatcgtaagtaaacctgcgaacctgtttgccactacggtcatacacca       <  1:2530126‑M2/100‑60 (MQ=255)
                                                                     aCCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccac      <  2:2769316‑M2/100‑61 (MQ=255)
                                                                      ccGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaggtcaagctgcgcacctgtttgccactacggtcatacaccacc     <  1:2043878‑M2/100‑62 (MQ=255)
                                                                        ggggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttctcatccgcagtaaagctgcgcacctgtttgccactacggtcatacacccccg    <  1:2416062‑M2/99‑63 (MQ=255)
                                                                        ggggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgc   <  2:2332322‑M2/100‑64 (MQ=255)
                                                                        ggggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgc   <  2:2484532‑M2/100‑64 (MQ=255)
                                                                         gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacacccccgcc  <  1:1037179‑M2/100‑65 (MQ=255)
                                                                         gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgcc  <  2:1085244‑M2/100‑65 (MQ=255)
                                                                         gggCCGGATGGTGGCGCACCGTCACACGGGCCGGCCccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgcc  <  2:2397692‑M2/100‑65 (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.