New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3621534NA (NA)5 (0.040) 4/162 NT 100% coding (2343/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 = 3621592 0 (0.000)coding (2401/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

GCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCATGAAACTCGGCGACACACCGCTGGTGGAG  >  NC_000913/3621486‑3621633
                                                |                                                                                                   
gCTTTGTCAGCATGAGACCAGACATGCGTACAACGCGCAGTGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGC                                                  >  2:2703956‑M2/1‑100 (MQ=255)
gCTCTGGCAGCATGAGACCATACATTCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGC                                                  >  1:1453208‑M2/1‑100 (MQ=255)
gCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGC                                                  >  2:2480983‑M2/1‑100 (MQ=255)
gCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGC                                                  <  1:3830484‑M2/100‑1 (MQ=255)
gCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGC                                                  <  2:358058‑M2/100‑1 (MQ=255)
gCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGC                                                  >  1:1833595‑M2/1‑100 (MQ=255)
gCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGC                                                  >  1:162044‑M2/1‑100 (MQ=255)
gCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAACCTGCCCGCCGTGGAATGGCTGACCTACGGC                                                  <  2:2107299‑M2/100‑1 (MQ=255)
gCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACAGCTCTATACCGGACCGCCTGCCCGCCGTGGAATGGCTGACCTACGGC                                                  >  2:3275135‑M2/1‑100 (MQ=255)
gCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGACTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGC                                                  <  1:2761473‑M2/100‑1 (MQ=255)
 ctctGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTCGAATGGCTGACCTACGGC                                                  <  2:1395869‑M2/99‑1 (MQ=255)
  tctGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAG                                                >  2:598256‑M2/1‑100 (MQ=255)
  tctGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAG                                                >  1:3835226‑M2/1‑100 (MQ=255)
   ctGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGC                                               <  1:2390445‑M2/100‑1 (MQ=255)
   ctGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGC                                               >  1:85840‑M2/1‑100 (MQ=255)
   ctGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGC                                               >  2:2580944‑M2/1‑100 (MQ=255)
    ttgcagcatgAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCG                                              >  1:2113409‑M2/3‑100 (MQ=255)
    tGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCG                                              >  2:3409557‑M2/1‑100 (MQ=255)
    tGGCAGCATGAGACCAGACATGCGTACAACGCGAAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCTCGCCGTGGAATGGCTGACCTACGGCAGCG                                              <  2:2348358‑M2/100‑1 (MQ=255)
     ggCAGCATGAGACCAGACATGCGTATAACGCGCAGTGTTTGGCGAATCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                             >  1:337282‑M2/1‑100 (MQ=255)
     ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                             >  1:1406859‑M2/1‑100 (MQ=255)
     ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                             <  2:892440‑M2/100‑1 (MQ=255)
     ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                             >  1:917521‑M2/1‑100 (MQ=255)
     ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                             >  1:841741‑M2/1‑100 (MQ=255)
     ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                             >  1:3092163‑M2/1‑100 (MQ=255)
     ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                             <  1:3063701‑M2/100‑1 (MQ=255)
     ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCAGTGGAATGGCTGACCTACGGCAGCGG                                             <  2:996787‑M2/100‑1 (MQ=255)
     ggCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTATATACCGGACAGCCTACCCGCCGTGGAATGGCTGACCTACGGCAGCGG                                             <  1:457828‑M2/100‑1 (MQ=255)
      gcagcaTGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGT                                            >  2:3303413‑M2/1‑100 (MQ=255)
      gcagcaTGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGT                                            <  2:3804969‑M2/100‑1 (MQ=255)
       cagcaTGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTT                                           <  2:1917081‑M2/100‑1 (MQ=255)
       cagcaTGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCATGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTT                                           <  2:2522360‑M2/100‑1 (MQ=255)
                                              gCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTcgccagcccctgcgcgttgtacgcatgtctggtctcatg    <  1:1333467‑M2/100‑40 (MQ=255)
                                               cGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTcgccagcccctgcgcgttgtacgcatgtctggtctcatgc   <  2:2377482‑M2/100‑41 (MQ=255)
                                               cGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTcgccagcacctgagcgatgtacgcatgtctggtctcatgc   <  1:1464265‑M2/100‑41 (MQ=255)
                                                gAACCGCTGTATACCGGACAGCATGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTGcgccagcccctgagcgttgtacgcatgtgtggtcgcatgct  <  2:535118‑M2/100‑42 (MQ=255)
                                                |                                                                                                   
GCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCATGAAACTCGGCGACACACCGCTGGTGGAG  >  NC_000913/3621486‑3621633

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