process_fasta: Read and Process FASTA Files.

Description Usage Arguments Value Author(s) Examples

Description

Function to read and process FASTA files into single-line sequences, R compatible objects.

Usage

1
process_fasta(path_to_fasta, min_seq_len)

Arguments

path_to_fasta

Path to a gzipped FASTA file.

min_seq_len

(int) Minimum sequence length. Sequences smaller than this value will not be included in the sequences object. Default: (5000).

Value

Returns vector containing per-Chromosome/Scaffold sequences.

Example: "ACGTACGT...", "ACGTACGT...", "ACGTACGT...", ...

Author(s)

Angel G. Rivera-Colon

Examples

1
2
3
4
5
6
7
# Path to FASTA file
myfasta <- system.file("extdata",
                       "test_geno.fa.gz",
                       package = "RADseqTools")

# Create a sequences object
mySeqs <- process_fasta(myfasta, 10000)

angelgr2/radseq_tools documentation built on May 15, 2019, 3:59 a.m.