新しいfile.txt
. ここでの私の目標は、複数のペアを開始して終了するかどうかを確認し、ペアごとに個別のファイルを作成することです。そうでない場合は、そのままにしておきます。
def searchPFAM(fname):
with open(fname,'rb') as f:
root = etree.parse(f)
for lcn in root.xpath("/protein/match[@dbname='PFAM']/lcn"):#find dbname =PFAM
try:
start = int(lcn.get("start"))#if it is PFAM then look for start value
end = int(lcn.get("end"))#if it is PFAM then also look for end value
yield start, end
except (TypeError , ValueError) as e:
pass
with open('newfile.txt','w') as fileinput:
for start, end in searchPFAM(fname):
print start, end
if start <= end:
text=''.join(makeList[(start-1):(end-1)])
fileinput.write(text)