-
Notifications
You must be signed in to change notification settings - Fork 0
/
select_seq_by_ID.py
33 lines (28 loc) · 1013 Bytes
/
select_seq_by_ID.py
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
28
29
30
31
32
33
# Python script to recover sequences from a fasta file from a list of IDs.
# Run the script as following: python select_seq_by_ID.py <fasta file> <ID file> <output file>
import sys
header = []
sequence = []
count = -1
file_open = open(sys.argv[1],"r") #fasta file to be analyzed
for line in file_open:
line = line.rstrip()
if ">" in line:
temp=line.split(" ")
header.append(temp[0])
count+=1
sequence.insert(count,'')
else:
sequence[count]+=line
con=count+1
print ("The number of sequences in this file is:") #for checking
print (con)
result=""
file3 = open(sys.argv[2], "r") #file contianing the IDs of sequences, one per line
file4 = open(sys.argv[3], "w") # output file
for line in file3:
line=line.rstrip()
for i in range (count):
if (line in header[i]) & (len(line) == len(header[i])-1):
result=str(header[i])+"\n"+ str(sequence[i]) + "\n"
file4.write(str(result))