fastBlockSearch profile read error
I downloaded the current PRINTS profiles from their ftp and converted them to profiles using the given prints2prfl.pl script.
I also concatenated the .prfl files from the Arthropod BUSCO (http://busco.ezlab.org/) dataset.
I then concatenated the PRINTS and BUSCO profiles together to scan my genome with fastBlockSearch:
fastBlockSearch --cutoff=1.1 genome.fa aug.prfl > fbs.out
I get a couple of warnings about profile being not significant enough which I assume to be normal. Then I get an error:
Warning: Block unknown_B is not significant enough, removed from profile.
Warning: Block unknown_F is not significant enough, removed from profile.
Warning: Block unknown_G is not significant enough, removed from profile.
Warning: Block unknown_M is not significant enough, removed from profile.
terminate called after throwing an instance of 'PP::ProfileReadError'
Aborted (core dumped)
What am I doing wrong? Should I not concatenate .prfl files using cat? Maybe the .prfl in BUSCO aren't in expected format?
Not really related to my problem, but the current PRINTS profile .dat file actually contained an error in one of the lines that causes prints2prfl.pl to throw an error. The error references a line in the .data file that needs to be modified for the script to run properly. Basically you just have to change the row header "fd" to "id".