Orthofinder:How to obtain single copy homologous sequences from result files
orthofinder -fg RESULTS_DIR -M msa -oa
The command described above does not yield a single copy of the homologous gene.
I also have a question here, that is, we use protein sequences for alignment and clustering analysis, and if I finally want to get the corresponding nucleotide sequence, what should I do?
Thanks in advance!
• 8,324 views
•
link
updated
by
sunnykevin97
•
written
by
Zhenpeng Yu •
•
written
by
Zhenpeng Yu •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
pal2nal -nogap -nomismatch options not working
written by pkfsantos •
-
Sequence identity between two protein sequences
written by Gene_MMP8
-
How to select multiple domains from protein sequences to align them and build a phylogenetic tree
written by Michela •
-
MSA of single copy gene sequences using OrthoVenn2 output
written by shanewidanagama •
-
multiple sequence alignment using HMM and simulated annealing
written by Rojhin •
-
How can a MSA be useful for sub-cellular location identification?
written by fabbrl •
-
Multiple Alignment of several homologous gene
written by kaihami •
-
[Noob Doubt] Mapping and Annotation.
written by Kato •
-
Multiple Alignment: Protein Or Nucleotide Sequence?
written by tomas.iesmantas •
-
Novel Strategy For Improving Sequence Alignment
written by Pals
We could do with a little more info here. Eg. what species are you working with? what data do you have at hand? What is the ultimate goal?
what do you mean with
thx
I'm interested in eight marine mammals, and now I've found more than 8000 single-copy homologous genes through orthofinder, and I want to extract the results through programs or scripts, but I'm not very good at programming myself. So ask if there is a corresponding software or script.
Thanks!
OK, so you want to obtain all the single-copy genes that are present in each of the 8 mammals , or simply all single copy genes (which can be only present in a few species) ?
Without going into detail: all this info can be parsed from the CSV output files you should get. Do you get those CSV files?
Can you also explain why you use that specific orthofinder cmdline? It's not wrong but it's not the 'default' one either
The command line I described above is from the orthofinder description document.
This issue has also been discussed, with the following links https://github.com/davidemms/OrthoFinder/issues/154
As you already used protein sequences of 8 marine species and found 8000 single copy genes (proteins) in result directory with orthogroups assigned by orthofinder, you can pick gene ids of each species from
SingleCopyOrthogroups.txtfile,like this:
then you can retrieve nucleotide sequences of
OrthologsIDS.txtfrom species`s nucleotide files.You're amazing. The
grepcommand is so cool that the problem is solved instantly. I sincerely thank you.The grep command is not working I tried for my data. Help is needed @ Mehmet.
Hello, I want to know how do you sort out the sequence IDs contained in each homology group after you get the OrthologsIDS.txt file? thank you very much
Can you define your question with more details? What do you want to do? In OrthologsIDS.txt file there are ortholog families with one gene in each ortholog family from each species.
I have obtained single copy orthogroups and corresponding amino acid sequence files containing orthogroups. How can I obtain the corresponding nucleotide sequence based on the id of these amino acid sequence files (such as XPxxxx of species A)? Thanks in advance.
You need something like this?
extract sequences based on ids file
Thank you very much. It solved my problem.