read_base_content {qckitfastq} | R Documentation |
Compute nucleotide content per position for a single base pair. Wrapper function around seqTools.
read_base_content(fseq, content)
fseq |
a seqTools::fastqq object |
content |
nucleotide. Options are "A","T","G","C","N"(either capital or lower case) |
Nucleotide sequence content per position.
Wenyue Xing, wenyue_xing@brown.edu, August Guang august_guang@brown.edu
infile <- system.file("extdata", "10^5_reads_test.fq.gz", package = "qckitfastq") fseq <- seqTools::fastqq(infile,k=6) read_base_content(fseq,"A")