Hi,
When I run annotate_gene_catalog command after gene prediction, the step completes without error but the output file is empty with only the headers present. This is the command I used:
MOCAT.pl -sf sample_file -annotate_gene_catalog -blasttype blastp
I have checked the gene prediction output and it looks good with many complete genes present so I am confused as to why I do not get any annotation results.
I checked the temp folder and the file sample_file-eggNOG.blastp has many entries but the subsequent files such as sample_file-eggNOG.filtered are empty. What can happen in the filtering step which removes all of the blasp results?
I have attached the log file here from which it appears all the steps ran correctly.
I have run the command with different samples and get the same result. Could there be some other issue with the way I am running the tool?
Would appreciate any help in fixing this problem.
Thanks,
The issue has been resolved. It was due to perl DBI module missing which was required by filterBlastReport.pl script though it is strange that this error did not show up in the logs. fixed it by installing the module.