This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Remove Duplicate Reads From Fasta File

Hi all,

I want to remove duplicate reads from my fasta file. I tried to use fastx_collapser. But since my reads contains lowercase letters and hyphens it failed.

Please help.

Thanks,D.

fasta read

It's like everybody wants to remove duplicates here!

3 answers

Try the sequniq tool from the GenomeTools suite:

gt sequniq -o output.fasta input.fasta

i tried this command, plz could u tell how this command applied...

Try CD-hit or Uclust

You can remove unwanted hyphens and convert to uppercase using sed:

echo FaSta-TEst | sed "s/-//g ; s/(.*)/\U&/g"

Or just tr: echo FaSta-TEst | tr -d - | tr 'a-z' 'A-Z'

Here is my free program on Github Sequence database curator (https://github.com/Eslam-Samir-Ragab/Sequence-database-curator)

It is a very fast program and it can deal with:

  1. Nucleotide sequences
  2. Protein sequences

It can work under Operating systems:

  1. Windows
  2. Mac
  3. Linux

It also works for:

  1. Fasta format
  2. Fastq format

Best Regards

Log in to answer this question.