BLAST to specify orthologs among genomes,
0
0
Entering edit mode
5.3 years ago

Greeting all,

I want to align all the orthologous genes of some genomes, so i could obtain the genome tree.

I had generated the fasta files, each file contained all the genes of one genome.

I used BLAST for each pair of genomes, and generated the lists.

From each genome I selected one gene, s.t. each gene was the best hit of the other genes.

But I obtained only less than 20 orthologous genes with just 3 genomes,

while more than 1000 genes were obtained with 2 genomes.

When more genomes was included, almost no orthologs could be found.....

Any one have efficient methods? Is the method I used too strict?

btw, this is the option I used for blastn :

-task blastn -> I tried blastn-short, but the result was somehow worst

-evalue 5 -> not sure what value to be

-max_target_seqs 5 -> when larger, or not to set this option, I got a little more orthologs ...... why?

Thanks.

blast ortholog evalue • 1.4k views
ADD COMMENT
0
Entering edit mode

You should try it on protein level instead..

ADD REPLY

Login before adding your answer.

Traffic: 1865 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