Mercurial > repos > mheinzl > hd
comparison hd.py @ 3:82eaf30dd089 draft
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]), |