Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Remove duplicate sequences with same id from a fasta file
1
gravatar for CB
2.9 years ago by
CB • 10
US
CB • 10 wrote:

Dear all, there are many posts about remove duplicate sequences in a fasta file (https://www.biostars.org/p/3003/), but I want to remove only the duplicate sequences with the same ids.

I have many duplicate sequences in my fasta file, but with different ids and I want to keep them.

How to remove only same id sequence duplicates? I have protein sequences and my sequences are split in different lines.
duplicate sequence remove fasta • 3.5k views
ADD COMMENT • link •
modified 2.9 years ago by Alex Reynolds ♦ 29k • written 2.9 years ago by CB • 10

BBMap's Dedupe utility has a "requirematchingnames" flag. This will make it only remove duplicates that have identical sequence and identical names. For example:

dedupe.sh in=file.fasta out=deduped.fasta ac=f requirematchingnames

One copy of each duplicate set will remain, unless you add the "uniqueonly" flag.

ADD REPLY • link modified 2.9 years ago • written 2.9 years ago by Brian Bushnell ♦ 17k
3
gravatar for shenwei356
2.9 years ago by
shenwei356 ♦ 5.0k
China
shenwei356 ♦ 5.0k wrote:

http://bioinf.shenwei.me/seqkit/usage/#rmdup
ADD COMMENT • link written 2.9 years ago by shenwei356 ♦ 5.0k

It worked very well. It is very easy to use. Thanks!

ADD REPLY • link written 2.9 years ago by CB • 10
2
gravatar for Alex Reynolds
2.9 years ago by
Alex Reynolds ♦ 29k
Seattle, WA USA
Alex Reynolds ♦ 29k wrote:

Make your FASTA files single-line (one line for header, one line for sequence): Multiline Fasta To Single Line Fasta

Then:

$ awk '{ if (($0 ~ /^>/) && (!seen[$0]++)) { print $0; printSeq=1; } else if (($0 ~ /^[^>]/) && printSeq) { print $0; printSeq=0; }  }' in.fa > out.fa
ADD COMMENT • link modified 2.9 years ago • written 2.9 years ago by Alex Reynolds ♦ 29k
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 944 users visited in the last hour