This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Blast plus command line filter !

Hello All,

I am doing blast with "NR" database which is taking very long time and I am interested in searching only specific species. Is there a way to submit in the blast option to restrict it to search for specific species. I am using blast plus. In older blast there is command blastall -u "organism name". Is there any similar option in blast plus other than creating separate database from "NR".

Thanks

unix blast command-line
• 4,815 views • link
updated by Ram Ram • written by empyrean999 empyrean999

1 answer

Check this older post about setting up masked database for BLAST (Vertebrate Subset Nr Database? Build My Own?). In short, you download GI list for the species you're interested, use blastdb_aliastool with the -gilist flag to create a restricted database, and use this database to blast your sequences.

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

In addition, blast+ programs (blastn, blastp, etc) accept -gilist, -seqidlist or -negative_gilist to restrict the search to subsets of your database.

• 0 views • link
updated by Ram Ram • written by h.mon h.mon

Thank you. I will look in to that

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

Log in to answer this question.

More posts like this