New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3621553NA (NA)4 (0.040) 4/244 NT 100% coding (2362/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 = 3621603 0 (0.000)coding (2412/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

AGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCC  >  NC_000913/3621449‑3621679
                                                                                                        |                                                                                                                              
aGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTG                                                                                    >  4:172835‑M2/1‑149 (MQ=255)
aGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTG                                                                                    >  3:495171‑M2/1‑149 (MQ=255)
 gCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTG                                                                                    >  1:564389‑M2/1‑148 (MQ=255)
   gTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCA                                                                                 >  4:742894‑M2/1‑149 (MQ=255)
   gTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGC                                                                                  >  1:178634‑M2/1‑148 (MQ=255)
   gTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGC                                                                                  >  1:17811‑M2/1‑148 (MQ=255)
    tCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                                                <  4:166731‑M2/149‑1 (MQ=255)
    tCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                                                >  3:226967‑M2/1‑149 (MQ=255)
    tCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                                                <  3:460842‑M2/149‑1 (MQ=255)
    tCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCA                                                                                 <  1:332908‑M2/148‑1 (MQ=255)
     cAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                                                <  2:164546‑M2/148‑1 (MQ=255)
     cAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                                                <  2:501117‑M2/148‑1 (MQ=255)
     cAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                                                <  2:618816‑M2/148‑1 (MQ=255)
      aGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGC                                                                              >  4:8910‑M2/1‑149 (MQ=255)
      aGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGC                                                                              >  3:709181‑M2/1‑149 (MQ=255)
       gACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGC                                                                              >  2:87368‑M2/1‑148 (MQ=255)
       gACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGC                                                                              <  1:544460‑M2/148‑1 (MQ=255)
              cATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                                                >  3:660607‑M2/1‑139 (MQ=255)
              cATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                                                <  4:660607‑M2/139‑1 (MQ=255)
                             gAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCA                                                                                 >  2:402308‑M2/1‑123 (MQ=255)
                             gAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCA                                                                                 <  1:402308‑M2/123‑1 (MQ=255)
                                                                                                 aCCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCtgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagt       >  3:71418‑M2/1‑58 (MQ=255)
                                                                                                 aCCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCtgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagt       <  4:71418‑M2/129‑72 (MQ=255)
                                                                                                      aCAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCtgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagtgcttc  <  2:360269‑M2/129‑77 (MQ=255)
                                                                                                      aCAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCtgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagtgcttc  >  1:360269‑M2/1‑53 (MQ=255)
                                                                                                        |                                                                                                                              
AGCGTCAGACGGTGCATCACCCGCAGACGGAAGCACTGCTCTGGCAGCATGAGACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCC  >  NC_000913/3621449‑3621679

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