OrthoFinder
OrthoFinder copied to clipboard
Running with pre-blast results gets stuck
Hi,
I put 3587 protein seqs into orthofinder with -op
option.
Then I finished 12866569 times BLASTP via parallel
by my own.
It may continue to finish the rest of the pipeline using -b
option pointing to the BLAST directory but the program seems to be stuck in some progress.
JOBID NAME USER CPUS MIN_MEMORY TIME NODELIST(REASON)
15417698 orthofinder_358 jzye 32 0 2-03:37:43 cu315
orthofinder -b Results_Aug28/WorkingDirectory -t ${threads} -a ${threads} -M msa -T fasttree
The log file only presents:
OrthoFinder version 2.5.5 Copyright (C) 2014 David Emms
2024-09-03 11:33:45 : Starting OrthoFinder 2.5.5
32 thread(s) for highly parallel tasks (BLAST searches etc.)
32 thread(s) for OrthoFinder algorithm
Checking required programs are installed
----------------------------------------
Test can run "mcl -h" - ok
Test can run "mafft" - ok
Test can run "fasttree" - ok
strace
this job id the it turns out to be the program is reading all BLAST files repeatly:
stat("/lustre/home/jzye/E.faecalis/Phylogeny/all_ncbi_accessions/3587_Samples/Results_Aug28/WorkingDirectory/Blast454_563.txt", 0x7fff5250c9a0) = -1 ENOENT (No such file or directory)
stat("/lustre/home/jzye/E.faecalis/Phylogeny/all_ncbi_accessions/3587_Samples/Results_Aug28/WorkingDirectory/Blast454_563.txt.gz", {st_mode=S_IFREG|0644, st_size=155324, ...}) = 0
stat("/lustre/home/jzye/E.faecalis/Phylogeny/all_ncbi_accessions/3587_Samples/Results_Aug28/WorkingDirectory/Blast454_564.txt", 0x7fff5250c9a0) = -1 ENOENT (No such file or directory)
stat("/lustre/home/jzye/E.faecalis/Phylogeny/all_ncbi_accessions/3587_Samples/Results_Aug28/WorkingDirectory/Blast454_564.txt.gz", {st_mode=S_IFREG|0644, st_size=145655, ...}) = 0
Can anyone give me some advice? Thanks!