makeList[start]のインデックスをmakeList[end]に出力しようとしています。開始値と終了値はすべて正しく、合計で約 11 ペアあります。ただし、可能なすべての出力が出力されないため、正しい出力が得られません。なぜそれが起こっているのですか?
def searchPFAM(fname):
with open(fname,'rb') as f:
root = etree.parse(f)
for lcn in root.xpath("/protein/match[@dbname='PFAM']/lcn"):
try:
start = int(lcn.get("start"))
end = int(lcn.get("end"))
yield start, end
except (TypeError , ValueError) as e:
pass
uniprotID ='Q14591'
def createRef():
fileInput = open(uniprotID +'.txt','r')
readAA = fileInput.read()
store_value = (readAA.partition('\n'))
aaSequence = store_value[2].replace('\n', '')
makeList = list(aaSequence)
return makeList
makeList = createRef()
for start, end in searchPFAM(fname):
print start, end
if start <= end:
with open('newfile.txt','w') as fileinput:
while start <= end:
result = makeList[start]
fileinput.write(result)
start += 1