Finding 1 pfam in metagenome

[ElmavWieren]

Hi,

I am currently using this tool to find one gene in a metagenome (knowing that is not exactly what this tool is meant to do). The following command seems to work for this particular pfam: cblaster search -m hmm -qp PF08545 -db place/where/dmnd/fasta/sqlite3/are -pfam place/of/pfam/db -p test.html -s test.json -o test.csv -mh 1 -u 1. But when I try to use a different pfam (PF10520) or try to use the -qf option, the resulting files do not contain any hits even though the output is this:

[10:18:06] INFO - Starting hmmer search [10:18:06] INFO - Pfam database found [10:18:06] INFO - Fetching profiles from Pfam-A file [10:18:06] INFO - Profiles found: ['PF10520.12'] [10:18:23] INFO - Saved profiles to /work/elma/pfam/cblaster_20220308101823.hmm [10:18:23] INFO - Performing hmmsearch [10:18:37] INFO - Found 19 hits meeting score thresholds for hmm search [10:18:37] INFO - Fetching genomic context of hits [10:18:37] INFO - Querying local SQLite3 database: /work/elma/cblaster_dbs/all_bins_db.sqlite3 [10:18:38] INFO - Searching for clustered hits across 19 organisms [10:18:38] INFO - Writing current search session to /work/elma/cblaster_output/carF_all_bins_db.json [10:18:38] INFO - Writing summary to /work/elma/cblaster_output/carF_all_bins_db.csv [10:18:38] INFO - Saving cblaster plot HTML to: /work/elma/cblaster_output/carF_all_bins_db.html [10:18:38] INFO - Done.

Do you have any suggestions about what I could do, to make this work? Or is this not possible with cblaster?

Thank you!