Entering edit mode
3.2 years ago
devhimd
▴
10
# Program to find the indexes of Cys in the given mutlifasta sequences
fasta = open('out.fa', 'r+')
for line in fasta.read().split('\n'):
if line.startswith(">"):
header = line
print(header)
else:
indexes = []
for i in range(0, len(line)-1):
if line[i] == 'C':
indexes.append(i+1)
print("Cys : ", indexes)
Is there a question?