@files( (("../data/znf.data", "profile.fasta" ) , ))
def createProfile( infile, outfile ):
    '''convert mali to profile
    '''
    outfile = open(outfile, "w")
    for line in open(infile):
        if line.startswith("#"): continue
        data = re.split("\s+", line[:-1],1)
        print data
        pid, sequence = data
        outfile.write(">%s\n%s\n" % (pid, sequence ))
    outfile.close()
dataset/ETHPy150Open CGATOxford/cgat/obsolete/pipeline_prdm9.py/createProfile/original