Can Orthofinder produce FastA?
0
0
Entering edit mode
5.5 years ago
bsp017 ▴ 50

Is it possible to get a FastA file of orthologs using Orthofinder? I ran orthofinder (v2.2.7) like this:

orthofinder.py -f sane_files/ -t 40 -M msa

I got around 9000 orthogroups from my genomes. I would like to Blast these against some databases e.g. CAZy database. But I need a FastA file.

Thanks,

James

Orthofinder CAZy sequence alignment • 1.3k views
ADD COMMENT
1
Entering edit mode

Would the following option be relevant in your case?

os               Stop after writing sequence files for orthogroups
                   (requires '-M msa')
ADD REPLY
0
Entering edit mode

Adding that command gives the proteins that make up the orthogroup. I was wanting a single consensus sequence for each orthogroup. Not sure if there is an easy way to do that? For example, I got 9000 orthogroups, so I wanted a FastA file with 9000 proteins.

ADD REPLY

Login before adding your answer.

Traffic: 1826 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6