New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3620119 =NA (NA)5 (0.040) 3/266 NT 100% coding (928/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 3620146 = 0 (0.000)coding (955/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

AACCGTCCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCC  >  NC_000913/3619997‑3620263
                                                                                                                                                     |                                                                                                                     
taacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggc                                                                                                                       <  2:1729385‑M2/1‑1 (MQ=255)
taacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggc                                                                                                                       >  1:2048933‑M2/150‑150 (MQ=255)
taacggatatccggcgggcccgtgtgccggtgcgccaccatccggccccggtatttatcctcgtcagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagtacgccgcgtggcgtccagccatagcgcaccagcggc                                                                                                                       <  1:280254‑M2/1‑1 (MQ=255)
 aacggatttccggccggcccgtgtgacggtgcgccaccatccggcccgggtatttatcatcgtaagtaaagctgcgcgcctgtttgccacgacggtcatgcgccaccgccggttcgccgcgtggcgtccaggcatggcgcaccaggggcg                                                                                                                      >  2:878950‑M2/149‑150 (MQ=255)
 aacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcg                                                                                                                      <  1:1591696‑M2/2‑1 (MQ=255)
 aacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcg                                                                                                                      >  1:1894484‑M2/149‑150 (MQ=255)
  acggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgcccctacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgc                                                                                                                     >  1:1939525‑M2/148‑150 (MQ=255)
  acggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcataccccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgc                                                                                                                     >  1:207167‑M2/148‑150 (MQ=255)
  acggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacacccccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgc                                                                                                                     >  1:1603061‑M2/148‑150 (MQ=255)
  acggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgccccagcggcgc                                                                                                                     >  1:914784‑M2/148‑150 (MQ=255)
   cggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcg                                                                                                                    >  1:1629752‑M2/147‑150 (MQ=255)
   cggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcg                                                                                                                    <  2:1650747‑M2/4‑1 (MQ=255)
   cggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcg                                                                                                                    >  2:115510‑M2/147‑150 (MQ=255)
   cggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcaccagtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcg                                                                                                                    >  1:1608096‑M2/147‑150 (MQ=255)
    ggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatgatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcgg                                                                                                                   <  1:1666887‑M2/5‑1 (MQ=255)
    ggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcgg                                                                                                                   >  1:1003579‑M2/146‑150 (MQ=255)
     gatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggc                                                                                                                  <  1:1736906‑M2/6‑1 (MQ=255)
                                                                                                         cgaaTTTACCTGCCGCGCCGCGGGTGCGCTAGGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTATACTTACGATGAAAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCC              >  1:1615366‑M2/2‑150 (MQ=255)
                                                                                                         agaATTTCCCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTGACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCC              >  1:1100399‑M2/1‑150 (MQ=255)
                                                                                                         agaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCC              <  2:2151245‑M2/150‑1 (MQ=255)
                                                                                                            aTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC           >  1:1459707‑M2/1‑150 (MQ=255)
                                                                                                            aTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC           >  1:1104326‑M2/1‑150 (MQ=255)
                                                                                                            aTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGC           <  2:319526‑M2/150‑1 (MQ=255)
                                                                                                             tttACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCTGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC          <  1:252894‑M2/150‑1 (MQ=255)
                                                                                                             tttACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC          <  1:515319‑M2/150‑1 (MQ=255)
                                                                                                             tttACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC          <  1:5644‑M2/150‑1 (MQ=255)
                                                                                                             tttACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCC          >  2:1505639‑M2/1‑150 (MQ=255)
                                                                                                             tttACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGGGTATGCCCGTAGTGGCAAACAGGTGCGCAGCGTTACTTACGATGCTAAAGACCGGGGCCGGATGGTGGCGCACCGTCACCCGGGCCGGCC          <  1:837575‑M2/150‑1 (MQ=255)
                                                                                                                gttcgccgcgtggcgtccagccatagcgcaccagcggCGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGa       >  1:1497943‑M2/38‑150 (MQ=255)
                                                                                                                gttcgccgcgtggcgtccagccatagcgcaccagcggCGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGa       >  2:469937‑M2/38‑150 (MQ=255)
                                                                                                                   cgccgcgtggcgtccagccatagcgcaccagcggCGAGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAAt    >  1:1658960‑M2/35‑150 (MQ=255)
                                                                                                                    gttcgccgcgtggcgtccagccatagcgcaccaGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATcc  >  2:2004408‑M2/34‑150 (MQ=255)
                                                                                                                    gccgcgtggcgtccagccatagcgcaccagcggCGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAGTc   >  2:420978‑M2/34‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                     
AACCGTCCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCC  >  NC_000913/3619997‑3620263

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