OrthoFinder icon indicating copy to clipboard operation
OrthoFinder copied to clipboard

Running with pre-blast results gets stuck

Open yyliang12 opened this issue 5 months ago • 0 comments

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
image This job has been running for 2 days long and nothing output to the log file. The command I used was:
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!

yyliang12 avatar Sep 05 '24 07:09 yyliang12