This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error: Duplicate seq_ids are found

Hello everyone, I am using the web-based Galaxy tool, not the command line version. I merged FASTA files into one and I'm trying to construct a BLAST database with these local sequences through the makeblastdb function. I get an error message that reads, "Error: Duplicate seq_ids are found: GNL|BL_ORD_ID|9650923".

Can anyone assist in finding a way to remove the duplicate seq IDs using the web-based Galaxy tool preferentially?

Thank you!

galaxy blast
• 2,026 views • link
updated by cpad0112 cpad0112 • written by rupjandu • rupjandu
  1. Check if there are any fasta ID / header repairing tools in galaxy toolshed and see if they are not installed on galaxy instance.
  2. If they are not, try to install them if you are admin or if you have access to admin account. If not request admin to install such a tool.
  3. If 1 and 2 are not possible, download fasta file.
  4. Run this function sed -nr '/^>/p' <input.fa> |sort -V | uniq -D | uniq -c on download file (input.fa). This should print duplicated/identical headers and their count.
  5. Download seqkit tool and run seqkit rename -n <input.fa> -o <output.fa>. This would generate a new file output.fa and append numbers serially at the end of fasta IDs/headers if they are identical.
  6. Run this function sed -nr '/^>/p' <output.fa> | sort -V | uniq -D | uniq -c on new file (output.fa). This should not print any line.
  7. Upload the new fasta file and run blast on it.
• 0 views • link
written by cpad0112 cpad0112

0 answers

No answers yet.

Log in to answer this question.

More posts like this