# PICRUSt command
**URL:** https://forum.qiime2.org/t/picrust-command/3985
**Category:** User Support
**Created:** [April 30, 2018, 11:35pm UTC](https://forum.qiime2.org/t/picrust-command/3985 "2018-04-30T23:35:32Z")
**Posts on this page:** 7
**Page:** 1
### Author:  [@Harry](https://forum.qiime2.org/u/Harry)
#### Post date: [April 30, 2018, 11:35pm UTC](https://forum.qiime2.org/t/picrust-command/3985/1 "2018-04-30T23:35:32Z")
Hi,
May I know after this ([Picrust support?](https://forum.qiime2.org/t/picrust-support/1376)), what command should be used to get the PICRUSt format biom file?
I use the command from this discussion ([How to create a feature table with qiime2 for PICRUST with the taxonomic assignment? - #12 by jairideout](https://forum.qiime2.org/t/how-to-create-a-feature-table-with-qiime2-for-picrust-with-the-taxonomic-assignment/2526/12)), specifically: biom add-metadata -i ExportedTableForPICRUSt/feature-table.biom -o table-with-taxonomy.biom --observation-metadata-fp 97\_otu\_taxonomy.txt --sc-separated taxonomy
But it does not work, and gives me the errors:
Traceback (most recent call last):
File "/Users/ranranhu/miniconda3/envs/qiime2-2018.4/bin/biom", line 6, in
sys.exit(biom.cli.cli())
File "/Users/ranranhu/miniconda3/envs/qiime2-2018.4/lib/python3.5/site-packages/click/core.py", line 722, in **call**
return self.main(\*args, \*\*kwargs)
File "/Users/ranranhu/miniconda3/envs/qiime2-2018.4/lib/python3.5/site-packages/click/core.py", line 697, in main
rv = self.invoke(ctx)
File "/Users/ranranhu/miniconda3/envs/qiime2-2018.4/lib/python3.5/site-packages/click/core.py", line 1066, in invoke
return \_process\_result(sub\_ctx.command.invoke(sub\_ctx))
File "/Users/ranranhu/miniconda3/envs/qiime2-2018.4/lib/python3.5/site-packages/click/core.py", line 895, in invoke
return ctx.invoke(self.callback, \*\*ctx.params)
File "/Users/ranranhu/miniconda3/envs/qiime2-2018.4/lib/python3.5/site-packages/click/core.py", line 535, in invoke
return callback(\*args, \*\*kwargs)
File "/Users/ranranhu/miniconda3/envs/qiime2-2018.4/lib/python3.5/site-packages/biom\_format-2.1.6-py3.5-macosx-10.7-x86\_64.egg/biom/cli/metadata\_adder.py", line 107, in add\_metadata
float\_fields, sample\_header, observation\_header)
File "/Users/ranranhu/miniconda3/envs/qiime2-2018.4/lib/python3.5/site-packages/biom\_format-2.1.6-py3.5-macosx-10.7-x86\_64.egg/biom/cli/metadata\_adder.py", line 174, in \_add\_metadata
header=observation\_header)
File "/Users/ranranhu/miniconda3/envs/qiime2-2018.4/lib/python3.5/site-packages/biom\_format-2.1.6-py3.5-macosx-10.7-x86\_64.egg/biom/parse.py", line 527, in from\_file
raise BiomParseException("No header line was found in mapping "
biom.exception.BiomParseException: No header line was found in mapping file
By the way, I used dada2 to denoise.
Thank you so much!
---
### Author:  [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [May 1, 2018, 12:11pm UTC](https://forum.qiime2.org/t/picrust-command/3985/2 "2018-05-01T12:11:13Z")
---
### Author:  [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [May 1, 2018, 12:30pm UTC](https://forum.qiime2.org/t/picrust-command/3985/3 "2018-05-01T12:30:29Z")
> [@Harry](#):
>
> By the way, I used dada2 to denoise.
It is fine to use dada2 for initial denoising — but if performing a picrust analysis you must then perform closed-reference OTU picking agains the picrust database using vsearch, as described in the thread that you linked to. Since you mention dada2, I just want to make sure that you saw the vsearch step in that linked protocol.
> [@Harry](#):
>
> biom.exception.BiomParseException: No header line was found in mapping file
Check out this thread, it looks like you are having the same error:
> [@Adding qiime2 taxonomy data to a .biom file](https://forum.qiime2.org/t/adding-qiime2-taxonomy-data-to-a-biom-file/1516):
>
> Hello! I've been trying these steps in order to import the biom file to Qiime1 and summarize the taxa barplot by categories. The export parts worked fine, I think, but when trying to add the metadata values to the biom file this error is returned: \_ biom.exception.BiomParseException: No header line was found in mapping file. \_ I have changed the biom-taxonomy.tsv file header line in order to reproduce what you wrote on your answer, like this: #OTUID taxonomy confidence This way, when …
The error seems fairly clear: there is a formatting issue with your metadata file (i.e., 97\_otu\_taxonomy.txt) that needs to be fixed.
Note that biom-format and picrust are external software and we cannot really provide support for them on this forum, though we will do what we can to help.
Let me know if fixing the header line in your metadata file helps!
---
### Author:  [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [May 1, 2018, 12:30pm UTC](https://forum.qiime2.org/t/picrust-command/3985/4 "2018-05-01T12:30:31Z")
---
### Author:  [@Harry](https://forum.qiime2.org/u/Harry)
#### Post date: [May 2, 2018, 11:49am UTC](https://forum.qiime2.org/t/picrust-command/3985/5 "2018-05-02T11:49:47Z")
Thank you so much Nicholas, it works.
BUT, PICRUSt galaxy server ([http://galaxy.morganlangille.com/](http://galaxy.morganlangille.com/)) seems down, I have not got a chance to test the biom file. I will update here as soon as the server resumes working.
---
### Author:  [@Harry](https://forum.qiime2.org/u/Harry)
#### Post date: [May 2, 2018, 4:51pm UTC](https://forum.qiime2.org/t/picrust-command/3985/6 "2018-05-02T16:51:14Z")
Hi @Nicholas_Bokulich,
The galaxy resumes working, and I tested the biom file. IT WORKS!
The taxonomy file needs to be modified
The 97\_otu\_taxonomy.txt file needs to add "#OTUID taxonomy" at the top, but should be Tab instead of space.
---
### Author:  [@system](https://forum.qiime2.org/u/system)
#### Post date: [June 2, 2018, 10:51pm UTC](https://forum.qiime2.org/t/picrust-command/3985/7 "2018-06-02T22:51:17Z")
This topic was automatically closed 31 days after the last reply. New replies are no longer allowed.