Thanks a lot
• 0 views
•
link
I have a header like this-
>sp|A0A0H2Z7X0|TPBB_PSEAB Diguanylate cyclase TpbB OS=Pseudomonas aeruginosa (strain UCBPP-PA14) OX=208963 GN=tpbB PE=1 SV=1
MNRRRRYTGSNPSLRRVLYRAHLGVALVAVFTAGLAVTLVGLLTLRAYADPNQQLIARSI
SYTVEAAVVFGDAQAAEESLALIASSEEVSSAIVYDRQGQPLASWHRESTGPLHLLEQQL
And I want my header like this
>A0A0H2Z7X0
MNRRRRYTGSNPSLRRVLYRAHLGVALVAVFTAGLAVTLVGLLTLRAYADPNQQLIARSI
SYTVEAAVVFGDAQAAEESLALIASSEEVSSAIVYDRQGQPLASWHRESTGPLHLLEQQL
Thank you!
Under the assumption that everything starts with sp:
sed 's/>sp|/>/;s/|.*//' YourSequences.fa > YourNewSequences.fa
Thanks a lot
You can use awk too for this propose-
cat yourdata.fasta
>sp|A0A0H2Z7X0|TPBB_PSEAB Diguanylate cyclase TpbB OS=Pseudomonas aeruginosa (strain UCBPP-PA14) OX=208963 GN=tpbB PE=1 SV=1
MNRRRRYTGSNPSLRRVLYRAHLGVALVAVFTAGLAVTLVGLLTLRAYADPNQQLIARSI
SYTVEAAVVFGDAQAAEESLALIASSEEVSSAIVYDRQGQPLASWHRESTGPLHLLEQQL
cat yourdata.fasta |awk 'BEGIN { FS=OFS="|" } /^>/ { print ">" $2; next }1'
>A0A0H2Z7X0
MNRRRRYTGSNPSLRRVLYRAHLGVALVAVFTAGLAVTLVGLLTLRAYADPNQQLIARSI
SYTVEAAVVFGDAQAAEESLALIASSEEVSSAIVYDRQGQPLASWHRESTGPLHLLEQQL
Log in to answer this question.