New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3621556NA (NA)7 (0.080) 4/170 NT 100% coding (2365/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 = 3621604 0 (0.000)coding (2413/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

ACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACC  >  NC_000913/3621502‑3621655
                                                      |                                                                                                   
aCCAGACATGCGTACCACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                        >  1:2564218‑M2/1‑100 (MQ=255)
aCCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                        >  1:588153‑M2/1‑100 (MQ=255)
aCCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                        >  2:836835‑M2/1‑100 (MQ=255)
aCCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                        >  2:1426917‑M2/1‑100 (MQ=255)
aCCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAG                                                        >  1:2295891‑M2/1‑100 (MQ=255)
 ccAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGG                                                       <  1:2234049‑M2/100‑1 (MQ=255)
 ccAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGG                                                       <  1:256150‑M2/100‑1 (MQ=255)
  cAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACTGCAGCGGTTACCTGGCAGGC                                                      <  2:1551332‑M2/100‑1 (MQ=255)
  cAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGC                                                      >  2:1956685‑M2/1‑100 (MQ=255)
   aGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCA                                                     >  1:1524458‑M2/1‑100 (MQ=255)
   aGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCA                                                     >  1:1453262‑M2/1‑100 (MQ=255)
   aGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCA                                                     >  2:709108‑M2/1‑100 (MQ=255)
   aGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTATATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCA                                                     >  1:307661‑M2/1‑100 (MQ=255)
           gTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGGCCTACGGCAGCGGCTACCTGGCAGGCAggc                                                  <  2:1870768‑M2/95‑4 (MQ=255)
                                                  cAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCAggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgca      <  1:1507273‑M2/100‑48 (MQ=255)
                                                    gCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCAggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatg    <  2:163966‑M2/100‑50 (MQ=255)
                                                    gCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCAggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatg    <  2:1925438‑M2/100‑50 (MQ=255)
                                                    gCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCAggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatg    <  1:1351022‑M2/100‑50 (MQ=255)
                                                      cTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCAggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtc  <  1:1099983‑M2/100‑52 (MQ=255)
                                                      cTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCAggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtc  <  2:2567592‑M2/100‑52 (MQ=255)
                                                      |                                                                                                   
ACCAGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACC  >  NC_000913/3621502‑3621655

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