Mercurial > repos > iuc > humann2_join_tables
diff generate_test_data @ 0:14d5bf0d63a0 draft
planemo upload for repository https://github.com/ASaiM/galaxytools/tree/master/tools/humann2/ commit b46aa969c01b7e5f4f133192899fa4da286ecf89-dirty
author | iuc |
---|---|
date | Mon, 13 Mar 2017 12:36:38 -0400 |
parents | |
children | 9ada8e95bdd0 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/generate_test_data Mon Mar 13 12:36:38 2017 -0400 @@ -0,0 +1,124 @@ +#/usr/bin/env bash + +# humann2 +humann2 \ + --input 'test-data/input_sequences.fasta' \ + -o 'humann2_output' \ + --annotation-gene-index 8 \ + --taxonomic-profile 'test-data/taxonomic_profile.tabular' \ + --metaphlan-options="-t rel_ab" \ + --nucleotide-database 'test-data/' \ + --protein-database 'test-data/' \ + --evalue '1' \ + --search-mode 'uniref50' \ + --prescreen-threshold '0.01' \ + --identity-threshold '50' \ + --translated-subject-coverage-threshold '50' \ + --translated-query-coverage-threshold '50' \ + --translated-alignment 'diamond' \ + --xipe 'off' \ + --minpath 'on' \ + --pick-frames 'on' \ + --gap-fill 'off' \ + --output-format 'tsv' \ + --output-max-decimals '10' \ + --output-basename 'humann2' \ + --pathways 'metacyc' + +humann2 \ + --input 'test-data/input_sequences.fasta' \ + -o 'humann2_output' \ + --annotation-gene-index 8 \ + --metaphlan-options="-t rel_ab" \ + --nucleotide-database 'test-data/' \ + --protein-database 'test-data/' \ + --evalue '1' \ + --search-mode 'uniref90' \ + --prescreen-threshold '0.01' \ + --identity-threshold '50' \ + --translated-subject-coverage-threshold '50' \ + --translated-query-coverage-threshold '50' \ + --translated-alignment 'rapsearch' \ + --xipe 'off' \ + --minpath 'on' \ + --pick-frames 'on' \ + --gap-fill 'off' \ + --output-format 'tsv' \ + --output-max-decimals '10' \ + --output-basename 'humann2' \ + --pathways 'unipathway' + +# humann2_regroup_table +humann2_regroup_table \ + --input 'test-data/demo_genefamilies.tsv' \ + -o 'test-data/regrouped_gene_families_to_infogo1000.tsv' \ + --function 'sum' \ + --groups 'uniref90_infogo1000' \ + --ungrouped 'Y' \ + --protected 'Y' + +humann2_regroup_table \ + --input 'test-data/demo_genefamilies.tsv' \ + -o 'test-data/regrouped_gene_families_to_ko.tsv' \ + --function 'mean' \ + --groups 'uniref90_ko' \ + --ungrouped 'Y' \ + --protected 'Y' + +humann2_regroup_table \ + --input 'test-data/demo_genefamilies.tsv' \ + -o 'test-data/regrouped_gene_families_to_rxn.tsv' \ + --function 'sum' \ + --groups 'uniref90_rxn' \ + --ungrouped 'Y' \ + --protected 'Y' + +# humann2_renorm_table +humann2_renorm_table \ + --input 'test-data/demo_pathabundance.tsv' \ + -o 'test-data/cpm_community_renormalized_pathway_abundance.tsv' \ + --units 'cpm' \ + --mode 'community' \ + --special 'n' + +humann2_renorm_table \ + --input 'test-data/demo_pathabundance.tsv' \ + -o 'test-data/relab_levelwise_renormalized_pathway_abundance.tsv' \ + --units 'relab' \ + --mode 'levelwise' \ + --special 'y' + +# humann2_join_tables +mkdir join_table_tmp_dir +cp 'test-data/demo_pathabundance.tsv' join_table_tmp_dir +cp 'test-data/demo_pathcoverage.tsv' join_table_tmp_dir +humann2_join_tables \ + -i 'join_table_tmp_dir' \ + -o 'test-data/joined_pathway_coverage_abundance.tsv' + +# humann2_merge_abundance_tables +humann2_merge_abundance_tables \ + --input-genes 'test-data/demo_genefamilies.tsv' \ + --input-pathways 'test-data/demo_pathabundance.tsv' \ + -o 'test-data/merged_gene_families_pathways_abundances.tsv' + +# humann2_reduce_table +humann2_reduce_table \ + --input 'test-data/demo_genefamilies.tsv' \ + -o 'test-data/max_reduced_gene_family_abundance.tsv' \ + --function 'max' \ + --sort-by 'name' + +# humann2_rename_table +humann2_rename_table \ + --input 'test-data/demo_genefamilies.tsv' \ + -o 'test-data/renamed_genefamilies.tsv' \ + --names 'metacyc-pwy' + +# humann2_split_table +mkdir split_table_tmp_dir +humann2_split_table \ + --input 'test-data/joined_pathway_coverage_abundance.tsv' \ + -o 'split_table_tmp_dir' +cp 'split_table_tmp_dir/demo_Abundance.tsv' 'test-data/' +cp 'split_table_tmp_dir/demo_Coverage.tsv' 'test-data/'