Blaster implements an efficient BLAST-like sequence comparison algorithm, written in C++11 and using native R datatypes. Blaster is light-weight, fast and dependency-free. The code base of Blaster is adapted from nsearch. An implementation of nsearch for Python is available at npysearch.
conda install -c conda-forge r-blaster
::install_github("tamminenlab/blaster") devtools
# Read a query file into DataFrame
<- read_fasta("inst/extdata/query.fasta")
query
# Read a database file into a DataFrame
<- read_fasta("inst/extdata/db.fasta")
db
# BLAST the query against the database
<-
blast_table blast(query, db)
# BLAST protein sequence file against itself using filenames as blast function arguments
<-
prot_blast_table blast(query = "inst/extdata/prot.fasta",
db = "inst/extdata/prot.fasta",
alphabet = "protein")
# Filter the sequences containing motif GAGACTT
<- read_fasta("query.fasta", "GAGACTT") query
Details available at https://anaconda.org/conda-forge/r-blaster/files.