I am trying to subset a fasta file like described in another post (How To Extract Multiple Fasta Sequences At A Time From A File Containing Sequences Ids Using R). I simply want to be able to extract a subset of sequences from a fasta file.
I load the package and the data as described in the other thread:
library("seqinr")
fastafile<- read.fasta(file = "proteins.fasta",
seqtype = "AA",as.string = TRUE, set.attributes = FALSE)
as well as the list of IDs I was to use for subsetting (the column in subsetlist is called id).
subsetlist<-read.table("~/scripts/subsetfasta/test.txt", header=TRUE)
when I attempt to use the solution from the previous thread:
fastafile[names(fastafile) %in% subsetlist$id]
I get the following:
named list()
What am I doing wrong or missing?
best regards
Henrik
•
written
by
henrikkjeldal •
3 answers
Looks like none of your fasta record names are in the subset list. You can check by just calling
names(fastafile) %in% subsetlist$id
Which will be a vector of FALSEs if that's the case.
You probably want to check exactly how the subset list and fasta file IDs are formatted, and use (g)sub, paste or related functions to get them to match.
R solution:
fastafile[c(which(names(fastafile) %in% subsetlist$id))]
•
written
by
Tanvir Ahamed
Not a R Solution:
$ ./faSomeRecords main.fasta id.txt output.fa
option -exclude will output sequences not present in main.fasta
•
written
by
Tanvir Ahamed
Log in to answer this question.
-
BLAST options error: File homo_ref.faa does not exist
written by logbio •
-
errors in multiple sequencing alignment in r by use msa package!!
written by kasemmohmmed287 •
-
Fasta Files - How to search for peptide sequences?
written by janainamace •
-
how to create a fasta file with subset of sequences using seqinr
written by Camilla •
-
read.fasta() function couldn't found
written by sherlockyz •
-
how to subset sequences based on IDs, but order shall remain preserved?
written by majeedaasim •
-
Amino Acid Substitution Loop Generates Incorrect Output: No Mutations in the Output, but the same s…
written by jguevar •
-
fasta file loading error
written by tcf.hcdg
-
How To Extract Multiple Fasta Sequences At A Time From A File Containing Sequences Ids Using R
written by catagui •
-
Seqinr Cannot Recognize ">" Character In File
written by Partha •
Dear all,
Does someone know how to do partial matching based on a column in a dataframe (here column id in dataframe subsetlist)? in case id thus only contains part of the names of the fasta file? (since %in% needs a complete match)
Thanks! ellen
Please open a new question and add a reproducible example.