Commit f37b7117 authored by Lukas Jelonek's avatar Lukas Jelonek
Browse files

Include evaule in hmmer result. Use name as fallback if no accession is available

parent 0f9b5e49
......@@ -48,21 +48,28 @@ with open(filename) as f:
print('HMMER converter: Query ID "' + query_id + '" not found among initial queries!')
results = documents[query_id]['computations'][0]["results"]
dbxref = None
if split[1] == "-":
dbxref = args.dbxref + ":" + split[0]
else:
dbxref = args.dbxref + ":" + split[0]
results.append({
'target': {
'name': split[0],
'dbxref': args.dbxref + ':' + split[1],
'length': split[2],
'score': float(split[13]),
'start': int(split[15]),
'end': int(split[16]),
'description': split[22].rstrip()
'name': split[0], # target name
'dbxref': dbxref,
'length': split[2], # tlen
'score': float(split[13]), # this domain score
'evalue': float(split[12]), # this domain i-Evalue
'start': int(split[15]), # hmm coord from
'end': int(split[16]), # hmm coord to
'acc': float(split[21]), # acc
'description': split[22].rstrip() # description of target
},
'query': {
'start': int(split[17]),
'end': int(split[18]),
'envelop_start': int(split[19]),
'envelop_end': int(split[20])
'start': int(split[17]), # ali coord from
'end': int(split[18]), # ali coord to
'envelop_start': int(split[19]), # env coord from
'envelop_end': int(split[20]) # env coord to
}
})
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment