New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3620123 =NA (NA)3 (0.030) 3/166 NT 100% coding (932/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 3620147 = 0 (0.000)coding (956/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

CAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACAC  >  NC_000913/3620049‑3620246
                                                                                                  |                                                                                                   
ccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagc                                                                                                    <  1:1974165‑M2/2‑1 (MQ=255)
ccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagc                                                                                                    <  2:1310794‑M2/2‑1 (MQ=255)
 cggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcg                                                                                                   >  1:1279181‑M2/98‑100 (MQ=255)
 cggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcg                                                                                                   <  1:139830‑M2/3‑1 (MQ=255)
  ggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcgg                                                                                                  >  1:1229424‑M2/97‑100 (MQ=255)
  ggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcgg                                                                                                  <  2:629590‑M2/4‑1 (MQ=255)
  ggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcgg                                                                                                  >  1:2379941‑M2/97‑100 (MQ=255)
  ggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcatggcgtccagccatagcgcaccagcgg                                                                                                  >  1:169667‑M2/97‑100 (MQ=255)
   gtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggc                                                                                                 <  2:124942‑M2/5‑1 (MQ=255)
   gtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggc                                                                                                 <  2:2134677‑M2/5‑1 (MQ=255)
   gtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagtacgccgcgtggcgtccagccatagcgcaccagcggc                                                                                                 >  2:2043739‑M2/96‑100 (MQ=255)
    tatttattatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcg                                                                                                >  1:1981621‑M2/95‑100 (MQ=255)
    tatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcg                                                                                                <  2:2464724‑M2/6‑1 (MQ=255)
      tttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagatcggaa                                                                                              >  2:1692554‑M2/93‑93 (MQ=255)
                                                      caccaccgccagttcgccgcgtggcgtccagccatagcgcaccaGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTttact                                              >  1:2490804‑M2/45‑100 (MQ=255)
                                                                gccgcgtggcgtccagccatagcgcaccagcggcGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATa                                    >  1:1115102‑M2/35‑100 (MQ=255)
                                                                 gttcgccgcgtggcgtccagccatagcgcaccaGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATaa                                   >  2:585426‑M2/34‑100 (MQ=255)
                                                                                                  GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCGAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGAGGGTGGCGCACCGTcaca   >  1:2978786‑M2/1‑99 (MQ=255)
                                                                                                  GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCACCCCGGTGCGCAGCTTTACTTACGATGATAACTCCCGGGGCCGGATGGTGGCGCACCGTCacac  <  2:902953‑M2/100‑1 (MQ=255)
                                                                                                  GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCacac  >  2:1986283‑M2/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACAC  >  NC_000913/3620049‑3620246

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