This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Fasta header contain multiple string

I have large fasta file. As you see below there are > sign present in some fasta header like

>exon2_ENST00000218032|>exon2_ENST00000218032
>exon17_ENST00000253024|>exon17_ENST00000253024

I want to remove the >sign from the header sequence, after remove the header is then look like this

>exon2_ENST00000218032|exon2_ENST00000218032
>exon17_ENST00000253024|exon17_ENST00000253024

It is an actual fasta header that has only one > sign not more than 1.

>exon2_ENST00000218032|>exon2_ENST00000218032
GAAGCTTTTGGTCTATATTGTTAATTGCCATTGCTGTAAATCTTAAAATGAATGAATAAAAATGTTTCATTTTACAAAAAACATGTTCCTTCAGTCGTCAATGCTGACCTGCATTTTCCTGCTAATATCTGGTTCCTGTGAGTTATGCGC
>exon1_ENST00000218032|>exon1_ENST00000218032
GCTGCTGCAAGTTACGGAATGAAAAATTAGAACAACAGAAACATGGTTTCTCTTCTCGGCCACCTCCTGCATAGAGGGTACCATTCTGC
>exon1_ENST00000218032|exon2_ENST00000218032
AAGCTTTTGGTCTATATTGTTAATTGCCATTGCTGTAAATCTTAAAATGAATGAATAAAAATGTTTCATTTTACAAGTTTCTCTTCTCGGCCACCTCCTGCATAGAGGGTACCATTCTGCGCTGCTGCAAGTTACGGAATGAAAAATTAG
>exon17_ENST00000253024|>exon17_ENST00000253024
TCCCTGGTGGCCCCATCCCCCAGTTCCTCACGATATGGTTTTTACTTCTGTGGATTTAATAAAAACTTCACCAGTTACAAGGCAGACGTGCAGTCCATCATCGGCCTGCAGCGCTTCTTCGAGACGCGCATGAACGAGGCCTTCGGTGAC
>exon16_ENST00000253024|>exon16_ENST00000253024
TACAGCTCCCCACAGGAGTTTGCCCAGGATGTGGGCCGCATGTTCAAGCAATTCAACAAGTTAACTGAGACCAGCCCGGTGGCACCCTGGATCTGACCCTGATCCGTGCCCGCCTCCAGGAGAAGTTGTCACCTCCC

So please can anyone tell me how I remove the multiple > sign in my fasta header. Thanks in advance

fasta
• 1,824 views • link
updated by cpad0112 cpad0112 • written by harry • harry

1 answer

Removing all > and then inserting a new one at the beginning.

sed '/^>/ { s/>//g; s/^/>/ }'  seqs.fasta > result.fasta

Ref:

• 0 views • link
written by shenwei356 shenwei356

Thanks it working fine can you just help me how do i remove duplicate header with sequence as you can see ENST00000368129_1|ENST00000368129_3, It contain 2 sequence: I tried this command to remove duplicate sequence:

paste -d $'\t' - - <fastaFileWithNoLinebreaksInSeq | sort -t $'\t' -uk1,1 | awk 'BEGIN{FS="\t";OFS="\n"}{print $1,$2}'

but it does not give right output: It give output like this

>_10|_10
AAAGCAAAACACAGAATTGACGGAAAGACTTACGTTATTAAACGTGTTAAATATAATAACGAGTTTGGCATGGATTTTAAAGAAATAGAATTAATTGGCTCAGGTGGATTTGGCCAAGTTTTC
>_10|_11
ATTATGATCCTGAGACCAGTGATGATTCTCTTGAGAGCAGTGATTATGATCCTGAGAACAGCAAAAATAGTTCAAGGTTTGGCATGGATTTTAAAGAAATAGAATTAATTGGCTCAGGTGGATTTGGCCAAGTTTTCAAAGCAAAACACA
>_10|_12
TTTGGAACTCTTTGAACAAATAACAAAAGGGGTGGATTATATACATTCAAAAAAATTAATTCATAGAGATCTTAAGGTTTGGCATGGATTTTAAAGAAATAGAATTAATTGGCTCAGGTGGATTTGGCCAAGTTTTCAAAGCAAAACACA
>_10|_13
ACTTGTAACATCTCTGAAAAATGATGGAAAGCGAACAAGGAGTAAGGGAACTTTGCGATACATGAGCCCAGAACAGGTTTGGCATGGATTTTAAAGAAATAGAATTAATTGGCTCAGGTGGATTTGGCCAAGTTTTCAAAGCAAAACACA
>_10|_14
AACGAATTTCTTCGCAAGACTATGGAAAGGAAGTGGACCTCTACGCTTTGGGGCTAATTCTTGCTGAACTTCTTCATGTTTGGCATGGATTTTAAAGAAATAGAATTAATTGGCTCAGGTGGATTTGGCCAAGTTTTCAAAGCAAAACAC

I don't want to lose any header information, so can you tell me what is wrong with this command.

>ENST00000368129_1|ENST00000368129_1
GAAAGTCCACAGAGGAGTTTAAAGCAGCCATGCCAAAAGTGCACTTGCACTCTAAGGAAGCTGAGGTGGGGGAGGCGG
>ENST00000368129_1|ENST00000368129_3
CATATTATTTGACCTAAGTGACAACACTGGGAAAATGGAAGTACTGGGGGTTAGAAACGAGGACACAATGAAATGTTGCACTTGCACTCTAAGGAAGCTGAGGTGGGGGAGGCGGGAAAGTCCACAGAGGAGTTTAAAGCAGCCATGCCA
>ENST00000368129_1|ENST00000368129_2
TGAAACCCCGAAGATCAACACGCTTCAAACTCAGCCCCTTGGAACAATTGTGAATGGTTTGTTTGTAGTCCAGAAGTGCACTTGCACTCTAAGGAAGCTGAGGTGGGGGAGGCGGGAAAGTCCACAGAGGAGTTTAAAGCAGCCATGCCA
>ENST00000368129_1|ENST00000368129_3
CATATTATTTGACCTAAGTGACAACACTGGGAAAATGGAAGTACTGGGGGTTAGAAACGAGGACACAATGAAATGTTGCACTTGCACTCTAAGGAAGCTGAGGTGGGGGAGGCGGGAAAGTCCACAGAGGAGTTTAAAGCAGCCATGCCA
>ENST00000368130_6|ENST00000368130_6
TTAGGATAGAATAATTGCTGGATAAACAAATTCAGAATATCAACAGATGATCACAATAAACATCTGTTTCTCATTCAGTTATTAAGGCCAAAAAAAAAACATAGAGAAGTAAAAAGGACCAATTCAAGCCAACTGGTCTAAGCAGCATTT

Thanks in advance for your instant reply.

• 0 views • link
written by harry • harry

There is no problem with the command. Said command doesn't trim any thing. Check your input if it has any spaces in headers. Please post this as separate post with headers from fastaFileWithNoLinebreaksInSeq.

• 0 views • link
written by cpad0112 cpad0112

based on OP format, it can be further simplified:

$ sed '/^>/ s/>//2' input.fasta

2g in case of more > in header.

• 0 views • link
written by cpad0112 cpad0112

Log in to answer this question.

More posts like this
  • remove the fasta sequence which doesn't contain ATGC
    written by harry • harry
    I have a fasta sequence in which some sequence doesn't contain "ATGC" like below so I want to remove all fasta sequences from my file …
  • grep command for fasta header
    written by harry • harry
    I used this command:-- grep -Fw -A 1 -f header.txt test.fa >test_result.fa But it extracts only 1 header, not the whole which are present in …
  • I want to reverse the fasta sequence
    written by harry • harry
    I have a large file and I want to reverse the fasta sequence. Please can anyone tell how to do it Thanks in advance
  • Remove fasta sequence on the basis of header name
    written by harry • harry
    >exon18_ENST00000194900|exon19_ENST00000194900 GCTATCAAGAGACTGCAGCAAGCACAACTTTACCCCATTGCCATTTTCATCAAGCCCAAGTCCATTGAAGCCCTTATTACTACCCGGCCTCGACGTGATAATGAGGTGGATGGACAAGACTACCACTTTGTGGTGTCCCGAGAACAAATG >exon18_ENST00000194900|exon21_ENST00000194900 TCTAACAGAGGAAGAACTGTATTATCAAACCTTATGGTCTACCCATGAATAAATAAAATATTTGTTCAAGCACAAATACTACCCGGCCTCGACGTGATAATGAGGTGGATGGACAAGACTACCACTTTGTGGTGTCCCGAGAACAAATGG >exon18_ENST00000194900|exon20_ENST00000194900 CAGACATATGAACAAGCAAATAAGATCTATGACAAAGCCATGAAACTGGAGCAGGAATTTGGAGAGTACTTTACAGTACTACCCGGCCTCGACGTGATAATGAGGTGGATGGACAAGACTACCACTTTGTGGTGTCCCGAGAACAAATGG >exon18_ENST00000194900|exon21_ENST00000194900 CTCTAACAGAGGAAGAACTGTATTATCAAACCTTATGGTCTACCCATGAATAAATAAAATATTTGTTCAAGCACAAATACTACCCGGCCTCGACGTGATAATGAGGTGGATGGACAAGACTACCACTTTGTGGTGTCCCGAGAACAAATG >exon17_ENST00000194900|>exon17_ENST00000194900 ATGATGACCTGATCTCCGAATTTCCACATAAATTTGGATCCTGTGTGCCACTCACTATGCAAGGCCTGTGATCATCCTGGGCCCAATGAAGGACCGAGTCA >exon19_ENST00000194900|exon21_ENST00000194900 TCTAACAGAGGAAGAACTGTATTATCAAACCTTATGGTCTACCCATGAATAAATAAAATATTTGTTCAAGCACAAAGCAAGCACTGCATCTTAGATGTTTCCGGCAATGCTATCAAGAGACTGCAGCAAGCACAACTTTACCCCATTGCC >exon19_ENST00000194900|exon20_ENST00000194900 CAGACATATGAACAAGCAAATAAGATCTATGACAAAGCCATGAAACTGGAGCAGGAATTTGGAGAGTACTTTACAGGCAAGCACTGCATCTTAGATGTTTCCGGCAATGCTATCAAGAGACTGCAGCAAGCACAACTTTACCCCATTGCC >exon19_ENST00000194900|exon21_ENST00000194900 CTCTAACAGAGGAAGAACTGTATTATCAAACCTTATGGTCTACCCATGAATAAATAAAATATTTGTTCAAGCACAAAGCAAGCACTGCATCTTAGATGTTTCCGGCAATGCTATCAAGAGACTGCAGCAAGCACAACTTTACCCCATTGC In this example of fasta sequence, you see there …
  • Want to reverse the fasta sequence
    written by harry • harry
    >ENST00000014935 TAGGTTTTTAATTACCTTTATTTTATTTTGCCAAACATACCTGGGAATACCTTTTATTTTTTTTTTACCTTGGGGTGATGGTTCCAAACCATAAATGTGATTATAGTTAACACATGACCCTTCTAGCGTCCCAGCCAGTGTTTTTCCTGACCTCTCTTCTTTGGAGAGGAGGATGGAAGGGAGGGGTCCGGCATGCTGCTGGCATTTTGCTGTGTCCTGCAGCCCCTTTCCGGGACACCTGGGTTCACACAGCTTTTTAGCTTACATAACTGGTGCAGATTTTCTGTGTGGAGATGTTGCCTTGACCAGCCTTGGCTGGACTTTACCAGGCATGCAGAAGCCTGTACCAACACAGACTACAGCACCCAGGAGGTGCGAGTGTGGCTGCTCAGCGGTTATAACAGGCCTGACTGCATTGTTCACCGGATTATAATGAGCCAAAATGTTTCCCGGTGTTTGCTGGTTTCAGGGAAGGAGTTTGATATAGCAGATTAACCACCCTCCTTGTAGCTATTGGGGCTTAATGGTTTCCTGGTGATTCTTACCAATCCACAATAAACATGGCCCATTGGCATATCTGCTGCACAAGTGTCCTATCTCACCAATCTGGGTTTTTGTTCTCAGTAACTTTCCTTCTTGTCATACAACATCTTCATTCCTCTTTCTGAACCCTCCCTTCCCCTACCCCAACCCAGAGCCCACTTTGTCTCCACTCCTGATACTACACTACCTGGCAGGTGGCATGAGTGCAGGGCCCCTGGCTTCCTCTCCTAATCTAGGCACAAGCCCAACCAAAGAACAAGAGCCAAATCAAACAAGGCAGGCAGGGGTGGACTACAGTCACAGGGCAACTATAGTTGAAGCCCCCCAGCCCCAGGGCTGGATGGACGGGGGAGGCTGGGGTTTAAGTCCCAAAAGGCAGCAGGCCCTGGGGGGGTAGGGGGACGCTCAGGCAGCAGGGCACAGCTGAGGGGACAGGAGTGATAGCAGCAACAGAACAGTGAGGCTGAGAGGCTGGACGCTGTGCGCCTGGCTCAGCTTCAGCTCCACCTCCACGGGGTAGTGGTCACTGATGTTGAGGGC >ENST00000014944 TCCTCCTCGGTGAGCTGGAAGCTCGTGGGGAAGTCAAAGGCAGCCGCAGTGTGCAGCAGACTCCGGCAGCGCTCCCCGTGCAGCACGACGCGGTCATAGGTGCAGTGGGTGCTGGCCCGCACTGTGGTGTCCTCCCCATCGGCAATCACCCAGTGGAAGCCTGGCTCAGTCCGCAGCTCCAGCTTGTCCAGGCGCTTTTTGGTCAGTGAAGCGCAGTCAGCATTGAAGTCCCCAAGCAGGATCACGTC Likewise, I have a huge fasta sequence and I want to reverse the sequence of those fasta sequences like below:- >ENST00000014935 …
  • trimming of fasta file
    written by harry • harry
    As I have multiple fasta sequences like see below: >read1 CAGGTCTGGCTGGATGAAGGGCACGGCATAGGTCTGACCTGCCAGGGAGTGCTGCATCCTCACAGGAGTCATGGTGCTGCTGAAGATGTCTCCAGAGACCTTCTGCAGGTACTGCAGGGCATCCGCCATCTGCTGGACGGCCTCCTCTCGCCG >read2 CATCTGCGGAGGCTGCCGTGACGTAGGGTATGGGCCTAAATAGGCCATTGTGAGTCATGAGCTTGGTCTGTAGAGGCTGACTGGAGAAAGTTCTGGGCCTGGAGAGGCTGCCGGGAGGTAGGAGTGGTGAGGTCGACTTGAGAAAGTTCAGGGCCTGGAGAGAAGGCTGGGAGGCAGGAGCTGGGTCTAAAGAGGCCATTGTAACGATGGAGCTGTGCCTGTGGAGGCTGTTGTGAGGCAGTAGCCT >read3 TTGAGGTGGGAGGATCGCTTCAGCCTGGAAGGTTGAGGCTGCAGTCAGCTGCGATAGCACTACTACACTCCAGCCTTGGACAACAGAGGGAGACCTTTCGCTGTCACCCCTCTAGAATCCACGTATACGAAAATTCCAAATGTTAGTTGGGCATAGTGGCAAGCACCTGTAGTCTCAGCCACGTGGGAGG These are in different lengths so I want to isolate …
  • Script for making exon file
    written by harry • harry
    Hi everyone, I have an exon.fasta file which looks like --- >ENSE00002234944 TTAACTTGCCGTCAGCCTTTTCTTTGACCTCTTCTTTCTGTTCATGTGTATTTGCTGTCTCTTAGCCCAGACTTCCCGTGTCCTTTCCACCGGGCCTTTGAGAGGTCACAGGGTCTTGATGCTGTGGTCTTCATCTGCAGGTGTCTGACTTCCAGCAACTGC TGGCCTGTGCCAGGGTGCAAGCTGAGCACTGGAGTGGAGTTTTCCTGTGGAGAGGAGCCATGCCTAGAGTGGGATGGGCCATTGTTCATCTTCTGGCCCCTGTTGTCTGCATGTAACTTAATACCACAACCAGGCATAGGGGAAAGATTGGAGGAAAGATGA GTGAGAGCATCAACTTCTCTCACAACCTAGGCCA >ENSE00001746346 CAGGGCCATCAGGCACCAAAGGGATTCTGCCAGCATAGTGCTCCTGGACCAGTGATACACCCGGCACCCTGTCCTGGACACGCTGTTGGCCTGGATCTGAGCCCTGGTGGAGGTCAAAGCCACCTTTGGTTCTGCCATTGCTGCTGTGTGGAA >ENSE00002312635 CAGGGCCATCAGGCACCAAAGGGATTCTGCCAGCATAGTGCTCCTGGACCAGTGATACACCCGGCACCCTGTCCTGGACACGCTGTTGGCCTGGATCTGAGCCCTGGTGGAGGTCAAAGCCACCTTTGGTTCTGCCATTGCTGCTGTGTGGAAGTTCACTCC TGCCTTTTCCTTTCCCTAGAGCCTCCACCACCCCGAGATCACATTTCTCACTGCCTTTTGTCTGCCCAGTTTCACCAGAAGTAGGCCTCTTCCTGACAGGCAGCTGCACCACTGCCTGGCGCTGTGCCCTTCCTTTGCTCTGCCCGCTGGAGACGGTGTTTG TCATGGGCCTGGTCTGCAGGGATCCTGCTACAAAGGTGAAACCCAGGAGAGTGTGGAGTCCAGAGTGTTGCCAGGACCCAGGCACAGGCATTAGTGCCCGTTGGAGAAAACAGGGGAATCCCGAAGAAATGGTGGGTCCTGGCCATCCGTGAGATCTTCCCA GGGCAGCTCCCCTCTGTGGAATCCAATCTGTCTTCCATCCTGCGTGGCCGAGGGCCAGGCTTCTCACTGGGCCTCTGCAGGAGGCTGCCATTTGTCCTGCCCACCTTCTTAGAAGCGAGACGGAGCAGACCCATCTGCTACTGCCCTTTCTATAATAACTAA AGTTAGCTGCCCTGGACTATTCACCCCCTAGTCTCAATTTAAGAAGATCCCCATGGCCACAGGGCCCCTGCCTGGGGGCTTGTCACCTCCCCCACCTTCTTCCTGAGTCATTCCTGCAGCCTTGCTCCCTAACCTGCCCCACAGCCTTGCCTGGATTTCTAT CTCCCTGGCTTGGTGCCAGTTCCTCCAAGTCGATGGCACCTCCCTCCCTCTCAACCACTTGAGCAAACTCCAAGACATCTTCTACCCCAACACCAGCAATTGTGCCAAGGGCCATTAGGCTCTCAGCATGACTATTTTTAGAGACCCCGTGTCTGTCACTGA AACCTTTTTTGTGGGAGACTATTCCTCCCATCTGCAACAGCTGCCCCTGCTGACTGCCCTTCTCTCCTCCCTCTCATCCCAGAGAAACAGGTCAGCTGGGAGCTTCTGCCCCCACTGCCTAGGGACCAACAGGGGCAGGAGGCAGTCACTGACCCCGAGACG …
  • Re-naming just part of header in fasta file
    written by seta seta
    Hi all, I usually use some command to rename the header, but they changed the whole header. I want to re-name just part of the …
  • How to remove some fasta sequences by header information from a large fasta file, any command and s…
    written by seta seta
    Dear all, I would like to remove some fasta sequences within a large fasta file in term of header information (sequence name), could anybody please …
  • removing sequencing gaps
    written by cdwilliam524 • cdwilliam524
    Could anyone tell me how to remove sequencing gaps from a *.sam file? I use **bwa** to map the reference genome (fasta) to the metatranscriptome …