fa_percent_GC: Get GC percent for multi-fasta file

In sethiyap/fastaR: Fasta Sequence Manipulation

Description

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

Usage

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

## Not run: 
 ref_fasta <- system.file("exdata", "Sc_nucl_R64-2-1.fasta", package = "fastaR")
 fastaR::fa_percent_GC(fasta_file=ref_fasta)

## End(Not run)

sethiyap/fastaR documentation built on June 16, 2020, 1:41 a.m.