Get GC percent for multi-fasta file

Source: R/fa_utils.R

for given multi-fasta file get GC percent for each sequence

fa_percent_GC(fasta_file)

Arguments

fasta_file

Either a path or a connection to multi-fasta file. The input sequence file should have extention .fa or .fasta In the sequence header: only string before first space and/or first colon (:) will be considered for futher processes. **Important consideration when header have big names.

Value

A tibble of GC percent and a barplot of GC percent distribution if the fasta file contains sequences less than 20

Examples

if (FALSE) {
 ref_fasta <- system.file("exdata", "Sc_nucl_R64-2-1.fasta", package = "fastaR")
 fastaR::fa_percent_GC(fasta_file=ref_fasta)
}