1 answer
I don't know about orthoMCL, but if you just want to change the header and make them unique, do the following (in linux, or install GnuWin32 from here for Windows to get gawk command)
gawk '{if ($0 ~/^>/) {h[$1]++; $1=$1 "_" h[$1]} print}' myfasta.fa > updatedIDs_myfasta.fa
# myfasta.fa is your fasta file.
•
written
by
biocyberman
That's weird, my gawk-fu can't be failing :-) Could you post an excerpt of the fasta file with sequences trimmed to about 10 bases?
>c10047_g1_i1|m.4145 c10047_g1_i1|g.4145 ORF c10047_g1_i1|g.4145 c10047_g1_i1|m.4145 type:complete len:387 (-) c10047_g1_i1:511-1671(-)</p>
>c10047_g2_i1|m.4146 c10047_g2_i1|g.4146 ORF c10047_g2_i1|g.4146 c10047_g2_i1|m.4146 type:5prime_partial len:589 (+) c10047_g2_i1:2-1768(+)
These are headers of my fasta file. The same headers I want to merge or remove for my next work. The headers have different sequence.
Oh, this is different from what you gave in the question. In your fasta file, the tools that generated it form unique headers like this: one header: c10047_g1_i1|m.4145; another header: c10047_g2_i1|m.4146 but orthoMCL propably only consider header before the pipe '|' signs. Therefore you can make this:
gawk 'BEGIN{FS=" "}{if ($0 ~/^>/){gsub("\\|", "pp", $1)} print}' myfasta.fa >updatedIDs_myfasta.fa
Change "pp" to anything you like, but keep it distinguishable.
•
written
by
biocyberman
Log in to answer this question.
-
Header conflicts in concatenated FASTQ file?
written by Dunois •
-
how to add > to fasta file header and merge two line headers into one line(see below)
written by yaqinguo629 •
-
Extracting fasta sequence from CDS matching with sequence ID
written by metallicasster •
-
Merge multiple overlapping consensus sequences
written by kspata
-
Combining two fasta sequences into one
written by Lille My •
-
merge two mulitfasta files and eliminate fasta with duplicated headers from the first
written by Pawel Osipowski •
-
fasta - reverse complement sequence
written by mbk0asis
-
Creating Unique Fasta Headers
written by peerlesshm •
-
finding orthologs in four genomes
written by Mehmet
-
Frequency of a unique sequence in a fasta file
written by Tanvir Ahamed
It seems that your upstream tool spit out different fragments of the same sequence. Merge them with same padding 'N' may work, but the quicker and better method is to make the headers unique.