Select human transcripts on the Diamond aligner


I’m new to RNA seq and I’m trying to use blastx from the Diamond aligner on galaxy to blast a fasta file of possible novel lncRNAs using the ncbi nr database against human transcripts. However, I am not able to select the human taxon only. I have tried entering ‘homo sapiens’ and ‘9606’ in the option to restrict the search taxonomically, but both do not work and give me an error: ‘Error: Output format requires taxonomy mapping information built into the database (use --taxonmap parameter for the makedb command).’ Does anyone know how to do this?

Hi @bart_joosten,
this option hasn’t been implemented correctly. I’m going to work on it (PR); I’ll try to make it available for the next Monday.