how to replace fasta header

Question

[Nelo]

Hi!

I have two files, ba.fasta and ba.header

ba.fasta                                                                               ba.header
>bas12                                                                                >seq1
MDNAVGYH..                                                                     MDNAVGYH..
>bas13                                                                                >seq2
MLSRTEQR..                                                                      MLSRTEQR..

So basically I wanted to replace the header of ba.fasta with that of ba.header while keeping the sequences

Any easy command for this please

Tnks in advance

[cpad0112]

$ tail -n+1 *.fa
==> test1.fa <==
>bas12
MDNAVGYH
>bas13
MLSRTEQR
>seq
MLR

==> test2.fa <==
>seq1
MDNAVGYH
>seq2
MLSRTEQR

output:

$ seqkit --quiet replace -K  -p '(.*)' -r '{kv}' -k  <(paste <(seqkit seq -n test1.fa) <(seqkit seq -n test2.fa)| awk '$2!="" {print}') test1.fa
>seq1
MDNAVGYH
>seq2
MLSRTEQR
>seq
MLR

Answer 1

[Pierre Lindenbaum]

awk '/^>/{printf(">seq%d\n",++N);next;} {print;}'

[Nelo]

ba.fasta ba.header

Bol0453 >Bol.11 MKDIRGHTWQP... MKDIRGHTWQP... Bol0358 >Bol.75 MKSREVBTW.... MKSREVBTW....

What if situation is something like this. I want the header of file ba.header in the place of ba.fasta header.

[Pierre Lindenbaum]

linearize the fasta and use paste to merge the data.

[Nelo]

First of all thank you for responding but I am not able to get you

[Nelo]

Ok... I got it !!!