# Problems with q2-quality-control using Mockrobiota

**URL:** https://forum.qiime2.org/t/problems-with-q2-quality-control-using-mockrobiota/10374
**Category:** Technical Support
**Tags:** mock-communities, quality-control
**Created:** [June 22, 2019, 2:27pm UTC](https://forum.qiime2.org/t/problems-with-q2-quality-control-using-mockrobiota/10374 "2019-06-22T14:27:00Z")
**Posts on this page:** 1
**Showing post:** 3

<div class="post-metadata">

### Author: ![xchromosome](https://forum.qiime2.org/user_avatar/forum.qiime2.org/xchromosome/32/13536_2.png) [@xchromosome](https://forum.qiime2.org/u/xchromosome)
#### Post date: [June 23, 2019, 12:03pm UTC](https://forum.qiime2.org/t/problems-with-q2-quality-control-using-mockrobiota/10374/3 "2019-06-23T12:03:38Z")

</div>

Hiya,

Thanks for the quick response, especially on a Saturday!

> [@Nicholas\_Bokulich](#):
>
> I think you took a wrong turn pretty early on. It sounds like you were taking the fastq sequence data from mockrobiota, processing it, and then attempting to use that as the “expected composition”

Yup, that's exactly what I was doing 😑 I was trying to follow the instructions I'd found from previous threads but I must have got mixed up!

> [@Nicholas\_Bokulich](#):
>
> Instead, you need to find the “expected taxonomy” file for that mock community and import that as a relative frequency table.
> 
> See this tutorial; it is using a fungal mock community but gives you a good idea of how the “expected taxonomy” from mockrobiota should be used:

Ok, I have got to the page where there are 2 things I can download: expected-sequences.fasta or taxonomy.tsv. I right-clicked on taxonomy.tsv and chose "save link as...".  
According to the fungal ITS analysis tutorial, the next step is as follows:

> [@Fungal ITS analysis tutorial](https://forum.qiime2.org/t/fungal-its-analysis-tutorial/7351/1):
>
> That file was annotated with an older version of the UNITE database, so let’s fix some of the annotations (to match the version we are using) before proceeding:
> 
> ```auto
> sed 's/s __Tylospora_asterophora/s__ unidentified/' expected-taxonomy.tsv | sed 's/p __Zygomycota;c__ Mortierellomycotina_cls_Incertae_sedis/p __Mortierellomycota;c__ Mortierellomycetes/' | sed 's/c __Chytridiomycetes/c__ Spizellomycetes/' | sed 's/g __Coprinopsis/g__ Coprinopsis;s __unidentified/' | sed 's/g__ Tricholoma/g __Tricholoma;s__ Tricholoma_vaccinum/'> expected-taxonomy-mod.tsv
> 
> ```

I don't really know what this means - is it relevant to me? I thought it might be something that only applied to older versions of Qiime 2 so I skipped it to see if I could. The next thing I ran was this:

```auto
biom convert \
  -i mockrobiota-taxonomy.tsv \
  -o mockrobiota-taxonomy.biom \
  --table-type="OTU table" \
  --to-json

```

This was the error message I got:

Traceback (most recent call last):  
File "/home/qiime2/miniconda/envs/qiime2-2019.1/lib/python3.6/site-packages/biom/parse.py", line 660, in load\_table  
table = parse\_biom\_table(fp)  
File "/home/qiime2/miniconda/envs/qiime2-2019.1/lib/python3.6/site-packages/biom/parse.py", line 412, in parse\_biom\_table  
t = Table.from\_tsv(fp, None, None, lambda x: x)  
File "/home/qiime2/miniconda/envs/qiime2-2019.1/lib/python3.6/site-packages/biom/table.py", line 4631, in from\_tsv  
t\_md\_name) = Table.\_extract\_data\_from\_tsv(lines, \*\*kwargs)  
File "/home/qiime2/miniconda/envs/qiime2-2019.1/lib/python3.6/site-packages/biom/table.py", line 4747, in \_extract\_data\_from\_tsv  
md\_name = header[-1]  
IndexError: list index out of range

During handling of the above exception, another exception occurred:

Traceback (most recent call last):  
File "/home/qiime2/miniconda/envs/qiime2-2019.1/bin/biom", line 11, in   
sys.exit(cli())  
File "/home/qiime2/miniconda/envs/qiime2-2019.1/lib/python3.6/site-packages/click/core.py", line 764, in **call**  
return self.main(\*args, \*\*kwargs)  
File "/home/qiime2/miniconda/envs/qiime2-2019.1/lib/python3.6/site-packages/click/core.py", line 717, in main  
rv = self.invoke(ctx)  
File "/home/qiime2/miniconda/envs/qiime2-2019.1/lib/python3.6/site-packages/click/core.py", line 1137, in invoke  
return \_process\_result(sub\_ctx.command.invoke(sub\_ctx))  
File "/home/qiime2/miniconda/envs/qiime2-2019.1/lib/python3.6/site-packages/click/core.py", line 956, in invoke  
return ctx.invoke(self.callback, \*\*ctx.params)  
File "/home/qiime2/miniconda/envs/qiime2-2019.1/lib/python3.6/site-packages/click/core.py", line 555, in invoke  
return callback(\*args, \*\*kwargs)  
File "/home/qiime2/miniconda/envs/qiime2-2019.1/lib/python3.6/site-packages/biom/cli/table\_converter.py", line 114, in convert  
table = load\_table(input\_fp)  
File "/home/qiime2/miniconda/envs/qiime2-2019.1/lib/python3.6/site-packages/biom/parse.py", line 662, in load\_table  
raise TypeError("%s does not appear to be a BIOM file!" % f)  
TypeError: mockrobiota-taxonomy.tsv does not appear to be a BIOM file!

Where have I screwed up now?! 😨

> [@Nicholas\_Bokulich](#):
>
> yeah, q2-quality-control will show you actual accuracy scores (excel will not), but it will not make a bar chart of the compositions (excel will).

It sounds like it would be beneficial to do both! That's my plan!

Massive thanks,  
Lindsay

---

_[View the full topic](https://forum.qiime2.org/t/problems-with-q2-quality-control-using-mockrobiota/10374)._