import sys,os,time
import multiprocessing as mp
bamfile = sys.argv[1]
#desDir = '/public0/nong/tcga_luad/bedgraph_test'
desDir = sys.argv[2]
os.system('mkdir -p ' + desDir)
bams = open(bamfile,'r').read().strip().split('\n')
cmd_filter = 'python3 '+ os.path.join(os.path.dirname(sys.argv[0]), '2.1_filter_chr.py')
def genomceCoverage(bam, outdir):
bamname = os.path.basename(bam)
graph = os.path.join(outdir,bamname ) + '.bedgraph'
cmd = 'genomeCoverageBed -bg -ibam '+bam + ' -g /public0/nong/source/hg38.chrom.chr.sizes -split | ' + cmd_filter + ' >' + graph
print(cmd)
os.system(cmd)
pool = mp.Pool(8) #多线程, 以及线程数!
for bam in bams:
t = pool.apply_async(genomceCoverage, (bam, desDir))
pool.close()
pool.join()
2018-03-10 python multiprocessing 多线程
最后编辑于 :
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 笔记如下 什么是多线程下载?图二 多线程下载的分析.jpg 首先要获得要下载文件的总大小并创建规定大小的空文件 每...
- python中的多线程就是在一个进程中存在着多个线程,在线程中,所有的线程都是共享资源的,线程之间的数据通信很简单...
- 进程:程序的一次执行(程序载入内存,系统分配资源运行)。每个进程有自己的内存空间,数据栈等,进程之间可以进行通讯,...
- IO 密集型应用 IO 密集型应用CPU等待IO时间远大于CPU 自身运行时间,太浪费;常见的 IO 密集型业务...
- 仅供交流学习 原文: https://rumenz.com/rumenbiji/python-requests-m...