官方文档:https://pyranges.readthedocs.io/en/latest/autoapi/pyranges/index.html
持续更新中
pyrange结构可以存储bed区域及其注释内容。
需求:
例如需要对2个bed文件区域取交集
import pandas as pd
import pyranges as pr
###需注意,intersect时,列名必须是'Chromosome','Start','End'
bed1 = pd.read_csv(bed1_path,sep='\t',names=['Chromosome','Start','End'])
bed2 = pd.read_csv(bed2_path,sep='\t',names=['Chromosome','Start','End'])
pr1 = pr.PyRanges(bed1)
pr2 = pr.PyRanges(bed2)
df_filter = pr1.intersect(pr2).as_df()
pr1
pr2