Mercurial > repos > mheinzl > hd

comparison hd.py @ 3:82eaf30dd089 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 5e42600b728963d8559753e8c82845372d9dc3cc
author mheinzl
date Tue, 15 May 2018 10:43:32 -0400
parents 316fbf91dd12
children ceb3a9b4bef5
comparison
equal deleted inserted replaced
2:316fbf91dd12 3:82eaf30dd089
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)
809 print(result2)
810 print(data_array)
808 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)
809 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)
810 proc_pool_b.close() 813 proc_pool_b.close()
811 proc_pool_b.join() 814 proc_pool_b.join()
812 diff = numpy.concatenate((numpy.concatenate([item[0] for item in diff_list_a]), 815 diff = numpy.concatenate((numpy.concatenate([item[0] for item in diff_list_a]),