View source: R/blast_protein_to_protein.R
blast_protein_to_protein | R Documentation |
Run protein to protein BLAST of reference sequences against a blast-able database or fasta file.
blast_protein_to_protein(
query,
subject,
output.path = NULL,
is.subject.db = FALSE,
task = "blastp",
db.import = FALSE,
postgres.user = NULL,
evalue = 0.001,
out.format = "csv",
cores = 1,
max.target.seqs = 10000L,
db.soft.mask = FALSE,
db.hard.mask = FALSE,
blast.path = NULL
)
query |
path to input file in fasta format. |
subject |
path to subject file in fasta format or blast-able database. |
output.path |
path to folder at which BLAST output table shall be stored.
Default is |
is.subject.db |
logical specifying whether or not the |
task |
protein search task option. Options are:
|
db.import |
shall the BLAST output be stored in a PostgresSQL database and shall a connection be established to this database? Default is |
postgres.user |
when |
evalue |
Expectation value (E) threshold for saving hits (default: |
out.format |
a character string specifying the format of the file in which the BLAST results shall be stored. Available options are:
|
cores |
number of cores for parallel BLAST searches. |
max.target.seqs |
maximum number of aligned sequences that shall be retained. Please be aware that |
db.soft.mask |
shall low complexity regions be soft masked? Default is |
db.hard.mask |
shall low complexity regions be hard masked? Default is |
blast.path |
path to BLAST executables. |
Hajk-Georg Drost
blast_nucleotide_to_nucleotide
, blast_nucleotide_to_protein
,
blast_protein_to_nucleotide
, blast_best_hit
## Not run:
blast_example <- blast_protein_to_protein(
query = system.file('seqs/qry_aa.fa', package = 'metablastr'),
subject = system.file('seqs/sbj_aa.fa', package = 'metablastr'),
output.path = tempdir(),
db.import = FALSE)
# look at BLAST results
blast_example
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.