Hi everyone, I am very new in environment R. I am trying to obtain a final fasta file selecting a subset of sequences from a downloaded fasta file, using seqinr. For this purpose I type the following script:
myfasta<- read.fasta(file = "MRP2.fas", seqtype = "AA",as.string = TRUE, set.attributes = FALSE)
subsetlist<-read.table("list.txt", header=TRUE)
myfasta[names(myfasta) %in% subsetlist$id]
write.fasta(sequences = myfasta, names = names(myfasta), nbchar = 80, file.out = "parsedMRP2.fasta")
parsedMRP2<- read.fasta("parsedMRP2.fasta", set.attributes = FALSE)"
The MRP2.fas contains 5039 sequences, while the list.txt is a list of 3000 sequences. When I run my script, I obtain final file (parsedMRP2.fasta) always with 5039 sequences instead of the 3000 sequences that I indicated in the list.txt.
Could you please help me? thank you very much in advance
1 answer
Because you are not saving data in new object. Before you write in the file, you need to save in another object. See the change below.
myfasta<- read.fasta(file = "MRP2.fas", seqtype = "AA",as.string = TRUE, set.attributes = FALSE)
subsetlist<-read.table("list.txt", header=TRUE)
my_fasta_sub <- myfasta[names(myfasta) %in% subsetlist$id]
write.fasta(sequences = my_fasta_sub, names = names(my_fasta_sub), nbchar = 80, file.out = "parsedMRP2.fasta")
Hi. I would like to ask a similar question: I have a Fasta alignment, that I have created a object with: myotis <- read.dna("epfu.fas",format= "fasta") I am interested in using bootstrap to estimate confidence intervals, if it is actually possible. I would like to resample with replacement sequences from that alignment, and then using the function tajimaD from the package strataG or the function tajima.test$D from the package pegas.
Log in to answer this question.
-
BLAST options error: File homo_ref.faa does not exist
written by logbio •
-
Writing a fasta file in R
written by margo •
-
Fasta Files - How to search for peptide sequences?
written by janainamace •
-
write.fasta command - losing attributes/annotations
written by c.scott •
-
Amino Acid Substitution Loop Generates Incorrect Output: No Mutations in the Output, but the same s…
written by jguevar •
-
Reading FASTA file with multiple sequences
written by Seot •
-
fasta file loading error
written by tcf.hcdg
-
Subsetting a fasta file using seqinr in R
written by henrikkjeldal •
-
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 •
If you need an alternative then using
faSomeRecordsfrom Jim Kent's utilities is a great/fast option. Linux version is linked but macOS available as well. After downloading do not forget to add execute permissions (chmod a+x faSomeRecords).