diff qiime2/merge_feature_table_and_taxonomy.xml @ 0:370e0b6e9826 draft

Uploaded
author florianbegusch
date Wed, 17 Jul 2019 03:05:17 -0400
parents
children 914fa4daf16a
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/qiime2/merge_feature_table_and_taxonomy.xml	Wed Jul 17 03:05:17 2019 -0400
@@ -0,0 +1,34 @@
+<?xml version="1.0" ?>
+<tool id="merge_feature_table_and_taxonomy" name="merge feature table and taxonomy" version="2019.4">
+	<description>- Merge features and taxonomy into a single biom file.</description>
+	<requirements>
+		<requirement type="package" version="2019.4">qiime2</requirement>
+	</requirements>
+	<command>
+mkdir out;
+
+qiime tools export --input-path --output-path out $taxonomy_qza;
+qiime tools export --input-path --output-path out $biom_qza;
+
+sed -i '1s;^;#;' out/*.tsv;
+sed -i 's/Confidence/confidence/g' out/*.tsv;
+sed -i 's/Taxon/taxonomy/g' out/*.tsv;
+sed -i 's/Feature ID/OTU ID/g' out/*.tsv;
+
+biom add-metadata -i out/*.biom -o table_with_sample_metadata.biom --observation-metadata-fp out/*.tsv
+;
+mv *.biom $biom_and_tax_merged</command>
+	<inputs>
+		<param label="Feature table" name="biom_qza" type="data" format="qza,no_unzip.zip" collection_type="list"/>
+		<param label="Taxonomy" name="taxonomy_qza" type="data" format="qza,no_unzip.zip"/>
+	</inputs>
+	<outputs>
+		<data format="biom2" label="${tool.name} on ${on_string}: biom_and_tax_merged.biom" name="biom_and_tax_merged"/>
+	</outputs>
+	<help>
+	</help>
+<macros>
+	<import>qiime_citation.xml</import>
+</macros>
+<expand macro="qiime_citation" />
+</tool>