Mercurial > repos > tiagoantao > clone_filter
view STACKS_clone_filter.py @ 7:2059c028c07a draft
planemo upload commit 4d00a0b7a4b859b642eceb3c7cd4adb1453903a0-dirty
author | tiagoantao |
---|---|
date | Thu, 04 Feb 2016 08:00:32 -0500 |
parents | 1a5ebb16aaea |
children | a8b87bc0c3ad |
line wrap: on
line source
#!/usr/bin/python # -*- coding: utf-8 -*- import argparse import glob import os import subprocess import shutil def __main__(): parser = argparse.ArgumentParser() parser.add_argument('-f', type=str) parser.add_argument('-fwd', type=str) parser.add_argument('-rev', type=str) #parser.add_argument('-k', action='store_true') options = parser.parse_args() # create the working dir os.mkdir('job_outputs') os.chdir('job_outputs') # create the populations command input line cmd_line=['clone_filter'] if options.f is not None: cmd_line.extend(['-f', options.f]) else: cmd_line.extend(['-1', options.fwd]) cmd_line.extend(['-2', options.rev]) print "[CMD]:"+' '.join(cmd_line) subprocess.call(cmd_line) # postprocesses clean_files = glob.glob('*.fq*') if options.f is not None: shutil.copy(clean_files[0], 'clean.fq') else: shutil.copy(clean_files[0], 'fwd.fq') shutil.copy(clean_files[1], 'rev.fq') if __name__ == '__main__': __main__()