This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Open source All vs All Blast script?

I am looking to an all vs all blast between two genomes, does anyone know of a script that does this?

Thank you

genome blast
• 4,400 views • link
updated by Ram Ram • written by lufiealf • lufiealf

Are you looking to do just the single-way comparison or reciprocal BLAST?

• 0 views • link
updated by Ram Ram • written by DG DG

Do you know of anything to All v All compare invertebrate genomes?

• 0 views • link
written by lufiealf • lufiealf

2 answers

You just have to have two fasta files of e.g. proteins of the two genomes, you run formatdb on one of them and then run blastall with one fasta as input file (-i) and the other as database (-d). Use -m8 to get the output in tabular format and you'll get a tab-delimited file with all against all proteins.

• 0 views • link
written by Asaf Asaf

Although if you want reciprocal best hits, which you often do, you have to do the other comparison and then compare the two results.

• 0 views • link
updated by Ram Ram • written by DG DG

LAST is able to compare two vertebrate genomes

• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

Do you know of anything to All v All compare invertebrate genomes?

• 0 views • link
written by lufiealf • lufiealf

It should work for invertebrates as well. It isn't restricted to vertebrate genomes. The only limiting factor is really the size of the genomes you are comparing to each other

• 0 views • link
written by DG DG

Log in to answer this question.

More posts like this