This is a test version of Biostars. For the public version, visit https://www.biostars.org.
fasta header with leading zeros

Hi, I would like to rename the sequences in my file to have leading zeros in front of the header name using unix or perl based on the number of sequences in the file. I am unable to open these files using a text editor since they are too large. Example.

>gmp-15 
>gmp-126
>gmp-127890 to
>gmp-000015
>gmp-000126
>gmp-127890

Thanks

sequence fasta
• 2,604 views • link
updated by 5heikki 5heikki • written by rrsowmya • rrsowmya

1 answer 


awk -F "-" '{if(/^>/)printf("%s-%06d\n",$1,$2); else print $0}' file.fasta > out.fasta

• 0 views • link
written by 5heikki 5heikki

That worked just great! Thanks so much! I really appreciate your help!

• 0 views • link
written by rrsowmya • rrsowmya

so validate 5heikiki's answer by clicking on the hand on the left please.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Please accept @5heikki's answer by clicking on the green check mark on the left of the answer.

• 1 views • link
written by Ram Ram

Log in to answer this question.

More posts like this