File: getSeq.py

package info (click to toggle)
ray 2.3.1-9
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 6,008 kB
  • sloc: cpp: 49,973; sh: 339; makefile: 281; python: 168
file content (27 lines) | stat: -rwxr-xr-x 403 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
 
#!/usr/bin/python3
"""Get a sequence from a fasta file by specifying the sequence name

Usage:
    getSeq.py fasta_file sequence_name
"""

# Importing modules
import sys

if len(sys.argv)!=3:
	print(__doc__)
	sys.exit(1)

file=sys.argv[1]
name=sys.argv[2]

dump=False

for line in open(file):
	if line[0]=='>':
		if line.find(name)>=0:
			dump=True
		elif dump:
			break
	if dump:
		print(line.strip())