view STACKS_clone_filter.py @ 5:1a5ebb16aaea draft

planemo upload commit 4d00a0b7a4b859b642eceb3c7cd4adb1453903a0-dirty

#!/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__()