This is a test version of Biostars. For the public version, visit https://www.biostars.org.
sequence renaming

I have a fasta file how to rename the sequence according to the name of the file, for example, the sequence name in the gene.fasta file is >1, >2, renamed >gene1, gene2. Thank you very much for telling me

renaming
• 1,473 views • link
updated by cpad0112 cpad0112 • written by zhichusun • zhichusun 
$ sed -r  '/^>/ s/^./>gene_/' test.fa
• 0 views • link
written by cpad0112 cpad0112

2 answers

Using Linux (any flavour), Unix (OSX)

perl -pi -e 's/^(>[0-9])+.*/\1gene/g' myfile.fa

myfile.fa

>1
AGTC
>2
AGTC
>3
AGTC

output

>1gene
AGTC
>2gene
AGTC
>3gene
AGTC

Its called perl pie (one liner) and is extremely quick at handling massive files. Happy to explain the reg-ex if needed. It will change the file in situ so there's not needed to pipe it, or make a copy, the -i takes care of that.

• 0 views • link
written by M__ • M__

This can be done with seqtk:

seqtk rename gene.fasta gene > renamed.fasta
• 0 views • link
written by Mensur Dlakic Mensur Dlakic

Also, a simple replace command with perl:

perl -p -e 's/\>/\>gene/g' gene.fasta > renamed.fasta
• 0 views • link
written by Mensur Dlakic Mensur Dlakic

Log in to answer this question.

More posts like this