Hi everyone,
I would like to build a human+bacteria+archaea subset of the nr database. So what I did is:
- Get the nr database
- Get the GI list of human+bacteria+archaea from Entrez Protein database
Use this GI list with blastdb_aliastool to build an aliased database of the subset proteins. I used the following command and get the error
blastdb_aliastool -db nr -gilist sequence_human_bacteria_archaea_26112013.gi -dbtype prot -out human_bacteria_archaea_nov2013 -title "Human plus bacteria and archaea entries"
Converted 74595073 GIs from sequence_human_bacteria_archaea_26112013.gi to binary format in human_bacteria_archaea_nov2013.p.gil BLAST Database error: BLASTDB alias file creation failed. Some referenced files may be missing
Does anyone have any idea what the problem is?
Thanks!