Could not open .../round-2/families.stk for reading!
Closed this issue · 1 comments
I was trying to use this new version but there is an error at round-2. The last several lines of the STDOUT are shown below:
Comparison Time: 00:00:13 (hh:mm:ss) Elapsed Time, 561 HSPs Collected
- RECON: Running imagespread..
RECON Elapsed: 00:00:00 (hh:mm:ss) Elapsed Time - RECON: Running initial definition of elements ( eledef )..
RECON Elapsed: 00:00:00 (hh:mm:ss) Elapsed Time - RECON: Running re-definition of elements ( eleredef )..
RECON Elapsed: 00:00:00 (hh:mm:ss) Elapsed Time - RECON: Running re-definition of edges ( edgeredef )..
RECON Elapsed: 00:00:00 (hh:mm:ss) Elapsed Time - RECON: Running family definition ( famdef )..
RECON Elapsed: 00:00:00 (hh:mm:ss) Elapsed Time - Obtaining element sequences
Number of families returned by RECON: 262
Processing families with greater than 15 elements
Could not open /mnt/tiger/osj/TE/Selaginella/S_lep_V4.2/RepeatModeler2/RM_46181.ThuMay111045252017/round-2/families.stk for reading!
In the round-2 directory, there was no families.stk file, but in the parent directory I could see the file families.stk was created and data were dumped in it.
I think the bug lies around line 534 in RepeatModeler where the directory of families.stk is assigned:
my $combFamiliesFile = "$tmpDir/families.stk";
But it was a bit complicated for me to debug it.
The dependent program versions:
RepeatMasker/4.0.6
RMBLAST/2.6.0
RECON/1.08
trf404.linux64
RepeatScout/1.0.5
I also installed RepeatModeler/1.0.8 and it worked well with the above dependent programs.
Thanks for the report. You are correct. There is a bug in the code that is exercised when no repeats are detected in a round. I have patched the MASTER branch and will make a release shortly.