BLAST to specify orthologs among genomes,
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?


blast ortholog evalue • 1.4k views
Entering edit mode

You should try it on protein level instead..


Login before adding your answer.

Traffic: 1865 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6