New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3620119 =NA (NA)11 (0.090) 7/166 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

GACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACA  >  NC_000913/3620047‑3620245
                                                                                                   |                                                                                                   
gtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggc                                                                                                     <  2:600167‑M2/1‑1 (MQ=255)
gtatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggc                                                                                                     <  1:2322327‑M2/1‑1 (MQ=255)
gatcttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggc                                                                                                     <  1:667724‑M2/1‑1 (MQ=255)
 tatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcg                                                                                                    >  2:1644688‑M2/99‑100 (MQ=255)
 tatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcg                                                                                                    >  2:827948‑M2/99‑100 (MQ=255)
 tatttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcg                                                                                                    <  1:892418‑M2/2‑1 (MQ=255)
   tttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcg                                                                                                  >  2:467531‑M2/97‑100 (MQ=255)
   tttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcg                                                                                                  >  2:2764115‑M2/97‑100 (MQ=255)
   tttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcg                                                                                                  <  1:2229855‑M2/4‑1 (MQ=255)
   tttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcg                                                                                                  <  1:2120569‑M2/4‑1 (MQ=255)
    ttatcatcgtaagtaaagctgcgcacctgtttgccactacgttcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcgg                                                                                                 >  2:2213312‑M2/96‑100 (MQ=255)
    ttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcgg                                                                                                 >  2:550709‑M2/96‑100 (MQ=255)
    ttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcgg                                                                                                 <  2:17093‑M2/5‑1 (MQ=255)
    ttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcgg                                                                                                 >  2:3010828‑M2/96‑100 (MQ=255)
    ttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcgg                                                                                                 <  2:2626969‑M2/5‑1 (MQ=255)
    ttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcgg                                                                                                 >  2:695254‑M2/96‑100 (MQ=255)
    ttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcgg                                                                                                 >  1:1162443‑M2/96‑100 (MQ=255)
    ttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcgg                                                                                                 <  2:844926‑M2/5‑1 (MQ=255)
    ttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcgg                                                                                                 <  1:467529‑M2/5‑1 (MQ=255)
    ttatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcgg                                                                                                 >  1:2421783‑M2/96‑100 (MQ=255)
     tatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggc                                                                                                >  2:131565‑M2/95‑100 (MQ=255)
     tatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggc                                                                                                >  2:1397569‑M2/95‑100 (MQ=255)
     tatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggc                                                                                                <  2:416677‑M2/6‑1 (MQ=255)
     tatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggc                                                                                                >  2:449396‑M2/95‑100 (MQ=255)
     tatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggc                                                                                                <  1:1206858‑M2/6‑1 (MQ=255)
     tatcatcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcagatc                                                                                                >  2:667727‑M2/95‑95 (MQ=255)
          tcgtaagtaaagctgcgcacctgtttgccactacggtcatacaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgaactggcgg                                                                                           <  1:1543108‑M2/11‑10 (MQ=255)
                gtaaagctgcgcacctgtttgccactacggtcatagaccaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggcgcggcgggtacatt                                                                                       >  1:1984566‑M2/84‑90 (MQ=255)
                                                    accaccgccagttcgccgcgtggcgtccagccatagcgcaccagcggCGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCtt                                                 >  2:1843187‑M2/48‑100 (MQ=255)
                                                         cgccagttcgccgcgtggcgtccagccatagcgcaccagcggCGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTttactt                                            >  1:1805450‑M2/43‑100 (MQ=255)
                                                            ccgccagttcgccgcgtggcgtccagccatagcgcaccaGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACga                                        >  2:1343965‑M2/40‑100 (MQ=255)
                                                              gttcgccgcgtggcgtccagccatagcgcaccagcggCGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACgat                                       >  2:2543146‑M2/38‑100 (MQ=255)
                                                              gttcgccgcgtggcgtccagccatagcgcaccagcggCGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGGGCGCAGCTTTACTTACga                                        >  1:3267807‑M2/38‑99 (MQ=255)
                                                               ttcgccgcgtggcgtccagccatagcgcaccagcggCGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACgatg                                      >  2:3443383‑M2/37‑100 (MQ=255)
                                                                 cgccgcgtggcgtccagccatagcgcaccagcggCGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACgatgat                                    >  2:3373698‑M2/35‑100 (MQ=255)
                                                                  gccgcgtggcgtccagccatagcgcaccagcggCGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATa                                   >  1:2605319‑M2/34‑100 (MQ=255)
                                                                  gccgcgtggcgtccagccatagcgcaccagcggCGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATa                                   >  2:1635211‑M2/34‑100 (MQ=255)
                                                                      gccgcgtggcgtccagccatagcgcaccaGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATAc                              >  2:2868419‑M2/30‑100 (MQ=255)
                                                                      gccgcgtggcgtccagccatagcgcaccaGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATAc                              >  2:2264055‑M2/30‑100 (MQ=255)
                                                                                                   CGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTcaca  >  2:2599598‑M2/1‑100 (MQ=255)
                                                                                                   CGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTcaca  <  1:900977‑M2/100‑1 (MQ=255)
                                                                                                   CGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGAGGATAAGTACCGGGGTCGGATGGTGGCGCACCGTcaca  >  1:2512956‑M2/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GACAACGGCATCCGTCTGTCTGCCGTGTGGCTGACGCACGACCCGGAATACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCACGCGGCGAACTGGCGGTGGTGTATGACCGTAGTGGCAAACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACA  >  NC_000913/3620047‑3620245

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