Change Fasta Id'S
Hello,
I've got multiple blast hits in a fasta file. But I want the description, instead of accessions. So when I make a tree, I get nice names.
For example:
>gi|255587543|ref|XM_002534260.1| Ricinus communis Protein kinase APK1B, chloroplast precursor, putative, mRNA
To:
>Ricinus_communis_APK1B
I know this can be done by hand, but there's just too many of them!
Hope you can help me.
2 answers
duplicate of
• 204 views
•
link
written
by
Pierre Lindenbaum
Fire up perl: If starts with > , regex to extract between pipe and comma, print >[regex] else print.
Further info: http://learn.perl.org/
Log in to answer this question.
More posts like this
-
BWA with multiple references
written by Aurélien Guy-Duché •
-
Running Batch Files From Command Line
written by bioinfo
-
How To Add An Attribute To A Sequence Header (Of Fasta) Based On A Match Off Another Fasta?
written by microbeatic
-
Heebo-->Gene Symbol Conversion, Or A Better Way To Get Geo Annotation Files
written by ellenmseb •
-
How To Retrieve A Gff With Sequence Identifier Different Than Gb With Bioperl
written by fbrundu
-
Getting 'Bus Error: 10' When Running The Prody Package On A Multiple Sequence Alignment
written by Harpalss
-
Creating Table From Fasta File In Python
written by Peiska •
-
How To Join Text Files By A Certain Column (Accession Id)
written by 2011101101 •
-
Get Fasta File With Protein Sequences Given Entrez Gene Ids
written by Stephen
-
Bioperl Has Different Behaviours In Parsing Blast And Blast+ Result
written by Yixf •
duplicate of
Any Script to parse fasta headers?
With what do you want help? This task and many like it is quite simple using e.g. a Perl or shell script. Is that what you would like to learn? Right now, the tone of the question suggests that you just want someone to do it for you. That does not help you, in the long run or the people with the answer, in the short term.