View source: R/blast_nucleotide_to_protein.R
blast_nucleotide_to_protein | R Documentation |
Run nucleotide to protein BLAST of reference sequences against a blast-able database or fasta file. Internally BLAST translates the nucleotide sequence into a protein sequence and then searches for hits.
blast_nucleotide_to_protein(
query,
subject,
strand = "both",
output.path = NULL,
is.subject.db = FALSE,
task = "blastx",
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. |
strand |
Query DNA strand(s) to search against database/subject. Options are:
|
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 |
nucleotide 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_protein_to_protein
,
blast_protein_to_nucleotide
, blast_best_hit
## Not run:
blast_test <- blast_nucleotide_to_protein(
query = system.file('seqs/qry_nn.fa', package = 'metablastr'),
subject = system.file('seqs/sbj_aa.fa', package = 'metablastr'),
output.path = tempdir(),
db.import = FALSE)
# look at results
blast_test
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.