This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Batch rename protein fasta headers

Hey guys,

I have tons of protein multi-fasta files and I would like to append the name of the file to the fasta-headers. For example, for a input file one.txt with the headers

>1
ATGC...
>2
ATGCAT...

I would like to have the output

>one_1
ATGC...
>one_2
ATGCAT...

I use bbrename for DNA sequences, but doesn't work for protein files. Thanks!

sequence
• 3,249 views • link
updated by Joe Joe • written by genomes_and_MGEs • genomes_and_MGEs

your example files also don't really looks to be protein either ...

• 0 views • link
written by lieven.sterck lieven.sterck

5 answers 

 awk '/^>/ {printf(">%s_%s\n",substr(FILENAME,1,length(FILENAME)-3),substr($1,2));next;} {print}' *.txt
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

As easy as:

for file in /path/to/files/*.fasta ; do
    sed "s/>/>$(basename $file .fasta)/gi" $file
done

You can tweak it if you want to keep the extension or whatever...

• 0 views • link
written by Joe Joe

You could use simple python script like:

from Bio import SeqIO

with open("one.txt", "r") as input:
    with open("output_filename.fasta", "w") as output:
        for record in SeqIO.parse(input, "fasta"):
            record.id = f"one_{record.id}"
            record.description = ""
            SeqIO.write(record, output, "fasta")

This require installation of biopython.

• 0 views • link
written by Hood • Hood

Thanks for the reply. The thing is that I have multiple files to use as input. So, I guess using a loop to create a renamed output for each file would be better. Do you think you can help me with this? Usually for DNA sequences, I use

for F in *.fasta; do N=$(basename $F .fasta) ; bbrename.sh in=$F out=${N}_mod.fasta prefix=$F addprefix=t ; done

I need to find an alternative that works with multi-fasta protein files as input

• 0 views • link
written by genomes_and_MGEs • genomes_and_MGEs 
for f in `ls *.fasta | sed 's/.fasta//g' `; do sed "/^>/ s/.*/&_$f/" "$f.fasta" >  "$f_new.fasta" ; done

.

• 0 views • link
written by Fatima Fatima

Using seqkit

seqkit replace -p '(.+)' -r 'one_$1' Filename.fasta

• 0 views • link
written by lakhujanivijay lakhujanivijay

Log in to answer this question.

More posts like this