Count sequences in a FASTA file

Counts the number of sequences in a FASTA file by counting header lines (lines starting with >). Accepts gzip files.

Usage

count_seq_db(file)

Arguments

file

(Character, required) Path to a FASTA file (plain or gzip).

Value

An integer, the number of sequences.

See also

count_pattern_db()

Author

Adrien Taudière

Examples

db <- system.file("extdata", "example_unite.fasta", package = "dbpq")
count_seq_db(db)
#> [1] 5