Skip to content

Commit

Permalink
ENH Use hmmscan to match AMPSphere website
Browse files Browse the repository at this point in the history
This is now the exact same code as on the website
  • Loading branch information
luispedro committed Jun 28, 2024
1 parent 564ac7b commit e60b0ff
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion macrel/ampsphere.py
Original file line number Diff line number Diff line change
Expand Up @@ -102,7 +102,7 @@ def get_ampsphere_hmmer_match_local(args, seqs):
for (query_name, seq) in seqs:
f.write(f'>{query_name}\n{seq}\n')
_logged_subprocess_call(
['hmmsearch', '--tblout', output_file, hmm, query_file])
['hmmscan', '--domtblout', output_file, hmm, query_file])
return pd.read_csv(output_file, sep='\s+', comment='#', header=None,
names=['target', 'target_accession', 'query_name',
'query_accession', 'evalue', 'score', 'bias',
Expand Down

0 comments on commit e60b0ff

Please sign in to comment.