Last chance! 50% off unlimited learning
Sale ends in
Take vectors and/or matrices arguments and combine them row-wise without
recycling them (as is the case with rbind).
seqbind(..., blank = "-")vectors, matrices, and/or alignment ‘fasta’ objects to combine.
a character to add to short arguments, to achieve the same length as the longer argument.
Returns a list of class "fasta" with the following components:
an alignment character matrix with a row per sequence and a column per equivalent aminoacid/nucleotide.
sequence names as identifers.
the matched call.
Grant, B.J. et al. (2006) Bioinformatics 22, 2695--2696.
# NOT RUN {
## Read two pdbs
a.pdb <- read.pdb("1bg2")
b.pdb <- read.pdb("1goj")
seqs <- seqbind(aa321(a.pdb$atom[a.pdb$calpha,"resid"]),
aa321(b.pdb$atom[b.pdb$calpha,"resid"]))
# seqaln(seqs)
# }
Run the code above in your browser using DataLab