Mercurial > repos > mheinzl > hd

comparison hd.py @ 4:ceb3a9b4bef5 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/monikaheinzl/galaxyProject/tree/master/tools/hd commit 61a78acfd6dc00ec3bb7f5550f4f81875c164907
author mheinzl
date Tue, 15 May 2018 10:47:06 -0400
parents 82eaf30dd089
children 445442a6aeb0
comparison
equal deleted inserted replaced
3:82eaf30dd089 4:ceb3a9b4bef5
803 # for h, tag in zip(ham, result1): 803 # for h, tag in zip(ham, result1):
804 # output_file1.write("{}\t{}\n".format(tag, h)) 804 # output_file1.write("{}\t{}\n".format(tag, h))
805 805
806 # HD analysis for chimeric reads 806 # HD analysis for chimeric reads
807 proc_pool_b = Pool(nproc) 807 proc_pool_b = Pool(nproc)
808 print(chunk_sample) 808 print(chunks_sample)
809 print(result2) 809 print(result2)
810 print(data_array) 810 print(data_array)
811 diff_list_a = proc_pool_b.map(partial(hamming_difference, array2=result2, mate_b=False), chunks_sample) 811 diff_list_a = proc_pool_b.map(partial(hamming_difference, array2=result2, mate_b=False), chunks_sample)
812 diff_list_b = proc_pool_b.map(partial(hamming_difference, array2=result2, mate_b=True), chunks_sample) 812 diff_list_b = proc_pool_b.map(partial(hamming_difference, array2=result2, mate_b=True), chunks_sample)
813 proc_pool_b.close() 813 proc_pool_b.close()