fastBlockSearch profile read error

Discussions about predicting genes with AUGUSTUS. Not covered here: WebAUGUSTUS and BRAKER1

Moderator: bioinf

Post Reply
User avatar
katharina
Site Admin
Posts: 531
Joined: Wed Nov 18, 2015 6:14 pm
Location: Greifswald
Contact:

fastBlockSearch profile read error

Post by katharina » Thu Nov 19, 2015 3:18 pm

Originally posted in the old forum by Damian Kao on 08.07.2015 - 11:47
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'
what(): std::exception
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".
Katharina

Post Reply