New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3622151NA (NA)7 (0.070) 4/168 NT 100% coding (2960/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 = 3622177 0 (0.000)coding (2986/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

CGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCT  >  NC_000913/3622078‑3622276
                                                                                                   |                                                                                                   
atctctcacCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACC                                                                                                     >  1:762290‑M2/5‑100 (MQ=255)
    cTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCagat                                                                                                 <  2:762287‑M2/100‑5 (MQ=255)
                                                  gcgGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgact                                                   <  1:439935‑M2/100‑51 (MQ=255)
                                                                    gCCGCTGGTCGAAAGCCGCTATCTTTACGACCagcggctctgcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtg                                 <  1:775983‑M2/100‑69 (MQ=255)
                                                                     ccGCTGGTCGAAAGTCGCTATCTTTACGACCagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtga                                <  1:635539‑M2/100‑70 (MQ=255)
                                                                     ccGCTGGTCGAAAGTCGCTATCTTTACGACCagcggctctgcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgg                                <  1:2273326‑M2/100‑70 (MQ=255)
                                                                      cGCTGGTCGAAAGTCGCTATCTTTACGACCagcggctctgcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtggg                               <  1:1630199‑M2/100‑71 (MQ=255)
                                                                       gCTGGTCGAAAGTCGCTATCTTTACGACCagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagt                              <  1:808975‑M2/100‑72 (MQ=255)
                                                                        cTGGTCGAAAGTCGCTATCTTTACGACCagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtc                             <  2:1337720‑M2/100‑73 (MQ=255)
                                                                                               cGACCagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgctcatcatccgtgcggataac      >  1:665070‑M2/1‑5 (MQ=255)
                                                                                               cGACCagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgctcatcatccgtgcggataac      >  1:2318750‑M2/1‑5 (MQ=255)
                                                                                               cGACCagcggctcttcatagtgggtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgaggccgctcatcatccgtgcggataac      >  2:1720932‑M2/1‑5 (MQ=255)
                                                                                                gACCagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgctcatcatccgtgcggataacc     <  2:2412619‑M2/100‑97 (MQ=255)
                                                                                                 aCCagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgctcatcatccgtgcggataaccc    >  1:2417835‑M2/1‑3 (MQ=255)
                                                                                                 aCCagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgctcatcatccgtgcggataaccc    >  2:1763682‑M2/1‑3 (MQ=255)
                                                                                                  ccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgctcatcatccgtgcggataacccc   <  1:1740454‑M2/100‑99 (MQ=255)
                                                                                                  ccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgctcatcatccgtgcggataaccc    >  1:135862‑M2/1‑2 (MQ=255)
                                                                                                   cagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgctcatcatccgtgcggataacccct  <  1:2378341‑M2/100‑100 (MQ=255)
                                                                                                   cagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgctcatcatccgtgcggataacccct  >  2:1529680‑M2/1‑1 (MQ=255)
                                                                                                   |                                                                                                   
CGGACTCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCT  >  NC_000913/3622078‑3622276

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