Online BLAST search using R
Is there a way to incorporate online BLAST search using an R script which takes FASTA sequences as input and returns the BLAST hits?
• 4,209 views
•
link
updated
by
SequenceServer •
•
written
by
accibio •
•
written
by
accibio •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Fastest way to perform BLAST search using a multi-FASTA file against a remote database
written by accibio •
-
Obtaining FASTA format from blastp output
written by mbrav005 •
-
Extract sequences which do not have blast hits
written by karthic
-
"Similarity search" for epitopes
written by gil.hornung
-
creating graph of distribution of blast hits on the query sequence
written by hersh •
-
Retrieve only BLAST hit sequences from a de novo transcriptome
written by Cristian Vergara •
-
tool: changes in binding affinities caused by SNP
written by tonja.r
-
Online blast and standalone blast return different hits?
written by ishengomae
-
A Tool For Flagging Blast Searches?
written by ishengomae
-
Python Program Accepting Fasta And Blast, Outputting In Xml Then Returning Ids Of All Sequences
written by te502 •
Why do you wish to automate a web crawl on BLAST when blast has a command line utility with online search enabled?
I'm trying to develop an R package for extracting annotations for proteins based on results from BLAST searches. So I need to incorporate BLAST inside the R package.
Why not take BLAST results as input instead? Why not accept a FASTA file instead of a sequence? Use R for data processing, not reinventing the wheel on existing processes. For example, while an R package to parse BAM files is useful, an R package to call variants from BAM is not - not everything needs to be re-implemented in every language.
Maybe this helps https://www.rdocumentation.org/packages/base/versions/3.6.2/topics/system2 or this https://www.r-bloggers.com/2020/01/working-with-windows-cmd-system-commands-in-r/
metablastris probably what you're looking for. Or you could try using theBLASTAPI.I agree with Ram though. You're probably better off telling your users to do their own sequence searches, and just bring your package the output in a standardized format (which you can prescribe). Especially if you're just going to end up telling them to run
BLASTlocally (in which case you should be telling them to useMMseqs2at this point because it's faster and lightweight).Check this R package (https://github.com/gschofl/blastr and https://rdrr.io/github/gschofl/blastr/man/). It looks like BLAST with R. But I failed to find any tutorials about this package.
It's possible with SequenceServer. We've made examples for blasting from within python. And for remote blasting in the unix command-line.
The same concept (i.e. server-side API) should be accessible from R too. But we have yet to make example code.