New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3622174NA (NA)7 (0.090) 3/168 NT 100% coding (2983/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 = 3622205 0 (0.000)coding (3014/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

AGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACA  >  NC_000913/3622106‑3622304
                                                                                                   |                                                                                                   
aGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGT                                                                                                     <  2:28049‑M2/100‑1 (MQ=255)
aGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGT                                                                                                     >  1:2223953‑M2/1‑100 (MQ=255)
aGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGT                                                                                                     >  1:847394‑M2/1‑100 (MQ=255)
aGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACCGGT                                                                                                     <  2:2024145‑M2/100‑1 (MQ=255)
                                                              cTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcggctttcgaccagcggctctgcatattgtgtccgcgtgtagtgcaccagc                                       <  2:1233619‑M2/100‑63 (MQ=255)
                                                              cTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcggctttcgaccagcggctctgcatattgtgtccgcgtgtagtgcaccagc                                       <  1:786579‑M2/100‑63 (MQ=255)
                                                               tttACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagcc                                      <  2:1597390‑M2/100‑64 (MQ=255)
                                                                 tACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccgg                                    <  2:634135‑M2/100‑66 (MQ=255)
                                                                 tACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccgg                                    <  1:528419‑M2/100‑66 (MQ=255)
                                                                  aCGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgctccagccggt                                   <  2:2249739‑M2/100‑67 (MQ=255)
                                                                    gACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcggctttcgaccagcggctctgcatattgtgtccgcgtgtagtgcaccagccggtgc                                 <  1:1474294‑M2/100‑69 (MQ=255)
                                                                                               cGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtc      >  1:2301021‑M2/1‑5 (MQ=255)
                                                                                               cGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACg       <  1:690852/99‑1 (MQ=11)
                                                                                                gggTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtcc     >  2:1225930‑M2/1‑4 (MQ=255)
                                                                                                gggTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGAt     <  1:1766229/100‑1 (MQ=32)
                                                                                                 ggTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccg    <  2:1531525‑M2/100‑98 (MQ=255)
                                                                                                  gtgtGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGaaa    >  2:1647723/1‑97 (MQ=14)
                                                                                                  gtgtGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATAc   >  1:2219503/1‑100 (MQ=32)
                                                                                                  gtgtGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATAc   <  2:1584690/100‑1 (MQ=32)
                                                                                                  gtgtGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATAc   >  1:1860770/1‑100 (MQ=32)
                                                                                                  gTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgc   >  2:498226‑M2/1‑2 (MQ=255)
                                                                                                   tgtgGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACa  <  2:442120/100‑1 (MQ=33)
                                                                                                   tgtgGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACa  >  1:1264713/1‑100 (MQ=32)
                                                                                                   tcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggtgctgactgtcgtaatggtaccggtgagtccgct  >  2:632890‑M2/1‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTACGGCTGGGACGGCGACCGGCTGACCACGATACA  >  NC_000913/3622106‑3622304

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.