reciprocal fasta script at Sanger this fasta script is too old it's need PSU module which is not fully available. But the idea of reciprocal fasta is good -I could make two protein blasts with key "-m 8" and "-b 1" and I'll have two tables which i should only compare
Ortholog Detection Tool
Dear all nowadays I want to detect orthologs in wide group of bacteria, and I've choose blast method BRH. I've found quite good topics here about it http://biostar.stackexchange.com/questions/5104/how-to-prove-an-ortholog-does-not-exist-in-a-given-genome
http://biostar.stackexchange.com/questions/6161/blast-reblast-to-verify-homologies
http://biostar.stackexchange.com/questions/6161/blast-reblast-to-verify-homologies
But a lot of tools described there aren't exist or don't work correctly on my workstation (debian testing 16+8 gb RAM)
So could I ask you about software tool of detection orthologs or way to write a script in bash or in python Thank You for attention
2 answers
You might want to check the basic approaches of "reciprocal FASTA" where usually hits are identified as amino acid identity >40% over at least 80% of the total length of the protein....
For the mapping you can basically use whatever mapper you find useful. Usually FASTA, MUMmer, or BLAST are used.... You missed this biostar thread about it.
A simple google gave me the reciprocal fasta script at Sanger.
A starter at least.
Instead of predicting them you can refer to the already predicted groups of ortholog proteins in the NCBI Protein Clusters database: http://www.ncbi.nlm.nih.gov/proteinclusters . I use this database to find orthologs of bacterial genes and it works quite good.
Log in to answer this question.
More posts like this
-
Cog Paralogous And Orthologous Classes
written by Zhe Liu •
-
Graphic Of Genomic Positions Of Protein Domains
written by Tonig
-
Whole Genome Phylogeny
written by Gorysko
-
How To Monitor Shared Computing Resources?
written by Jim Vallandingham
-
Bam To Wig With Fixed Step?
written by 2184687-1231-83-
-
How Do You Analyze Metagenomics Data For Species Content And Biological Function?
written by Larry_Parnell
-
Most Sensitive Method To Search 100-1000Bp Sequences To Genome At Human-Mouse Genetic Distances
written by 2184687-1231-83-
-
Blast/Reblast To Verify Homologies
written by Mawe •
-
Genomic Sequences For Drosophila In Refseq
written by Dejian
-
Short Read Aligners/ Snp Calling Comparison
written by Pierre Lindenbaum
I don't think bash/python will help you in this case. It sounds like you need to work on your system administration skills. What tools have you tried to install on your system and what do you mean by "don't work correctly"? Error messages are helpful.
OrthoMCL and Inparanoid Inparanoid gives only empty files OrthoMCL has problems with MySQL