Description Usage Arguments Value Author(s) See Also Examples
Mask out low-quality bases in a DNA sequence prior to downstream analyses.
1 | qualityMask(seq, max.err)
|
seq |
A QualityScaledDNAStringSet containing one or more sequences with quality strings. |
max.err |
Numeric scalar specifying the error probability threshold, above which bases are to be masked. |
A DNAStringSet containing masked DNA sequences.
All bases with error probabilities above max.err
will be replaced with N
s.
Aaron Lun
1 2 3 4 5 6 7 8 | X <- QualityScaledDNAStringSet(
DNAStringSet(paste0("ACAGCTGCTGCAGCTGTG")),
PhredQuality(rbeta(18, 1, 10))
)
qualityMask(X, max.err=0.5)
qualityMask(X, max.err=0.1)
qualityMask(X, max.err=0.01)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.