New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3621014 =25 (0.720)8 (0.240) 5/280 2.4 19.4% coding (1823/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 3621066 = 43 (1.310)coding (1875/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Rejected: Frequency below cutoff threshold.

CGGAAGGTGGTGTGGCTGCCGTTTTCACTGGTCAGGCGGATGACCCGTCCGGCAGCATCGTATTCCATACTGCGCGTCAGCCCGCCCTGCGTGGTACGGACGGCCTTTCCCCACGCATCGTACTGTGTCCCGTTTCTGCTGCCGTCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3621161‑3621014
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgtccgTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA  >  NC_000913/3621066‑3621215
                                                                                                                                                                                                                                                                                                          
CGGAAGGTGGTGTGGCTGCCGTTTTCACTGGTCAGGCGGATGACCCGTCCGGCAGCATCGTATTCCATACTGCGCGTTAGCCCGCCCTGCGTGGTACGGACGGCCTTTCCCCACGCATCGTACTGTGTCCCGTTTCTGCTGCCGTCCGTAC                                                                                                                                                     <  1:70021/151‑1
          TGTGGCTGCCGTTTTCACTGGTCAGGCGGATGACCCGTCCGGCAGCATCGTATTCCATACTGCGCGTCAGCCCGCCCTGCGTGGTACGGACGGCCTTTCCCCACGCATCGTACTGTGTCCCGTTTCTGCTGCCGTCCGTACCACGCAGGGC                                                                                                                                           <  2:402593/151‑1
                                   GCGGATGACCCGTCCGGCAGCATCGTATTCCATACTGCGCGTCAGCCCGCCCTGCGTGGTACGGACGGCCTTTCCCCACGCATCGTACTGTGTCCCGTTTCTGCTGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAAT                                                                                                                     <  2:331599/148‑1
                                      GATGACCCGTCCGGCAGCATCGTATTCCATACTGCGCGTCAGCCCGCCCTGCGTGGTACGGACGGCCTTTCCCCACGCATCGTACTGTGTCCCGTTTCTGCTGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGAT                                                                                                                <  2:262736/150‑1
                                                                                            GGTACGGACGGCCTTTCCCCACGCATCGTACTGTGTCCCGTTTCTGCTGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTC                                                          >  2:40287/1‑150
                                                                                                                    ATCGTACTGTGTCCCGTTTCTGCTGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCCGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGA                                 >  1:422418/1‑151
                                                                                                                       GTACTGTGTCCCGTTTCTGCTGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATA                               >  1:315637/1‑150
                                                                                                                            GTGTCCCGTTTCTGCTGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAA                         >  1:374067/1‑151
                                                                                                                                            GCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGC          >  1:2757/1‑150
                                                                                                                                              CGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCG        >  1:534388/1‑150
                                                                                                                                              CGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCG        >  1:73236/1‑150
                                                                                                                                              CGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCG        >  2:89372/1‑150
                                                                                                                                                TCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCA      >  1:219803/1‑150
                                                                                                                                                  CGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACAC   <  1:482203/151‑1
                                                                                                                                                   GTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA  >  2:62110/1‑151
                                                                                                                                                                                                                                                                                                          
CGGAAGGTGGTGTGGCTGCCGTTTTCACTGGTCAGGCGGATGACCCGTCCGGCAGCATCGTATTCCATACTGCGCGTCAGCCCGCCCTGCGTGGTACGGACGGCCTTTCCCCACGCATCGTACTGTGTCCCGTTTCTGCTGCCGTCCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3621161‑3621014
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gccgtccgTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACA  >  NC_000913/3621066‑3621215

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