Add Information to Protein Fasta Headers
Hi, I have protein fasta file whose headers look like '>evm.model.chr.9.52'. There are almost 30k+ proteins. I have performed functional annotations and also added every information to gene structure we get from EVM. The thing is, in that files I had columns so I basically merged information. Now, I al performing some analysis and I want to add atleast protein name or even GO term in fasta header so it would make things alot easier for me. I want something like;
>evm.model.chr.9.52 GO:1234678
Can I do it with grep? Or Seqkit or any other method please? Any help would be highly appreciated. Thanks!!!
• 1,753 views
•
link
written
by
ahmadjoyyia •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to limit fasta header to 40 characters?
written by fishingline37 •
-
To find amino acid substitution from fasta alignment
written by shwetamgr1 •
-
Transfer annotation from fasta headers to associated gff
written by EJB •
-
Find defence related and disease resistant genes in Rice genome
written by vedikaa96 •
-
GenSAS merged annotation GFF3 file to NCBI .sqn file for NCBI genome submission
written by zack.saud •
-
fasta nucleotide to fasta protein
written by sg927357 •
-
Splitting a fasta file based on the number of headers.
written by sicat.paolo20 •
-
gene functional analysis using GO and KO numbers
written by blooming.daisy333 •
-
How to get list of all GO terms in each ontology
written by sayonidas03
-
Functional Annotation Clustering On The Basis Of Go Terms
written by figo
There is not enough information about where these GO terms or protein names are going to come from. Do you have a separate file which has this information?
Thank you for the reply! I have GO terms in separated tab delimited file.
You may want to post example lines from that file.
I already gave you example how my protein fasta header looks like. I had done functional annotation from many sources and made a custom file which looks like this;
That is a tab-delimited file and has proper column names. Like 5th column is of GO_Description. Lets say I want to copy GO_Description or GO from this file to my protein fasta file with matching ID, which in my case he would be 'chr.1.1128'. How can I do this? Thanks!!
GenoMax Hi.. I can see you have edited my reply and it is something I am looking for. But can you please tell me how did you do it? And how can I do it please?