New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 730041 =0 (0.000)12 (0.140) 7/168 NT 100% coding (459/4194 nt) rhsC Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 3619682 = NA (NA)coding (491/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

CCTGTATTTTGAGCACCTGTTTCCCGGTGAGGACGGTTACAGCCGCAGCGAGTCACTGTGGCTGGTGCGCGGCGGCGTGGCGAAACTGGATGAAGGTCACCGGCTGGCCGCACTCTGGCAGGCGCTGCCGGAAGAACTCCGCTTAAGTCCGCATCGTTATCTGGCGACAAACAGTCCGCAGGGGCCGTGGTGGCTGCTC  >  NC_000913/729942‑730140
                                                                                                   |                                                                                                   
catccgcttccggcacccgctcacaccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttc                                                                                                     <  1:1133598‑M1/1‑1 (MQ=255)
catccgcttccggcacccgctcacaccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttc                                                                                                     <  1:74569‑M1/1‑1 (MQ=255)
 atccgcttccggcacccgctcacaccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttcc                                                                                                    <  2:1008162‑M1/2‑1 (MQ=255)
 atccgcttccggcacccgctcacaccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttc                                                                                                     <  1:2170121‑M1/1‑1 (MQ=255)
  tccgcttccggcacccgctcacaccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttcc                                                                                                    >  1:399297‑M1/98‑99 (MQ=255)
  tccgcttccggcacccgctcacaccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCg                                                                                                   >  2:1927399‑M1/98‑100 (MQ=255)
  tccgcttccggcacccgctcacaccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCg                                                                                                   <  2:1936185‑M1/3‑1 (MQ=255)
   ccgcttccggcacccgctcacacccaccgagcagccaccacggcccctgcggactgtttgccgccagataacgatgcggacttcagcggagttcttCCgg                                                                                                  <  2:314752‑M1/4‑1 (MQ=255)
   ccgcttccggcacccgctcacaccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCgg                                                                                                  <  1:1040869‑M1/4‑1 (MQ=255)
   ccgcttccggcacccgctcacaccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCgg                                                                                                  <  2:134825‑M1/4‑1 (MQ=255)
    cgcttccggcacccgctcacaccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGc                                                                                                 <  2:274939‑M1/5‑1 (MQ=255)
    cgcttccggcacccgctcacaccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGc                                                                                                 <  2:274959‑M1/5‑1 (MQ=255)
    cgcttccggcacccgctcacaccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGc                                                                                                 <  2:1102374‑M1/5‑1 (MQ=255)
                         ccaaccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGC‑‑‑‑CGCACTCTGGCAGGCGCTGcc                                                                        >  1:890523‑M1/75‑100 (MQ=255)
                            cttccgatctcaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGCTGGCCGCACTCTGGCAGGCGCTGc                                                                         >  1:214986‑M1/72‑100 (MQ=255)
                            accgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGC‑‑‑‑CGCACTCTGGCAGGCGCTGCCGga                                                                     >  1:2361255‑M1/72‑100 (MQ=255)
                             ccgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGCTGGCCGCACTCTGGCAGGCGCTGcc                                                                        >  1:909769‑M1/71‑100 (MQ=255)
                              cgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGCTGGCCGCACTCTGGCAGGCGCTGCCg                                                                       >  1:991520‑M1/70‑100 (MQ=255)
                              cgagcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGCTGGCCGCACTCTGGCAGGCGCTGCCg                                                                       >  2:1729527‑M1/70‑100 (MQ=255)
                                 gcagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGCTGGCCGCACTCTGGCAGGCGCTGCCGgaa                                                                    >  1:1067591‑M1/67‑100 (MQ=255)
                                  cagccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGCTGGCCGCACTCTGGCAGGCGCTGCCGgaag                                                                   >  1:1362994‑M1/66‑100 (MQ=255)
                                   agccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGCTGGCTGCACTCTGGCAGGCGCTGCCGgaaga                                                                  >  1:1562015‑M1/65‑100 (MQ=255)
                                   agccaccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGCTGGCCGCACTCTGGCAGGCGCTGCCGgaaga                                                                  >  2:24597‑M1/65‑100 (MQ=255)
                                      caccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGCTGGCCGCACTCTGgcaggcgctgcagatcggaag                                                               <  2:214984‑M1/39‑11 (MQ=255)
                                      caccacggcccctgcggactgtttgtcgccagataacgatgcggacttaagcggagttcttCCGGCTGGCCGCACTCTGGCAGGCGCTGCCGGAAGAACt                                                               >  1:114835‑M1/62‑100 (MQ=255)
                                                                                                   CCGGCTGGCCGCACTCTGGCAGGCGCTGCCGGAAGAACTCCGCTTAAGTCCGCATCGTTATCTGGCGACAAACAGTCCGCAGGGGCCGTGGTGGCTGCTc  >  1:245053‑M1/1‑100 (MQ=255)
                                                                                                   CCGGCTGGCCGCACTCTGGCAGGCGCTGCCGGAAGAACTCCGCTTAAGTCCGCATCGTTATCTGGCGACAAACAGTCCGCAGGGGCCGTGGTGGCTGCTc  >  2:320280‑M1/1‑100 (MQ=255)
                                                                                                   CCGGCTGGCCGCACTCTGGCAGGCGCTGCCGGAAGAACTCCGCTTAAGTCCGCATCGTTATCTGGCGACAAACAGTCCGCAGGGGCCGTGGTGGCTGCTc  <  2:970054‑M1/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCTGTATTTTGAGCACCTGTTTCCCGGTGAGGACGGTTACAGCCGCAGCGAGTCACTGTGGCTGGTGCGCGGCGGCGTGGCGAAACTGGATGAAGGTCACCGGCTGGCCGCACTCTGGCAGGCGCTGCCGGAAGAACTCCGCTTAAGTCCGCATCGTTATCTGGCGACAAACAGTCCGCAGGGGCCGTGGTGGCTGCTC  >  NC_000913/729942‑730140

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.