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.

GGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCG  >  NC_000913/3620864‑3621161
                                                                                                                                                      |                                                                                                                                                   
tgtgtgcggccgtcaaagccggtttcctgtatcagccggtcaagtacatcgtaacggaaggtggtgtggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgtcagcccgccctgcgtggtac                                                                                                                                                     >  2:62110‑M1/151‑151 (MQ=255)
 gtgtgcggccgtcaaagccggtttcctgtatcagccggtcaagtacatcgtaacggaaggtggtgtggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgtcagcccgccctgcgtggtacg                                                                                                                                                    <  1:482203‑M1/2‑1 (MQ=255)
    tgcggccgtcaaagccggtttcctgtatcagccggtcaagtacatcgtaacggaaggtggtgtggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgtcagcccgccctgcgtggtacgga                                                                                                                                                  >  1:219803‑M1/147‑150 (MQ=255)
      cggccgtcaaagccggtttcctgtatcagccggtcaagtacatcgtaacggaaggtggtgtggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgtcagcccgccctgcgtggtacggacg                                                                                                                                                >  2:89372‑M1/145‑150 (MQ=255)
      cggccgtcaaagccggtttcctgtatcagccggtcaagtacatcgtaacggaaggtggtgtggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgtcagcccgccctgcgtggtacggacg                                                                                                                                                >  1:73236‑M1/145‑150 (MQ=255)
      cggccgtcaaagccggtttcctgtatcagccggtcaagtacatcgtaacggaaggtggtgtggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgtcagcccgccctgcgtggtacggacg                                                                                                                                                >  1:534388‑M1/145‑150 (MQ=255)
        gccgtcaaagccggtttcctgtatcagccggtcaagtacatcgtaacggaaggtggtgtggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgtcagcccgccctgcgtggtaCGGACGgc                                                                                                                                              >  1:2757‑M1/143‑150 (MQ=255)
                       ttcctgtatcagccggtcaagtacatcgtaacggaaggtggtgtggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgtcagcccgccctgcgtggtaCGGACGGCAGCAGAAACGGGacac                                                                                                                              >  1:374067‑M1/128‑151 (MQ=255)
                             tatcagccggtcaagtacatcgtaacggaaggtggtgtggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgttagcccgccctgcgtggtaCGGACGGCAGCAGAAACGGGACACAGTAc                                                                                                                         >  1:315637‑M1/122‑150 (MQ=255)
                               tcagccggtcaagtacatcgtaacggaaggtggtgtggctgccgttttcacgggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgttagcccgccctgcgtggtaCGGACGGCAGCAGAAACGGGACACAGTACGAt                                                                                                                      >  1:422418‑M1/120‑151 (MQ=255)
                                                        gaaggtggtgtggctgccgttttcactggtcaggcggatgacccgtccggcagcatcgtattccatactgcgcgtcagcccgccctgcgtggtaCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTAcc                                                                                              >  2:40287‑M1/95‑150 (MQ=255)
                                                                           ccGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAcgc                                                                           <  1:200824/150‑1 (MQ=31)
                                                                                aaaGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTa                                                                     >  1:301877/1‑151 (MQ=31)
                                                                                    acacGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGa                                                                 <  1:38443/151‑1 (MQ=31)
                                                                                    acacGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGa                                                                 <  1:517889/151‑1 (MQ=31)
                                                                                      ccGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAAt                                                               <  1:358857/150‑1 (MQ=14)
                                                                                      acGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAAt                                                               <  2:278024/151‑1 (MQ=31)
                                                                                      acGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAAt                                                               <  1:172651/151‑1 (MQ=31)
                                                                                      acGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAAt                                                               <  1:424381/151‑1 (MQ=31)
                                                                                           gggCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGAt                                                          >  2:4989/1‑151 (MQ=31)
                                                                                            ggCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGAt                                                          <  2:353750/150‑1 (MQ=31)
                                                                                            ggCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGAt                                                          <  1:121606/150‑1 (MQ=31)
                                                                                              ccATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGAtgc                                                        <  2:201123/150‑1 (MQ=31)
                                                                                                     aCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGAc                                                <  2:331045/151‑1 (MQ=31)
                                                                                                        cgGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACggg                                             >  2:469887/1‑151 (MQ=31)
                                                                                                              gAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATc                                        >  1:123872/1‑150 (MQ=31)
                                                                                                              atcgtattccatactgcgcgtcagcccgccctgcgtggtaCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATc                                        <  2:262736‑M1/110‑1 (MQ=255)
                                                                                                                aTACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCg                                      >  2:350842/1‑150 (MQ=31)
                                                                                                                   attccatactgcgcgtcagcccgccctgcgtggtaCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGc                                     <  2:331599‑M1/113‑1 (MQ=255)
                                                                                                                     acaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGa                                 >  1:496186/1‑150 (MQ=31)
                                                                                                                        tCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAg                             >  2:380864/1‑151 (MQ=31)
                                                                                                                           ccGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTg                           <  1:220231/150‑1 (MQ=31)
                                                                                                                           ccGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGa                          >  2:249595/1‑151 (MQ=31)
                                                                                                                                gacctgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACg                     >  2:329927/1‑151 (MQ=31)
                                                                                                                                      accGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGcc               <  2:262783/151‑1 (MQ=31)
                                                                                                                                      accGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGcc               <  2:201784/151‑1 (MQ=31)
                                                                                                                                         gccctgcgtggtaCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCcaca            <  2:402593‑M1/138‑1 (MQ=255)
                                                                                                                                             tCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcaccac        <  1:186003/151‑1 (MQ=31)
                                                                                                                                                   gtaCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTAACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCg  <  1:70021‑M1/148‑1 (MQ=255)
                                                                                                                                                   ccccGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTcc   <  2:338269/150‑1 (MQ=31)
                                                                                                                                                      |                                                                                                                                                   
GGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCG  >  NC_000913/3620864‑3621161

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