def exon_graph(depth, refseq, prefix, outdir=os.getcwd(), trans=None, genes=None, threads=0): beds = HGVS(refseq, trans=trans, genes=genes) messages = beds.get_exons() threads = cpu_count() if not int(threads) else int(threads) pool = multiprocessing.Pool(processes=threads) for transcript, trans_messages in messages.iteritems(): pool.apply_async(graph, (depth, trans_messages, outdir, prefix)) pool.close() pool.join()
