Commit 58d65692 authored by Christina Floristean's avatar Christina Floristean
Browse files

Exclude uniprot hits from regular msa feature processing

parent 96f46b08
...@@ -725,7 +725,8 @@ class DataPipeline: ...@@ -725,7 +725,8 @@ class DataPipeline:
) )
# The "hmm_output" exception is a crude way to exclude # The "hmm_output" exception is a crude way to exclude
# multimer template hits. # multimer template hits.
elif(ext == ".sto" and not "hmm_output" == filename): # Multimer "uniprot_hits" processed separately.
elif(ext == ".sto" and filename not in ["uniprot_hits", "hmm_output"]):
msa = parsers.parse_stockholm(read_msa(start, size)) msa = parsers.parse_stockholm(read_msa(start, size))
else: else:
continue continue
...@@ -741,7 +742,7 @@ class DataPipeline: ...@@ -741,7 +742,7 @@ class DataPipeline:
if(ext == ".a3m"): if(ext == ".a3m"):
with open(path, "r") as fp: with open(path, "r") as fp:
msa = parsers.parse_a3m(fp.read()) msa = parsers.parse_a3m(fp.read())
elif(ext == ".sto" and not "hmm_output" == filename): elif(ext == ".sto" and filename not in ["uniprot_hits", "hmm_output"]):
with open(path, "r") as fp: with open(path, "r") as fp:
msa = parsers.parse_stockholm( msa = parsers.parse_stockholm(
fp.read() fp.read()
......
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