New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 52601048 (0.120)4 (0.010) 4/316 19.8 8.5% coding (2750/4281 nt) rhsD Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 = 526022 40 (0.100)coding (2762/4281 nt) rhsD Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Rejected: Position hash score below cutoff.

GACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/525883‑526010
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccggcAGCCGGTTGCCCGCCGGGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCCTGCCCGTGGCGCTGTAGCCGTATTCCCGCGTCTGTCGCGGGCCGCTGATGC  <  NC_000913/526022‑525893
                                                                                                                                                                                                                                                                  
GACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACC                                                                                                                                        <  1:860184/124‑1
GACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACC                                                                                                                                        >  2:860027/1‑124
                                                             GGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACC                                                                                                                                        >  1:1912398/1‑63
                                                             GGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACC                                                                                                                                        <  2:1912061/63‑1
                                                                               CCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGCCCCGGCAGCCGGTTGCCCGCCGGGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCCTGCCCGTGGCGCTGTAGCCGTATTCCCGCGTCTGTCGCGGGCCGCTGATGC  <  1:3965006/179‑1
                                                                               CCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCAGCCGGTTGCCCGCCGGGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCCTGCCCGTGGCGCTGTAGCCGTATTCCCGCGTCTGTCGCGGGCCGCTGATGC  >  2:3964347/1‑179
                                                                                          GACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCAGCCGGTTGCCCGCCGGGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCC                                                     >  1:3867124/1‑117
                                                                                          GACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCAGCCGGTTGCCCGCCGGGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCC                                                     <  2:3866493/117‑1
                                                                                                                             GGCGGCCGGTTGCCCGCCGGGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCCTGCCCGTGGCGCTGTAGCCGTATTCCCGCGTCTGTCGCGGGC           <  1:5132034/124‑1
                                                                                                                             GGCAGCCGGTTGCCCGCCGGGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCCTGCCCGTGGCGCTGTAGCCGTATTCCCGCGTCTGTCGCGGGC           >  2:5131161/1‑124
                                                                                                                              GCAGCCGGTTGCCCGCCGGGTCCGTGGCATACGGGATGCGGA                                                                                            <  1:1058509/42‑1
                                                                                                                              GCAGCCGGTTGCCCGCCGGGTCCGTGGCATACGGGATGCGGA                                                                                            >  2:1058301/1‑42
                                                                                                                                                                                                                                                                  
GACCTGGTGCGCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/525883‑526010
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccggcAGCCGGTTGCCCGCCGGGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCCTGCCCGTGGCGCTGTAGCCGTATTCCCGCGTCTGTCGCGGGCCGCTGATGC  <  NC_000913/526022‑525893

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