read_base_content: Compute nucleotide content per position for a single base...

Description Usage Arguments Value Author(s) Examples

View source: R/read_content.R

Description

Compute nucleotide content per position for a single base pair. Wrapper function around seqTools.

Usage

1
read_base_content(fseq, content)

Arguments

fseq

a seqTools::fastqq object

content

nucleotide. Options are "A","T","G","C","N"(either capital or lower case)

Value

Nucleotide sequence content per position.

Author(s)

Wenyue Xing, wenyue_xing@brown.edu, August Guang august_guang@brown.edu

Examples

1
2
3
infile <- system.file("extdata", "10^5_reads_test.fq.gz", package = "qckitfastq")
fseq <- seqTools::fastqq(infile,k=6)
read_base_content(fseq,"A")

qckitfastq documentation built on Nov. 8, 2020, 5:24 p.m.