New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3765531 =NA (NA)9 (0.120) 3/168 NT 100% coding (3349/4134 nt) rhsA Rhs protein with putative toxin 55 domain; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor
?NC_000913 3765547 = 0 (0.000)coding (3365/4134 nt) rhsA Rhs protein with putative toxin 55 domain; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor

TGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTG  >  NC_000913/3765448‑3765625
                                                                                                   |                                                                              
gtcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcgg                                                                                >  2:999316‑M2/100‑100 (MQ=255)
gtcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcgg                                                                                >  1:1475267‑M2/100‑100 (MQ=255)
gtcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcgg                                                                                >  2:1010128‑M2/100‑100 (MQ=255)
gtcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcgg                                                                                >  1:253121‑M2/100‑100 (MQ=255)
gtcgcagtggtacaggtagatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcgg                                                                                >  2:1005269‑M2/100‑100 (MQ=255)
 tctcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggc                                                                               <  2:2228506‑M2/2‑1 (MQ=255)
 tcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggc                                                                               >  1:1113484‑M2/99‑100 (MQ=255)
 tcgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggc                                                                               <  2:1276868‑M2/2‑1 (MQ=255)
 ccgatctggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggc                                                                               <  1:1357151‑M2/2‑1 (MQ=255)
  cgcagtggtacaggtggatttttcgcgcctgcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcg                                                                              <  2:1420373‑M2/3‑1 (MQ=255)
  cgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcg                                                                              <  1:1118512‑M2/3‑1 (MQ=255)
  cgcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcg                                                                              <  1:37992‑M2/3‑1 (MQ=255)
   gcagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccctctggttttgcatctgctccccggtgaggccgcacgatgccagccagcggctg                                                                             >  2:1418176‑M2/97‑98 (MQ=255)
    cagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcggc                                                                            >  2:2009123‑M2/96‑100 (MQ=255)
    cagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcggc                                                                            <  2:1025334‑M2/5‑1 (MQ=255)
    cagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcggc                                                                            <  1:1787138‑M2/5‑1 (MQ=255)
    cagtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcaga                                                                            >  1:2228521‑M2/96‑97 (MQ=255)
     agtggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcggc                                                                            <  1:1337161‑M2/5‑1 (MQ=255)
      gtggtataggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagcgagcggcggc                                                                            <  1:84922‑M2/5‑1 (MQ=255)
        ggtacaggtggatttttcgcgccggcgtgtacaccgggtccatctggttttgcatctgctccacggtcaggccgcacgatgccagccagcggcagatcag                                                                        >  2:1357145‑M2/92‑93 (MQ=255)
                         cgcgccggcgtgtacaccgggtccatctggttttgcatctgctccaccgtcaggccgcacgatgccagccaacgGCGGCCTGACGGTGGAGCAGATGCaa                                                       >  1:231138‑M2/75‑100 (MQ=255)
                                                                      accgtcaggccgcacgatgccagccagcgGCGGCCTGACGGTGGAGCAGATGCAAAACCAGATGGACCCGATGTACACGCCGGCGCGAAAAATCCACCTg          >  1:1174478‑M2/30‑100 (MQ=255)
                                                                       cggtcaggccgcacgatgccagccaaagTCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCAGGTGTACACGCCGGCTCGAAAAATCCACCTGt         >  2:875926‑M2/29‑100 (MQ=255)
                                                                            aggccgcacgatgccagccagcgGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCAc    >  2:540318‑M2/24‑100 (MQ=255)
                                                                             ggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt   >  2:337509‑M2/23‑100 (MQ=255)
                                                                             ggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACACCGGCGCGAAAAATCCACCTGTACCACt   >  1:1283684‑M2/23‑100 (MQ=255)
                                                                             ggccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGAAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACt   >  1:1612287‑M2/23‑100 (MQ=255)
                                                                              gccgcacgatgccagccagcgGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTTCCACTg  >  2:1714946‑M2/22‑100 (MQ=255)
                                                                              gccgcacgatgccagccagcgGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTg  >  1:213389‑M2/22‑100 (MQ=255)
                                                                                                   |                                                                              
TGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGAGCAGATGCAAAACCAGATGGACCCGGTGTACACGCCGGCGCGAAAAATCCACCTGTACCACTG  >  NC_000913/3765448‑3765625

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