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

AGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATA  >  NC_000913/3621505‑3621703
                                                                                                   |                                                                                                   
aGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCA                                                                                                     >  1:1453262‑M2/1‑100 (MQ=255)
aGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCA                                                                                                     >  1:1524458‑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                                                    <  1:1351022‑M2/100‑50 (MQ=255)
                                                 gCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCAggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatg                                                    <  2:1925438‑M2/100‑50 (MQ=255)
                                                 gCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCAggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatg                                                    <  2:163966‑M2/100‑50 (MQ=255)
                                                   cTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCAggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtc                                                  <  1:1099983‑M2/100‑52 (MQ=255)
                                                   cTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCAggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtc                                                  <  2:2567592‑M2/100‑52 (MQ=255)
                                                                       gACCTACGGCAGCGGTTACCTGGCAGGCAggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagc                              >  1:163968‑M2/1‑29 (MQ=255)
                                                                                                ggcaggctgtccggtatacagcggttcgccagcccctgcgtgttgtacgcatgtctggtctcatgctgccagagcagtgcttccgtctgcgggtgatgca     >  2:1082855‑M2/1‑4 (MQ=255)
                                                                                                ggcaggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagtgcttccgtctgcgggtgatgca     >  2:2363210‑M2/1‑4 (MQ=255)
                                                                                                ggcaggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagtgcttccgtctgcgggtgatgca     <  2:576108‑M2/100‑97 (MQ=255)
                                                                                                ggcaggctgtccggtatacagcggctcgccagcccctgcgcgttgtacgcaggtctggtctcaagctgccagagcagtgcttccgtctgcgggtgat        <  2:1551077‑M2/97‑94 (MQ=255)
                                                                                                 gcaggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagtgcttccgtctgcgggtgatgcac    <  2:854711‑M2/100‑98 (MQ=255)
                                                                                                 gcaggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagtgcttccgtctgcgggtgatgcac    >  2:227209‑M2/1‑3 (MQ=255)
                                                                                                  caggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtttcatgctgccagagcagtgcttccgtctgcgggtgatgcacc   <  1:308792‑M2/100‑99 (MQ=255)
                                                                                                  caggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagtgcttccgtctgcgggtgatgcacc   <  1:1658290‑M2/100‑99 (MQ=255)
                                                                                                   aggctgtgcggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagtgcttccgtctgcgggtgatgcaccg  >  1:1870777‑M2/1‑1 (MQ=255)
                                                                                                   aggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagtgcttccgtctgcgggtgatgcaccg  >  2:1959087‑M2/1‑1 (MQ=255)
                                                                                                   aggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagtgcttccgtctgcgggtgatgcaccg  >  2:850790‑M2/1‑1 (MQ=255)
                                                                                                   aggctgtccggtatacagcggttcgccagcccctgcgcgttgtacgcatgtctggtctcatgctgccagagcagtgcttccgtctgcgggtgatgcaccg  >  1:1335594‑M2/1‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGACATGCGTACAACGCGCAGGGGCTGGCGAACCGCTGTATACCGGACAGCCTGCCCGCCGTGGAATGGCTGACCTACGGCAGCGGTTACCTGGCAGGCATGAAACTCGGCGACACACCGCTGGTGGAGTACACCCGCGACCGCCTGCACCGGGAAACGCTGCGCAGCTTCGGCCGTTATGAACTCACCACCGCTTATA  >  NC_000913/3621505‑3621703

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