Dear BIostar community
I have around 50,000 sequence identifier in a text file and around 1Million sequences in a large fasta file.
I want to get the sequence of these 50,000 from the fasta file.
I checked the questions and answers already in the community pages. But most of them are using perl,linux,python.
How can it be done in R(fasomerecord/grep/something else)?
Any suggestion?
Thanks
2 answers
If you really want to do this in R, you could use the seqinr package. You end up loading the fasta file into memory, subsetting the results according to the names() accessor, and then writing the results of that to a file. Presumably you could alternatively pass only the names you want to write.fasta().
•
written
by
Devon Ryan
Log in to answer this question.
-
Rename fasta headers using a variable already present in the header?
written by FreshBio •
-
Batch Renaming Files
written by Tastulek •
-
Calculating p-distance for multiple fast files
written by Miguy •
-
How to remove duplicated sequences from multi fasta file?
written by tcf.hcdg
-
Extracting specific sequences from a big fasta file using ids of the sequences to be excluded
written by hasche89 •
-
How to retrive unique IDs from the txt file?
written by tcf.hcdg
-
id in text file for sequence retrival from multifasta
written by tcf.hcdg
-
sequence extraction from multifast file using IDs in another text file
written by tcf.hcdg
-
Biostar Participation As A Bioinformatics Course Assignment?
written by Obi Griffith
-
Appending Identifiers To The Fasta Header: Multiple Files
written by himanshu17 •
You want this done with
grep, but not in linux?I just want to know it this is possible in R?
Yes, it is (please see this post: How To Extract Multiple Fasta Sequences At A Time From A File Containing Sequences Ids Using R)
samtools faidx is pretty handy.