Created
November 7, 2013 16:14
-
-
Save samuell/7357232 to your computer and use it in GitHub Desktop.
Base complementer example in Papy, trying to add multple workers in to a single piper
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| __author__ = 'samuel' | |
| from papy.core import Worker, Plumber, Piper | |
| from papy.core import start_logger | |
| from numap import NuMap | |
| infile_name = "Homo_sapiens.GRCh37.67.dna_rm.chromosome.Y.fa.short" | |
| def base_complement( inbox ): | |
| def basecompl_nucl(nucl): | |
| if nucl == 'A': | |
| return 'T' | |
| elif nucl == 'T': | |
| return 'A' | |
| elif nucl == 'C': | |
| return 'G' | |
| elif nucl == 'G': | |
| return 'C' | |
| else: | |
| return '' | |
| sequence = inbox[0] | |
| new_sequence = "" | |
| for nucl in sequence: | |
| new_sequence += basecompl_nucl(nucl) | |
| return new_sequence | |
| def print_text( inbox ): | |
| text = inbox[0] | |
| print text | |
| local_computer = NuMap() # worker_num = 1 needed in order not to change ordering of lines | |
| base_complementer = Worker(base_complement) | |
| printer = Worker(print_text) | |
| base_complementer_node = Piper((base_complementer, base_complementer), parallel=local_computer) | |
| printer_node = Piper(printer, track=True) | |
| network = Plumber() | |
| network.add_pipe((base_complementer_node, printer_node)) | |
| end = network.get_outputs()[0] | |
| network.start([(line for line in open(infile_name))]) | |
| network.run() | |
| network.wait() | |
| print "Finished!" |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment