From 024d8a1fe8da4d1096cba30b333f8d61ee4343ac Mon Sep 17 00:00:00 2001 From: Eric Coissac Date: Tue, 15 Dec 2009 21:38:16 +0000 Subject: [PATCH] add two new scripts : fastaHead and fastaTail --- src/fastaHead.py | 37 +++++++++++++++++++++++++++++++++++++ src/fastaTail.py | 41 +++++++++++++++++++++++++++++++++++++++++ 2 files changed, 78 insertions(+) create mode 100644 src/fastaHead.py create mode 100644 src/fastaTail.py diff --git a/src/fastaHead.py b/src/fastaHead.py new file mode 100644 index 0000000..c700d33 --- /dev/null +++ b/src/fastaHead.py @@ -0,0 +1,37 @@ +#!/usr/local/bin/python +''' +Created on 15 dŽc. 2009 + +@author: coissac +''' +from obitools.fasta import fastaIterator,formatFasta +from obitools.options import getOptionManager + +import sys + +def addHeadOptions(optionManager): + optionManager.add_option('-n','--sequence-count', + action="store", dest="count", + metavar="###", + type="int", + default=10, + help="Count of first sequences to print") + + +if __name__ == '__main__': + optionParser = getOptionManager([addHeadOptions], + entryIterator=fastaIterator + ) + + (options, entries) = optionParser() + i=0 + + for s in entries: + if i < options.count: + print formatFasta(s) + i+=1 + else: + sys.exit(0) + + + diff --git a/src/fastaTail.py b/src/fastaTail.py new file mode 100644 index 0000000..83e818e --- /dev/null +++ b/src/fastaTail.py @@ -0,0 +1,41 @@ +#!/usr/local/bin/python +''' +Created on 15 dŽc. 2009 + +@author: coissac +''' +from obitools.fasta import fastaIterator,formatFasta +from obitools.options import getOptionManager +import collections + +def addHeadOptions(optionManager): + optionManager.add_option('-n','--sequence-count', + action="store", dest="count", + metavar="###", + type="int", + default=10, + help="Count of first sequences to print") + + +if __name__ == '__main__': + optionParser = getOptionManager([addHeadOptions], + entryIterator=fastaIterator + ) + + (options, entries) = optionParser() + i=0 + + queue = collections.deque() + + + for s in entries: + queue.append(s) + if len(queue) > options.count: + queue.popleft() + + while queue: + print formatFasta(queue.popleft()) + + + +