New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3620614NA (NA)7 (0.080) 6/166 NT 100% coding (1423/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 = 3620635 0 (0.000)coding (1444/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

ATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCC  >  NC_000913/3620531‑3620710
                                                                                   |                                                                                                 
aTGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCC                                                                                  >  1:2170080‑M2/1‑100 (MQ=255)
aTGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCC                                                                                  >  1:2435722‑M2/1‑100 (MQ=255)
 tgtgGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCT                                                                                 <  1:212673‑M2/100‑1 (MQ=255)
  gtggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTG                                                                                >  1:2320733‑M2/1‑100 (MQ=255)
   tggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGA                                                                               <  1:32768‑M2/100‑1 (MQ=255)
    ggtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGA                                                                               <  1:1611080‑M2/99‑1 (MQ=255)
     gtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACG                                                                             >  2:1333966‑M2/1‑100 (MQ=255)
     gtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGCCAGCCACCGGGCCTGACG                                                                             >  2:1150211‑M2/1‑100 (MQ=255)
     gtgACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAAACAAAACAACCAGTTAAAGTCAGCCACCGGGCCTGACG                                                                             <  1:869082‑M2/100‑1 (MQ=255)
          gggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGCCCGGTGGCTGACGttaac                                                                        >  1:940403‑M2/1‑95 (MQ=255)
           ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGCCCGGTGGCTGACGttaact                                                                       >  1:1576513‑M2/1‑94 (MQ=255)
           ggCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGCCCGGTGGCTGACGttaact                                                                       >  1:210629‑M2/1‑94 (MQ=255)
             ccTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGCCCGGTGGCTGACGttaactgg                                                                     >  2:2424627‑M2/1‑92 (MQ=255)
                                      gATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGGCCCGGTGGCTGACGttaactggttgtggtggttatagtaaaacgccg                                            >  2:1964786‑M2/1‑67 (MQ=255)
                                                                             cAGTTAACGTCAGCCACCGGGCCTGACGttaactggctgtggtggttatagtaaaacgccgatgccctgccatccggcgtggtgatgcgcgtgatgag       <  1:984331‑M2/98‑71 (MQ=255)
                                                                                  aaCGTCAGCCACCGGGCCTGACGttaaccggttgtggtggctatagtaaaacgccgatgccctgccatccggcgtggttatacgcgtgatcaggcccg  <  1:1635247‑M2/98‑76 (MQ=255)
                                                                                   |                                                                                                 
ATGTGGTGACGGGCCTCATCACGCGCATAACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCC  >  NC_000913/3620531‑3620710

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