This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Renaming Fasta ID's according to their number and on specific string length

Hi all,

I have a question about renaming my sequences in fasta format. I have 36,008 nucl. sequences and I need (want) rename them according to the following form. For every sequence, after > mark should follow string E_nip_trans_ (for all sequences) and after this I want place only digit string with length five, growing by one (from 00001 to 36008). But I don't have any idea how to do it.

>E_nip_trans_00001
>E_nip_trans_00002
>E_nip_trans_00003
> ...
>E_nip_trans_00010
> 
> ...
>E_nip_trans_00111  
> 
> ...
>E_nip_trans_36008

Thank you for reply!

sequence
• 2,326 views • link
updated by Ram Ram • written by Jiri • Jiri

1 answer 

awk '/^>/ {printf("%05i\n",i+1); i+=1; next;} {print $0;}' /path/to/your/fasta/file
• 0 views • link
written by cschu181 cschu181

Thank you very much. It's work fine.

• 0 views • link
written by Jiri • Jiri

If this answer was helpful it is appropriate to upvote it, and if this answer resolved your question completely you can 'accept' the answer, as such marking your question as solved.

• 0 views • link
written by WouterDeCoster WouterDeCoster

Log in to answer this question.

More posts like this