Extracting fasta sequence from CDS matching with sequence ID
I am trying to extract matching sequences from a header list file. My fasta sequence looks like this :
>TRINITY_DN0_c0_g1_i1.p1 TRINITY_DN0_c0_g1~~TRINITY_DN0_c0_g1_i1.p1 ORF type:complete len:141 (-),score=18.99,CCNL1_PONAB|68.750|6.95e-55,Cyclin_N|PF00134.23|6.5e-09 TRINITY_DN0_c0_g1_i1:85-507(-)
ATG <sequences in between>................TAG
I have a separate header file <header.list> that that has list of headers like TRINITY_DN0_c0_g1_i1
I now want to extract fasta file along with header that matches my header file.
The problem that i have is : those fasta headers are long enough and does not match with my header.list ( as these headers are just a part of headers in fasta file).
Does anyone have any idea how to overcome this problem
• 962 views
•
link
updated
by
Ram
•
written
by
metallicasster •
•
written
by
metallicasster •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to create a tax_id_map file for getting scientific names in BLASTN output
written by Maria •
-
Renaming fasta headers based on a list
written by Bertrand •
-
TransDecoder output row names different than Trinity.fasta
written by bioinfo_cat •
-
Using transdecoder to filter a fasta file with preserving whole transcripts
written by poecile.pal •
-
Help with extracting multiple sequences from a fasta file with a list of Ids (four counting per lin…
written by ahmed_bio82 •
-
Samtools extraction output is empty
written by spidcock01 •
-
extract headers of a file using a second file containing a list of IDs
written by paraskevopou •
-
extract fasta sequence using fasta header with python 3
written by ddowlin
-
Extract sequences based on keywords in headers (no perl or python code)
written by Crystal •
-
perl code to extract sequences from multi-line fasta works on all test files but not on research fi…
written by kdiaz17 •
This question has been answered multiple times on the forum. Please search the forum for "fasta extract id"