#! /usr/bin/env python # from __future__ import division # / operator performs float rather than int division import os, sys, string, cp import SeqIO def main(): """ %(filename)s A template for building scripts using the cp, SeqIO, and Seq modules. {{out name of outfile. Prints to stdout by default.}} {{h Prints documentation.}} {{v verbosity 1}} {{version print version info and exit}} %(version)s """ debug = 0 # copy from here... docstringdict = {'filename':os.path.split(sys.argv[0])[-1], 'version':'$Id: cp_template.py,v 1.1 2004/10/10 02:53:10 nghoffma Exp $'} optlist, formattedDocString, formattedSummary = cp.optStringParser(main.__doc__ % docstringdict, optWidth=15, lineWidth=60, valOffset = 4, putVals=1) dict = cp.commandparser(options=optlist, usage='Options:\n' + formattedSummary, debug = debug, exitWithUsage=1) if dict.status('version'): sys.exit( 'Version: %(version)s\n' % docstringdict ) #...to here if dict.status('h'): print formattedDocString sys.exit( 0 ) v = dict.value('v') #to read one fasta file as a file object fastafile = dict.value('infile', 'in','file') fasta = SeqIO.readFasta(fastafile, output='list') #read from a list of fasta files fasta_list = dict.value('infile_list', 'in','file') fasta = SeqIO.readFastaList(fasta_list, output='list') #write to stdout by default writeToFile = 0 if dict.status('out'): outfile = dict.value('out','out','file') writeToFile = 1 else: outfile = sys.stdout v = 0 if writeToFile: outfile.close() if __name__ == '__main__': main()