powered by
Crop genetic sequences with a set of primers
refdb_clean_seq_crop_primers( x, primer_forward, primer_reverse, max_error_in = 0.1, max_error_out = 0.1, include_primers = TRUE )
A reference database.
a reference database with a defined sequence field.
primer forward.
primer reverse.
maximum error for a match (frequency based on primer length).
a logical indicating whether the detected primers are included in the cropped sequences.
lib <- read.csv(system.file("extdata", "baetidae_bold.csv", package = "refdb")) lib <- refdb_set_fields_BOLD(lib) refdb_clean_seq_crop_primers(lib, "AGT", "TTTA")
Run the code above in your browser using DataLab