This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Add Specific Word To Fasta Header

I have more than 5000 fasta sequence in a file and want to add a word , for instance phosphate, to header of all sequence. please tell me a PERL solution for that.

fasta
• 19,123 views • link
updated by O.rka O.rka • written by Palu • Palu

4 answers

inplace: 

    perl -pi -e "s/^>/>phosphate-/g" your.fasta

or new file:

    perl -p -e "s/^>/>phosphate-/g" your.fasta > phosphate.fasta

to add it to the end, use this regexp

    's/^(>.*)$/$1-phosphate/g'
• 0 views • link
written by brentp brentp

THANKS brentp, but i want to in the last of my header..is there any trick

• 0 views • link
written by Palu • Palu

@palu I edited my answer, see the last line.

• 0 views • link
written by brentp brentp

thank you very much sir. for any layman person like me. final code will be like that

perl -p -e "s/^(>.*)$/$1-phosphate/g" your.fasta > phosphate.fasta

• 0 views • link
written by Palu • Palu

except you should use single quotes: perl -p -e 's/^(>.*)$/$1-phosphate/g' in.fasta > out.fasta

• 0 views • link
written by brentp brentp

palu, if you like Brent's answer the best, you should select it as such (hover over the votes to do that).

• 0 views • link
written by Neilfws Neilfws

@newlife thanks i do that

• 0 views • link
written by Palu • Palu

An easy way with sed:

sed 's/>.*/&_phosphate/' foo.in >bar.out
• 153 views • link
written by Daniel Daniel

A faster option, from the BBMap package:

bbrename.sh in=file.fasta out=renamed.fasta prefix=phosphate addprefix=t

• 0 views • link
written by Brian Bushnell Brian Bushnell

Is there a way to get this to work with protein fasta?

• 0 views • link
written by O.rka O.rka

I know there are lots of option and it can be easily done with many unix one liners, but here is another alternative (my favorite).

bioawk -c fastx '{ print ">PREFIX"$name; $seq }' input.fasta
bioawk -c fastx '{ print ">"$name"|SUFFIX"; $seq }' input.fasta
• 0 views • link
written by arnstrm arnstrm

Hi, I'm not sure I understand the (bio)awk syntax, but your command was not working for me (did not print sequences)...I put there a new line instead of a semicolon:

 bioawk -c fastx '{ print ">PREFIX" $name "\n" $seq }' input.txt >outupt.txt

which seems to work. Anyway thanks for pointing me towards the solution.

• 1 views • link
written by al-ash al-ash

Log in to answer this question.

More posts like this