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

AGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTA  >  NC_000913/3622102‑3622271
                                                                        |                                                                                                 
gtagtgcaccaGCCGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAAC                                                                        <  1:1529545‑M2/95‑1 (MQ=255)
aGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAAC                                                                        <  2:2170375‑M2/100‑1 (MQ=255)
aGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAAC                                                                        <  2:424261‑M2/100‑1 (MQ=255)
 tagtgcaccAGCCGGTGCACTACACGCGGACATAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCACGGTGGCAAAACG                                                                       <  1:2534804‑M2/96‑1 (MQ=255)
 ggcagCACCGTCTGGTGCACTCCACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTTGGCCGCAGGGTGGCAAAACG                                                                       >  1:2748799‑M2/3‑100 (MQ=255)
 gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACG                                                                       <  2:3573496‑M2/100‑1 (MQ=255)
 gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACG                                                                       >  2:2999346‑M2/1‑100 (MQ=255)
 gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACG                                                                       <  2:3598317‑M2/100‑1 (MQ=255)
 gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACG                                                                       <  2:2356813‑M2/100‑1 (MQ=255)
 gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACG                                                                       >  2:2118554‑M2/1‑100 (MQ=255)
 gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACG                                                                       <  2:1848997‑M2/100‑1 (MQ=255)
 gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACG                                                                       <  2:164417‑M2/100‑1 (MQ=255)
 gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACG                                                                       <  1:569804‑M2/100‑1 (MQ=255)
 gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACG                                                                       >  1:113342‑M2/1‑100 (MQ=255)
 gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACG                                                                       <  1:2586636‑M2/100‑1 (MQ=255)
  tCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGG                                                                      <  1:3775998‑M2/100‑1 (MQ=255)
  tCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGG                                                                      >  2:1382633‑M2/1‑100 (MQ=255)
  tCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGG                                                                      >  2:598877‑M2/1‑100 (MQ=255)
  agtgcacCAGCCGGTGCACTACACGCGGACGCAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGG                                                                      <  2:1299592‑M2/97‑1 (MQ=255)
  agtgcacCAGCCGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGG                                                                      <  2:2303768‑M2/97‑1 (MQ=255)
   cAGCACCGGCTGGTGCACTACACGCGGCCACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGG                                                                      >  1:3241633‑M2/1‑99 (MQ=255)
   cAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGCCCCGCTGGGCCGCAGGGTGGCTAAACGGG                                                                     >  1:388842‑M2/1‑100 (MQ=255)
   cAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGG                                                                     >  1:2164728‑M2/1‑100 (MQ=255)
   cAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGG                                                                     <  2:1976545‑M2/100‑1 (MQ=255)
   cAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTAGTCGGACCTCGCTATATTTACGACCCGCTGGGCCGCAGGGTGGCCAAACGG                                                                      <  1:3216403‑M2/99‑1 (MQ=255)
    aGCACCGGCTGGTGCACTACACGCGGACCCAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGT                                                                    <  1:270777‑M2/100‑1 (MQ=255)
    aGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGT                                                                    <  2:403707‑M2/100‑1 (MQ=255)
    aGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGT                                                                    <  1:3109758‑M2/100‑1 (MQ=255)
    aGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGT                                                                    <  1:1840067‑M2/100‑1 (MQ=255)
    aGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGT                                                                    <  1:1536569‑M2/100‑1 (MQ=255)
                                           agCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctcttcatatt                             <  2:1102977‑M2/100‑40 (MQ=255)
                                                                aTCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctctgcatattgtgtccgcgtgtagtgcacca        <  2:76649‑M2/100‑61 (MQ=255)
                                                                aTCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctctgcatattgtgtccgcgtgtagtgcacca        <  2:284328‑M2/100‑61 (MQ=255)
                                                                 tCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcacca        <  1:1958260‑M2/99‑61 (MQ=255)
                                                                   tttACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagcc     <  1:2444685‑M2/100‑64 (MQ=255)
                                                                    ttACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccg    <  1:2988400‑M2/100‑65 (MQ=255)
                                                                      aCGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTcgtaaagatagcgactttcgaccagcggctcttcatattgtgtccgcgtgtagtgcaccagccggt  <  1:989159‑M2/100‑67 (MQ=255)
                                                                        |                                                                                                 
AGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGAAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTGTGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCGCAAGTGACCTGGTA  >  NC_000913/3622102‑3622271

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.