New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 732888 =34 (0.710)4 (0.090) 4/274 6.7 10.5% coding (3306/4194 nt) rhsC Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 734728 = 35 (0.750)pseudogene (596/1521 nt) rhsO pseudogene, Rhs family
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGGTCGCAGTGGTACAGGTGAATTTTTCGCGCCGGCGTGTATACCGGGTCCATCTGGCTTTGCATCTGCGCCACCGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/733031‑732888
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttccCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAA  >  NC_000913/734728‑734795
                                                                                                                                                                                                                    
TGGTCGCAGTGGTACAGGTGACTTTTGCGCGCCGGCGTGTAGACCGGGTCCATCGGGCTTTGCAGCTGCGCCACCGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTC                                                                       <  2:33600/143‑1
           GTACAGGTGAATTTTTCGCGCCGGCGTGTATACCGGGTCCATCTGGCTTTGCATCTGCGCCACCGTCAGGCCGCACGATGCCAGCCCGCGGCGGCCCTCTTCACGCACCCGGTCGGCGAGGAGTTGACTTTCGCGCCGGgcgg                                                            >  1:368265/1‑139
                                    gcgtcTACCGGGTCCATCTGGCTTTGCATCTGCGCCACCGTCAGGCCGCGCTGTGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAA                                   <  2:9224/138‑1
                                     TGTATACCGGGTCCATCTGGCTTTGCATCTGCGCCACCGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAA                                  <  1:38977/143‑1
                                                                     tctactGTCAGGCCGCACGATTCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAA  <  1:374699/137‑1
                                                                                                                                                                                                                    
TGGTCGCAGTGGTACAGGTGAATTTTTCGCGCCGGCGTGTATACCGGGTCCATCTGGCTTTGCATCTGCGCCACCGTCAGGCCGCACGATGCCAGCCAGCGGCGGCTTTCCTCACTCACCCGGTCAGCCAGGATTTCACTTTCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/733031‑732888
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ttccCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAA  >  NC_000913/734728‑734795

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