dodir_perseus.py

dodir_perseus.py

import os
import sys
import glob

def system(cmd):
    print >>sys.stderr, cmd
    os.system(cmd)

# dodir_pyro.py directory

directory = sys.argv[1]

owd = os.getcwd()
os.chdir(directory)

for ln in sys.stdin:
    cols = ln.split("\t")

    name = cols[0]
    barcode = cols[1]

    print name, barcode

    cmd = "Perseus -sin %s_s60_c01_T350_s25_c08_cd.fa > %s_s60_c01_T350_s25_c08_cd.per" % (name, name)
    print cmd

    #system(cmd)

os.system("../Class.sh %s" % (directory,))