read_base_content {qckitfastq}R Documentation

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

Description

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

Usage

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

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

[Package qckitfastq version 1.4.0 Index]