sokrypton / sokrypton/ColabFold
Invalid database read error in colabfold_search
Nobody has claimed this yet.
- Dominant language
- Jupyter Notebook
- Stars
- 2.9k
- Forks
- 747
- PR merge metrics
- No merged PRs in 30d
Description
Expected Behavior
Hello,
I am trying to run batch searches against ColabFoldDB on a SLURM cluster, following the MSA instructions in the README.
Current Behavior
colabfold_search fails at the expandaln step with the error:
Invalid database read for database data file=[db_folder]/uniref30_2103_db.idx, database index=[db_folder]/uniref30_2103_db.idx.index
getData: local id (4294967295) >= db size (22)
Full log file: colabfold_search_output.txt
Steps to Reproduce (for bugs)
bash setup_databases.sh [db_folder]
Note:mmseqs createindexwas run with--split-memory-limit 128Gas mmseqs doesn't detect the SLURM job's memory limit otherwise.colabfold_search --db-load-mode 0 --mmseqs mmseqs_5185d3c/bin/mmseqs batch_1/input_sequences.fa [db_folder] batch_1/result_s8
Input sequences: input_sequences.fa
It looks like colabfold_search uses --split-memory-limit 0 in the prefilter steps and possibly later steps – I don't think this caused the issue as the job only reached 53 GB usage before it errored, but it would be nice to be able to set this to prevent the job from being killed.
Context
I'm looking to perform a batch search and the cluster jobs have a 250GiB limit, so I'm using --db-load-mode 0, but let me know if that isn't the best option.
Your Environment
- Git commit: 2a47c6f1459fbbdb5242cbc62173f9b513813cfa
- mmseqs commit: 5185d3cbb7af8a3122e202d47ddaaa785dc73890
- Server: Intel Xeon CPU with AVX, 256GiB memory (jobs limited to 250GiB, and lower limits can mean faster submission times)
- Operating system and version: CentOS 7
@thomashopf
Contributor guide
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Start with setup_databases.sh and the colabfold_search invocation, then inspect the expandaln step and the generated uniref30_2103_db.idx and .index files. Reproduce with mmseqs 5185d3c and --db-load-mode 0 under the stated memory limit, using the linked log and input sequence file; done means the batch search completes without the invalid database read error.
Written by the indexing model from the issue text.
Assessment
- Domain
- bioinformatics, cli
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Mostly clear
- Newbie friendliness
- 35/100