replicateReads: Replicate sequences from DNAStringSet object using counts...

In malnirav/hiReadsProcessor: Functions to process LM-PCR reads from 454/Illumina data

Description

Given a DNAStringSet object, the function replicates reads using counts=X marker at the end of definition line.

Usage

replicateReads(dnaSet, counts = NULL)

Arguments

dnaSet	DNAStringSet object to replicate.
counts	an integer or a numeric vector of length length(dnaSet) indicating how many times to repeat each sequence. Default is NULL, in which it uses counts=X notation from the definition line to replicate reads.

Value

DNAStringSet object.

See Also

dereplicateReads, removeReadsWithNs, findBarcodes, splitByBarcode

Examples

 
dnaSet <- c("CCTGAATCCTGGCAATGTCATCATC", "ATCCTGGCAATGTCATCATCAATGG", 
"ATCAGTTGTCAACGGCTAATACGCG", "ATCAATGGCGATTGCCGCGTCTGCA", 
"CCGCGTCTGCAATGTGAGGGCCTAA", "GAAGGATGCCAGTTGAAGTTCACAC", 
"CCTGAATCCTGGCAATGTCATCATC", "ATCCTGGCAATGTCATCATCAATGG", 
"ATCAGTTGTCAACGGCTAATACGCG", "ATCAATGGCGATTGCCGCGTCTGCA", 
"CCGCGTCTGCAATGTGAGGGCCTAA", "GAAGGATGCCAGTTGAAGTTCACAC") 
dnaSet <- dereplicateReads(dnaSet)
replicateReads(dnaSet)

malnirav/hiReadsProcessor documentation built on July 29, 2021, 6:33 a.m.