擅长:python、mysql、java
<p>看来你正试图找到它们的交叉点。Pybedtools将接受流作为输入。将数据读入bed格式的字符串。你知道吗</p>
<p>“chr,开始,停止”</p>
<p>我从一个python字典开始,循环浏览它。你知道吗</p>
<pre><code>bed_string += "{0} {1} {2} {3} {0}|{1}|{2}|{3}\n".format(chrom, coord_start, coord_stop, aberration)
# Now create your bedtools.
breakpoint_bedtool = pybedtools.BedTool(bed_string, from_string=True)
target_bedtool = pybedtools.BedTool(self.args.Target_Bed_File, from_string=False)
# Find target intersects for printing.
breakpoint_target_intersect = breakpoint_bedtool.intersect(target_bedtool, wb=True, stream=True)
</code></pre>