New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 52602740 (0.100)4 (0.010) 4/316 19.8 7.4% coding (2767/4281 nt) rhsD Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 = 1528459 61 (0.150)pseudogene (558/2037 nt) rhsE pseudogene, Rhs family
Rejected: Position hash score below cutoff.

GCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/525893‑526027
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccggGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCCTGCCCGTGGCGCTGTAGCCGTATTCCCGCGTCTGTCGCGGGCCGCTGATGCGCACCAGGTCGCCGTTGTCATTCCACCCGTAGTCACGGTCATACACCAGGCTGTTCAGGTGCTGGCTCTGTAACTGGCCTGCGGGGGTGTATGTGCTGGTCAGTTTATATGCGGCATTACTGCCTGCCATGCTGCCG  <  NC_000913/1528459‑1528210
                                                                                                                                                                                                                                                                                                                                                                                                 
GCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGG                                                                                                                                                                                                                <  2:3964347/179‑1
GCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGGGCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGG                                                                                                                                                                                                                >  1:3965006/1‑179
         GCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCC                                                                                                                                                                                                                                                              <  2:5131161/124‑1
         GCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCCGCC                                                                                                                                                                                                                                                              >  1:5132034/1‑124
                                                   GGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGGTCCGTGGCATACGGGATGCGGATGTCCAGGTC                                                                                                                                                                                                                           <  1:3867124/117‑1
                                                   GGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGGGTCCGTGGCATACGGGATGCGGATGTCCAGGTC                                                                                                                                                                                                                           >  2:3866493/1‑117
                                                                                          TCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGC                                                                                                                                                                                                                                                               >  1:1058509/1‑42
                                                                                          TCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGC                                                                                                                                                                                                                                                               <  2:1058301/42‑1
                                                                                                                                      GGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCCTGCCCGTGGCGCTGTAGCCGTATTCCCGCGTCTGTCGCGGGCCGCTGATGCGCACCAGGTCGCCGTTGTCATTCCACCCGTAGTCACGGTCATACACCAGGCTGTTCAGGTGCTGGCTCTGTAACTGGCCTGCGGGGGTGTATGTGCTGGTCCGTTTATATGCGGCATTACTGCCTGCCCTGCTGCCG  >  1:5935035/1‑251
                                                                                                                                      GGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCCTGCCCGTGGCGCTGTAGCCGTATTCCCGCGTCTGTCGCGGGCCGCTGATGCGCACCAGGTC                                                                                                                                 >  1:860184/1‑124
                                                                                                                                      GGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCCTGCCCGTGGCGCTGTAGCCGTATTCCCGCGTCTGTCGCGGGCCGCTGATGCGCACCAGGTC                                                                                                                                 <  2:860027/124‑1
                                                                                                                                      GGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCC                                                                                                                                                                                              <  1:1912398/63‑1
                                                                                                                                      GGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCC                                                                                                                                                                                              >  2:1912061/1‑63
                                                                                                                                                                                                                                                                                                                                                                                                 
GCATCAGCGGCCCGCGACAGACGCGGGAATACGGCTACAGCGCCACGGGCAGGCTGGAGAGTGTGCGCACCCTCGCACCAGACCTGGACATCCGCATCCCGTATGCCACGGACCCGGCGGGCAACCGGCTGCCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/525893‑526027
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccggGTCCGTGGCATACGGGATGCGGATGTCCAGGTCTGGTGCGAGGGTGCGCACACTCTCCAGCCTGCCCGTGGCGCTGTAGCCGTATTCCCGCGTCTGTCGCGGGCCGCTGATGCGCACCAGGTCGCCGTTGTCATTCCACCCGTAGTCACGGTCATACACCAGGCTGTTCAGGTGCTGGCTCTGTAACTGGCCTGCGGGGGTGTATGTGCTGGTCAGTTTATATGCGGCATTACTGCCTGCCATGCTGCCG  <  NC_000913/1528459‑1528210

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