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