New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 525435 =195 (1.940)3 (0.030) 3/270 6.8 2.0% coding (2175/4281 nt) rhsD Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 1527949 = 106 (1.080)pseudogene (48/2037 nt) rhsE pseudogene, Rhs family
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CGGCAACACGGTGGCCTTCGCTCAGGTGGCTGATGTCTGTCAGCCAGCCGTGGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/525490‑525435
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaCTACGGCTATGACGATAAAGGCCGCCTGACCGGCGAACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGAC  >  NC_000913/1527949‑1528087
                                                                                                                                                                                                   
CGGCAACACGGTGGCCTTCGCTCAGGTGGCTGATGTCTGTCAGCCAGCCGTGGCCACTACGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCT                                                                  >  1:193166/1‑131
CGGCAACACGGTGGCCTTCGCTCAGGTGGCTGATGTCTGTCAGCCAGCCGTGGCCACTACGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCT                                                                  <  2:193166/131‑1
                                                tcgGGCCACTACGGCTATGACGATAAAGGCCGCCTGACCGGCGAACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCAC         <  2:364586/137‑1
                                                      CACTACGGCTATGACGATAAAGGCCGCCTGACCGGCGAACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGA   <  1:1572969/140‑1
                                                      CACTACGGCTATGACGATAAAGGCCGCCTGACCGGCGAACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGA   <  2:1545365/140‑1
                                                      CACTACGGCTATGACGATAAAGGCCGCCTGACCGGCGAACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGA   <  2:1759217/140‑1
                                                      CACTACGGCTATGACGATAAAGGCCGCCTGACCGGCGAACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGA   <  2:751035/140‑1
                                                       ACTACGGCTATGACGATAAAGGCCGCCTGACCGGCGAACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGAC  <  2:272317/140‑1
                                                                                                                                                                                                   
CGGCAACACGGTGGCCTTCGCTCAGGTGGCTGATGTCTGTCAGCCAGCCGTGGCCA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/525490‑525435
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ccaCTACGGCTATGACGATAAAGGCCGCCTGACCGGCGAACGCCAGACGGTGGAGAACCCGGAGACGGGGGAACTGCTGTGGCACCATGAGACCGGACATGCGTACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGAC  >  NC_000913/1527949‑1528087

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