gravatar for anasofiamoreira94

2 hours ago by

I want to merge some sequences created from my lab with the nt databases.
So I did a makeblastdb command and then blastdb_aliastool.However the blast output doesn't give me a taxonomic name, it returns NA.
So now, I want to create a costume taxonomy file.
Here is an example of my fasta headers:

>se1 name_species gene

>se2 name_species gene

And then I created the custom taxonomy file like this:

se1 123452
se2 123445

I added the -taxid_map to the makeblastdb command, but results continues to be NA...
Thanks for the help



Source link