New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3620123 =NA (NA)4 (0.030) 3/266 NT 100% coding (932/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
?NC_000913 3620147 = 0 (0.000)coding (956/4236 nt) rhsB Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

ACCGTCCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAG  >  NC_000913/3619998‑3620296
                                                                                                                                                     |                                                                                                                                                     
acggtaacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccag                                                                                                                                                       <  2:933267‑M2/1‑1 (MQ=255)
 cggtaacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagc                                                                                                                                                      >  1:403340‑M2/149‑150 (MQ=255)
 cggtaacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagc                                                                                                                                                      >  2:1043057‑M2/149‑150 (MQ=255)
 cggtaacggatttccggccggcccgtgggacgttgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagc                                                                                                                                                      <  2:1334890‑M2/2‑1 (MQ=255)
 cggtaacggatttccggccggcacgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagc                                                                                                                                                      >  2:543257‑M2/149‑150 (MQ=255)
  ggtaacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcg                                                                                                                                                     >  1:222294‑M2/148‑150 (MQ=255)
   gtaacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcgg                                                                                                                                                    <  1:1145457‑M2/4‑1 (MQ=255)
   gtaacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcgg                                                                                                                                                    >  2:670811‑M2/147‑150 (MQ=255)
   gtaacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcgg                                                                                                                                                    <  1:585085‑M2/4‑1 (MQ=255)
   gtaacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcgg                                                                                                                                                    >  2:415496‑M2/147‑150 (MQ=255)
    taacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggc                                                                                                                                                   <  1:1262086‑M2/5‑1 (MQ=255)
    taacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggc                                                                                                                                                   <  2:544255‑M2/5‑1 (MQ=255)
    taacggatttccggccggcccgtgtgacggtgcgccaccatccggcccaggtatttataatcgtaagtaaagctgcgcacctgtttgccactacggtcataccccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggc                                                                                                                                                   >  1:1971113‑M2/146‑150 (MQ=255)
     aacggatttccggccggcccgtgtgacggtgcgccaccatccggccccggtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcg                                                                                                                                                  >  1:549426‑M2/145‑150 (MQ=255)
                                                                                                                  cagttcgccgcgtggcgtccagccatagcgcaccaGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAAt                                     >  2:1184633‑M2/36‑150 (MQ=255)
                                                                                                                    gttcgccgcgtggcgtccagccatagcgcaccaGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAACTACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATcc                                   >  2:971045‑M2/34‑150 (MQ=255)
                                                                                                                       cgccgcgtggcgtccagccatagcgcaccaGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATccgtt                                >  1:1632794‑M2/31‑150 (MQ=255)
                                                                                                                       cgccgcgtggcgtccagccatagcgcaccaGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATccgtt                                >  2:933672‑M2/31‑150 (MQ=255)
                                                                                                                                                     GCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACATGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAg  >  1:609475‑M2/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGTCCGCTTTTCCTGATACCCTGCCGGGTTACACCGAATATGGCCGGGACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGGCCGGAAATCCGTTACCGTTACGACAGCGACGGGCGGGTGACAG  >  NC_000913/3619998‑3620296

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