Difficulties in getting blast to work in assigning taxonomy.

Skip to first unread message


Mar 31, 2016, 8:18:34 AM3/31/16
to Qiime 1 Forum

Hi there,

Sorry to be a bother, but I am having difficulties in getting blast to work when assigning taxonomies. I get the following error message:

MacQIIME zn10508:~ $ assign_taxonomy.py -i /Users/christopherbarnes/Desktop/Postdoc_lab_results/Thapsia_project/ITSx_rep_set.fasta.ITS2.fasta -r /Users/christopherbarnes/Desktop/Sequencing\ resources/its_12_11_otus/rep_set/97_otus.fasta -t /Users/christopherbarnes/Desktop/Sequencing\ resources/its_12_11_otus/taxonomy/97_otu_taxonomy.txt -o assigned_taxonomy_old_unite_blast -m blast

Traceback (most recent call last):

  File "/macqiime/anaconda/bin/assign_taxonomy.py", line 417, in <module>


  File "/macqiime/anaconda/bin/assign_taxonomy.py", line 394, in main


  File "/macqiime/anaconda/lib/python2.7/site-packages/qiime/assign_taxonomy.py", line 500, in __call__

    abspath(reference_seqs_path), output_dir=blast_db_dir)

  File "/macqiime/anaconda/lib/python2.7/site-packages/bfillings/formatdb.py", line 121, in build_blast_db_from_fasta_path

    app_result = fdb(fasta_path)

  File "/macqiime/anaconda/lib/python2.7/site-packages/burrito/util.py", line 285, in __call__

    'StdErr:\n%s\n' % open(errfile).read())

burrito.util.ApplicationError: Unacceptable application exit status: 3


cd "/tmp/bltax-kTF6Ry/"; formatdb -l "97_otus.fasta.log" -o T -n "97_otus.fasta" -i "/Users/christopherbarnes/Desktop/Sequencing resources/its_12_11_otus/rep_set/97_otus.fasta" -p F > "/dev/null" 2> "/dev/null"



I have blast installed so that the output seems fine:

which blastall


which formatdb


And finally the output from print_qiime_config.py is:

System information


         Platform: darwin

   Python version: 2.7.10 |Anaconda 2.2.0 (x86_64)| (default, May 28 2015, 17:04:42)  [GCC 4.2.1 (Apple Inc. build 5577)]

Python executable: /macqiime/anaconda/bin/python

QIIME default reference information


For details on what files are used as QIIME's default references, see here:


Dependency versions


          QIIME library version: 1.9.1

           QIIME script version: 1.9.1

qiime-default-reference version: 0.1.2

                  NumPy version: 1.9.2

                  SciPy version: 0.15.1

                 pandas version: 0.16.1

             matplotlib version: 1.4.3

            biom-format version: 2.1.4

                   h5py version: 2.4.0 (HDF5 version: 1.8.14)

                   qcli version: 0.1.1

                   pyqi version: 0.3.2

             scikit-bio version: 0.2.3

                 PyNAST version: 1.2.2

                Emperor version: 0.9.51

                burrito version: 0.9.1

       burrito-fillings version: 0.1.1

              sortmerna version: SortMeRNA version 2.0, 29/11/2014

              sumaclust version: SUMACLUST Version 1.0.00

                  swarm version: Swarm 1.2.19 [Jun  2 2015 14:40:16]

                          gdata: Installed.

QIIME config values


For definitions of these settings and to learn how to configure QIIME, see here:



                     blastmat_dir: None

      pick_otus_reference_seqs_fp: /macqiime/anaconda/lib/python2.7/site-packages/qiime_default_reference/gg_13_8_otus/rep_set/97_otus.fasta

                         sc_queue: all.q

      topiaryexplorer_project_dir: None

     pynast_template_alignment_fp: /macqiime/anaconda/lib/python2.7/site-packages/qiime_default_reference/gg_13_8_otus/rep_set_aligned/85_otus.pynast.fasta

                  cluster_jobs_fp: start_parallel_jobs.py

pynast_template_alignment_blastdb: None

assign_taxonomy_reference_seqs_fp: /macqiime/anaconda/lib/python2.7/site-packages/qiime_default_reference/gg_13_8_otus/rep_set/97_otus.fasta

                     torque_queue: friendlyq

                    jobs_to_start: 1

                       slurm_time: None

            denoiser_min_per_core: 50

assign_taxonomy_id_to_taxonomy_fp: /macqiime/anaconda/lib/python2.7/site-packages/qiime_default_reference/gg_13_8_otus/taxonomy/97_otu_taxonomy.txt

                         temp_dir: /tmp/

                     slurm_memory: None

                      slurm_queue: None

                      blastall_fp: blastall

                 seconds_to_sleep: 60

Any help would be appreciated as I'm completely stumped!

Thanks in advance for any help you can offer.


Jenya Kopylov

Mar 31, 2016, 9:44:12 AM3/31/16
to Qiime 1 Forum
Hi Chris,

Could you check whether Blast taxonomy assignment tests are passing on your system?


$ python /macqiime/QIIME/tests/test_assign_taxonomy.py BlastTaxonAssignerTests
Ran 11 tests in 2.653s

Also, it seems you're using a 16S database to assign taxonomy to ITS sequences, am I interpreting your command line correctly?



Mar 31, 2016, 10:41:55 AM3/31/16
to Qiime 1 Forum
Hi Jenya,

Thanks for your response. 

Ah, so my command is definitely using fungal databases (UNITE databases), however the defaults from my print_qiime_configs.py are the greengenes database. Does this have an effect if I always specify databases I wish to assign to? I never noticed that before.

Regarding blast output:

MacQIIME zn10508:~ $ python /macqiime/QIIME/tests/test_assign_taxonomy.py BlastTaxonAssignerTests



Ran 11 tests in 0.634s


All tests of BLAST seem ok, could it be to do with the databases? Either the UNITE databases or something to do with these defaults within QIIIME? 

If it's of relevance, this database works fine using the default UCLUST taxonomic assignments.

Many thanks for your time,


Jenya Kopylov

Mar 31, 2016, 11:09:36 AM3/31/16
to Qiime 1 Forum
Hi Chris,

It's good to see the tests passing. One possible issue could be spaces between your folder name (ex. Sequencing resources), could you replace all spaces (ex. with '_') and retry your assign_taxonomy command?

From your command you are indeed using ITS so that shouldn't be an issue.



Mar 31, 2016, 12:30:05 PM3/31/16
to Qiime 1 Forum
Hi Jenya,

It worked, wow I feel a bit silly, totally didn't think about that. I know it's good form not to leave spaces, but I guess this really drums in the message!

Thanks so much for your time, hugely appreciated!

Reply all
Reply to author
0 new messages