This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Is there a way to change only the header of a fasta file?

Dear All,

I would like to change the header of the fasta file as follows, but I don't know how. The actual file is so long that I need a command. Could you tell me the command to change only the header.

Before

>TRINITY_DN22056_c0_g1_i1 len=300 path=[0:0-299] GTCCATTTGACTGATTCCCATTCCAAAAAAAACCACTCAACCAGGACTGACGCATAAAACCGTCAAATCGACGACGCAGGATAAGGTTCGATTATGACTTGTTACGGTCTCAAAGGCGCGTTTTTACACAATTATG >TRINITY_DN30979_c0_g1_i1 len=300 path=[0:0-299] TGACGAATGAGAACAAGGGCACGCAGTTCGTTTTATGTCAATTGAACTGCTGGTCACTAACCAGAATTCGTATTCCCGTTGGACGTCAATTTGCATAAAAAATGTAAAAGAGAAAATACGATCATAAAAAGTGGAAA

After

>AbTRINITY_DN22056_c0_g1_i1 len=300 path=[0:0-299] GTCCATTTGACTGATTCCCATTCCAAAAAAAACCACTCAACCAGGACTGACGCATAAAACCGTCAAATCGTACGACGCAGGATAAGGTTCGATTATGACTTGTTACGGTCTCAAAGGCGCGTTTTTACACAATTATG >AbTRINITY_DN30979_c0_g1_i1 len=300 path=[0:0-299] TGACGAATGAGAACAAGGGCACGCAGTTCGTTTTATGTCAATTGAACTGCTGGTCACTAACCAGAATTCGTATTCCCGTTGGACGTCAATTTGCATAAAAAATGTAAAAGAGAAAATACGATCATAAAAAGTGGAAA

This question may not exactly be a bioinformatics question. If this was inappropriate, I apologize. Thank you for your help.

awk seqkit fasta
• 1,987 views • link
updated by liorglic liorglic • written by Riku • Riku

A simple Linux command can do it:

sed 's/>/>Ab/' your.fasta > your.new.fasta

or modify in place:

sed -i 's/>/>Ab/' your.fasta
• 0 views • link
written by liorglic liorglic

Good point, thanks.

• 0 views • link
written by liorglic liorglic

2 answers 

sed 's/^>/>Ab/'
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

So simple! I've learned a lot. Thank you very much!

• 0 views • link
written by Riku • Riku 
$ bioawk -c fastx '{print ">Ab"$name;print $seq}' test.fa
• 0 views • link
written by cpad0112 cpad0112

Log in to answer this question.

More posts like this