Description Usage Arguments Value Author(s) References See Also Examples
Computes the nucleotide or amino acid frequencies at all polymorphic sites in the alignment.
1 | SummaryMuts(seqs, w = NULL, off = 0)
|
seqs |
DNAStringSet or AAStringSet with the haplotype sequences. |
w |
An optional numeric vector with the haplotype counts. When |
off |
Offset of first position in the alignment. |
Data frame with the polymorphic positions and nucleotide or amino acid frequencies.
Mercedes Guerrero-Murillo and Josep Gregori
Gregori J, Esteban JI, Cubero M, Garcia-Cehic D, Perales C, Casillas R, Alvarez-Tejado M, Rodríguez-Frías F, Guardia J, Domingo E, Quer J. Ultra-deep pyrosequencing (UDPS) data treatment to study amplicon HCV minor variants. PLoS One. 2013 Dec 31;8(12):e83361. doi: 10.1371/journal.pone.0083361. eCollection 2013. PubMed PMID: 24391758; PubMed Central PMCID: PMC3877031.
Ramírez C, Gregori J, Buti M, Tabernero D, Camós S, Casillas R, Quer J, Esteban R, Homs M, Rodriguez-Frías F. A comparative study of ultra-deep pyrosequencing and cloning to quantitatively analyze the viral quasispecies using hepatitis B virus infection as a model. Antiviral Res. 2013 May;98(2):273-83. doi: 10.1016/j.antiviral.2013.03.007. Epub 2013 Mar 20. PubMed PMID: 23523552.
1 2 3 4 5 6 | # Load haplotype alignment with abundances.
filepath<-system.file("extdata","ToyData_10_50_1000.fna", package="QSutils")
lst <- ReadAmplSeqs(filepath,type="DNA")
# Distribution of nucleotides at polymorphic sites.
SummaryMuts(lst$hseqs,lst$nr,off=0)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.