How can I extract sequence name from fast file in R

[xshang]

Here is my code:

setwd("/Users/xshang/Dropbox/project-ing/peptidome")
library("Biostrings")
human_protein <- readAAStringSet("HUMAN.fasta",use.names=TRUE)

When I use

names(human_protein[[1]])

, it returns "NULL". How can I extract sequence name from fast file in R?

Thanks!

[dpryan]

Try

names(human_protein)[1]

The names() function is an accessor for the names in the AAStringSet object and returns a character vector. Alternatively, you can also use

names(human_protein[1])

The [ function also returns an AAStringSet, rather than an AAString (returned by [[).

[xshang]

Try

names(human_protein)[1]

The names() function is an accessor for the names in the AAStringSet object and returns a character vector. Alternatively, you can also use

names(human_protein[1])

The [ function also returns an AAStringSet, rather than an AAString (returned by [[).

It works! Many thanks!