phyformat provides the ability to read and write sequence data in FASTA and PHYLIP formats, including the ability to interconvert between the two. It also includes a function to take a subset of a Newick-format phylogenetic tree.
1 2 3 4 5 6 | x <- ReadFasta("sequences.fa") #Creates a physeq object from an existing FASTA file
x$Sequences[[1]] <- "PQHLIRV" #Edit a sequence in the physeq object
WriteFasta(x, "sequences2.fa") #Write to a fasta file
WritePhylip(x, "sequences2.phy") #Write to a PHYLIP file
y <- physeq(c("Seq 1", "Seq 2"), c("RRARRPR", "PKKKRKV")) #Creates a physeq object
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.