Logging started at 11:14:24 on 03 Aug 2017 QIIME version: 1.9.1 qiime_config values: pick_otus_reference_seqs_fp /usr/lib/python2.7/site-packages/qiime_default_reference/gg_13_8_otus/rep_set/97_otus.fasta sc_queue all.q pynast_template_alignment_fp /usr/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 assign_taxonomy_reference_seqs_fp /usr/lib/python2.7/site-packages/qiime_default_reference/gg_13_8_otus/rep_set/97_otus.fasta torque_queue friendlyq jobs_to_start 1 denoiser_min_per_core 50 assign_taxonomy_id_to_taxonomy_fp /usr/lib/python2.7/site-packages/qiime_default_reference/gg_13_8_otus/taxonomy/97_otu_taxonomy.txt temp_dir /tmp/ blastall_fp blastall seconds_to_sleep 1 parameter file values: parallel:jobs_to_start 1 Input file md5 sums: otu_table_100.biom: 63b3be1644c48db79eff79d3e1a39b9f Diet_mapping.tsv: 3d4b6631f82380180efa3849de626af2 phylogeny.tre: 83992688af38a4a84a3785f986258c4e Executing commands. # Generate BIOM table summary command biom summarize-table -i otu_table_100.biom -o cdout_noless100/biom_table_summary.txt Stdout: Stderr: # Filter low sequence count samples from table (minimum sequence count: 2000) command filter_samples_from_otu_table.py -i otu_table_100.biom -o cdout_noless100/table_mc2000.biom -n 2000 Stdout: Stderr: # Rarify the OTU table to 2000 sequences/sample command single_rarefaction.py -i cdout_noless100/table_mc2000.biom -o cdout_noless100/table_even2000.biom -d 2000 Stdout: Stderr: Executing commands. # Beta Diversity (weighted_unifrac) command beta_diversity.py -i cdout_noless100/table_even2000.biom -o cdout_noless100/bdiv_even2000/ --metrics weighted_unifrac -t phylogeny.tre Stdout: Stderr: # Rename distance matrix (weighted_unifrac) command mv cdout_noless100/bdiv_even2000//weighted_unifrac_table_even2000.txt cdout_noless100/bdiv_even2000//weighted_unifrac_dm.txt Stdout: Stderr: # Principal coordinates (weighted_unifrac) command principal_coordinates.py -i cdout_noless100/bdiv_even2000//weighted_unifrac_dm.txt -o cdout_noless100/bdiv_even2000//weighted_unifrac_pc.txt Stdout: Stderr: /usr/lib64/python2.7/site-packages/skbio/stats/ordination/_principal_coordinate_analysis.py:107: RuntimeWarning: The result contains negative eigenvalues. Please compare their magnitude with the magnitude of some of the largest positive eigenvalues. If the negative ones are smaller, it's probably safe to ignore them, but if they are large in magnitude, the results won't be useful. See the Notes section for more details. The smallest eigenvalue is -0.00550270387315 and the largest is 0.176110397868. RuntimeWarning # Make emperor plots, weighted_unifrac) command make_emperor.py -i cdout_noless100/bdiv_even2000//weighted_unifrac_pc.txt -o cdout_noless100/bdiv_even2000//weighted_unifrac_emperor_pcoa_plot/ -m Diet_mapping.tsv Stdout: Stderr: # Beta Diversity (unweighted_unifrac) command beta_diversity.py -i cdout_noless100/table_even2000.biom -o cdout_noless100/bdiv_even2000/ --metrics unweighted_unifrac -t phylogeny.tre Stdout: Stderr: # Rename distance matrix (unweighted_unifrac) command mv cdout_noless100/bdiv_even2000//unweighted_unifrac_table_even2000.txt cdout_noless100/bdiv_even2000//unweighted_unifrac_dm.txt Stdout: Stderr: # Principal coordinates (unweighted_unifrac) command principal_coordinates.py -i cdout_noless100/bdiv_even2000//unweighted_unifrac_dm.txt -o cdout_noless100/bdiv_even2000//unweighted_unifrac_pc.txt Stdout: Stderr: # Make emperor plots, unweighted_unifrac) command make_emperor.py -i cdout_noless100/bdiv_even2000//unweighted_unifrac_pc.txt -o cdout_noless100/bdiv_even2000//unweighted_unifrac_emperor_pcoa_plot/ -m Diet_mapping.tsv Stdout: Stderr: Executing commands. # Alpha rarefaction command multiple_rarefactions.py -i cdout_noless100/table_mc2000.biom -m 10 -x 2000 -s 199 -o cdout_noless100/arare_max2000//rarefaction/ Stdout: Stderr: # Alpha diversity on rarefied OTU tables command alpha_diversity.py -i cdout_noless100/arare_max2000//rarefaction/ -o cdout_noless100/arare_max2000//alpha_div/ -t phylogeny.tre Stdout: Stderr: # Collate alpha command collate_alpha.py -i cdout_noless100/arare_max2000//alpha_div/ -o cdout_noless100/arare_max2000//alpha_div_collated/ Stdout: Stderr: # Removing intermediate files command rm -r cdout_noless100/arare_max2000//rarefaction/ cdout_noless100/arare_max2000//alpha_div/ Stdout: Stderr: # Rarefaction plot: All metrics command make_rarefaction_plots.py -i cdout_noless100/arare_max2000//alpha_div_collated/ -m Diet_mapping.tsv -o cdout_noless100/arare_max2000//alpha_rarefaction_plots/ Stdout: Stderr: /usr/lib/python2.7/site-packages/qiime/make_rarefaction_plots.py:50: MatplotlibDeprecationWarning: The set_axis_bgcolor function was deprecated in version 2.0. Use set_facecolor instead. ax.set_axis_bgcolor(background_color) Executing commands. # Sort OTU Table command sort_otu_table.py -i cdout_noless100/table_mc2000.biom -o cdout_noless100/taxa_plots/table_mc2000_sorted.biom Stdout: Stderr: # Summarize Taxonomy command summarize_taxa.py -i cdout_noless100/taxa_plots/table_mc2000_sorted.biom -o cdout_noless100/taxa_plots/ Stdout: Stderr: # Plot Taxonomy Summary command plot_taxa_summary.py -i cdout_noless100/taxa_plots/table_mc2000_sorted_L2.txt,cdout_noless100/taxa_plots/table_mc2000_sorted_L3.txt,cdout_noless100/taxa_plots/table_mc2000_sorted_L4.txt,cdout_noless100/taxa_plots/table_mc2000_sorted_L5.txt,cdout_noless100/taxa_plots/table_mc2000_sorted_L6.txt -o cdout_noless100/taxa_plots//taxa_summary_plots/ Stdout: Stderr: /usr/lib64/python2.7/site-packages/matplotlib/cbook.py:136: MatplotlibDeprecationWarning: The axisbg attribute was deprecated in version 2.0. Use facecolor instead. warnings.warn(message, mplDeprecation, stacklevel=1) Executing commands. # Collapse samples in OTU table by categories command collapse_samples.py -m Diet_mapping.tsv -b cdout_noless100/table_mc2000.biom --output_biom_fp cdout_noless100/taxa_plots_Diet/Diet_otu_table.biom --output_mapping_fp cdout_noless100/taxa_plots_Diet/Diet_map.txt --collapse_fields 'Diet' Stdout: Stderr: # Sort OTU Table command sort_otu_table.py -i cdout_noless100/taxa_plots_Diet/Diet_otu_table.biom -o cdout_noless100/taxa_plots_Diet/Diet_otu_table_sorted.biom Stdout: Stderr: # Summarize Taxonomy command summarize_taxa.py -i cdout_noless100/taxa_plots_Diet/Diet_otu_table_sorted.biom -o cdout_noless100/taxa_plots_Diet/ Stdout: Stderr: # Plot Taxonomy Summary command plot_taxa_summary.py -i cdout_noless100/taxa_plots_Diet/Diet_otu_table_sorted_L2.txt,cdout_noless100/taxa_plots_Diet/Diet_otu_table_sorted_L3.txt,cdout_noless100/taxa_plots_Diet/Diet_otu_table_sorted_L4.txt,cdout_noless100/taxa_plots_Diet/Diet_otu_table_sorted_L5.txt,cdout_noless100/taxa_plots_Diet/Diet_otu_table_sorted_L6.txt -o cdout_noless100/taxa_plots_Diet//taxa_summary_plots/ Stdout: Stderr: /usr/lib64/python2.7/site-packages/matplotlib/cbook.py:136: MatplotlibDeprecationWarning: The axisbg attribute was deprecated in version 2.0. Use facecolor instead. warnings.warn(message, mplDeprecation, stacklevel=1) Executing commands. # Collapse samples in OTU table by categories command collapse_samples.py -m Diet_mapping.tsv -b cdout_noless100/table_mc2000.biom --output_biom_fp cdout_noless100/taxa_plots_DayofDiet/DayofDiet_otu_table.biom --output_mapping_fp cdout_noless100/taxa_plots_DayofDiet/DayofDiet_map.txt --collapse_fields 'DayofDiet' Stdout: Stderr: # Sort OTU Table command sort_otu_table.py -i cdout_noless100/taxa_plots_DayofDiet/DayofDiet_otu_table.biom -o cdout_noless100/taxa_plots_DayofDiet/DayofDiet_otu_table_sorted.biom Stdout: Stderr: # Summarize Taxonomy command summarize_taxa.py -i cdout_noless100/taxa_plots_DayofDiet/DayofDiet_otu_table_sorted.biom -o cdout_noless100/taxa_plots_DayofDiet/ Stdout: Stderr: # Plot Taxonomy Summary command plot_taxa_summary.py -i cdout_noless100/taxa_plots_DayofDiet/DayofDiet_otu_table_sorted_L2.txt,cdout_noless100/taxa_plots_DayofDiet/DayofDiet_otu_table_sorted_L3.txt,cdout_noless100/taxa_plots_DayofDiet/DayofDiet_otu_table_sorted_L4.txt,cdout_noless100/taxa_plots_DayofDiet/DayofDiet_otu_table_sorted_L5.txt,cdout_noless100/taxa_plots_DayofDiet/DayofDiet_otu_table_sorted_L6.txt -o cdout_noless100/taxa_plots_DayofDiet//taxa_summary_plots/ Stdout: Stderr: /usr/lib64/python2.7/site-packages/matplotlib/cbook.py:136: MatplotlibDeprecationWarning: The axisbg attribute was deprecated in version 2.0. Use facecolor instead. warnings.warn(message, mplDeprecation, stacklevel=1) Executing commands. # Collapse samples in OTU table by categories command collapse_samples.py -m Diet_mapping.tsv -b cdout_noless100/table_mc2000.biom --output_biom_fp cdout_noless100/taxa_plots_Group/Group_otu_table.biom --output_mapping_fp cdout_noless100/taxa_plots_Group/Group_map.txt --collapse_fields 'Group' Stdout: Stderr: # Sort OTU Table command sort_otu_table.py -i cdout_noless100/taxa_plots_Group/Group_otu_table.biom -o cdout_noless100/taxa_plots_Group/Group_otu_table_sorted.biom Stdout: Stderr: # Summarize Taxonomy command summarize_taxa.py -i cdout_noless100/taxa_plots_Group/Group_otu_table_sorted.biom -o cdout_noless100/taxa_plots_Group/ Stdout: Stderr: # Plot Taxonomy Summary command plot_taxa_summary.py -i cdout_noless100/taxa_plots_Group/Group_otu_table_sorted_L2.txt,cdout_noless100/taxa_plots_Group/Group_otu_table_sorted_L3.txt,cdout_noless100/taxa_plots_Group/Group_otu_table_sorted_L4.txt,cdout_noless100/taxa_plots_Group/Group_otu_table_sorted_L5.txt,cdout_noless100/taxa_plots_Group/Group_otu_table_sorted_L6.txt -o cdout_noless100/taxa_plots_Group//taxa_summary_plots/ Stdout: Stderr: /usr/lib64/python2.7/site-packages/matplotlib/cbook.py:136: MatplotlibDeprecationWarning: The axisbg attribute was deprecated in version 2.0. Use facecolor instead. warnings.warn(message, mplDeprecation, stacklevel=1) Executing commands. # Collapse samples in OTU table by categories command collapse_samples.py -m Diet_mapping.tsv -b cdout_noless100/table_mc2000.biom --output_biom_fp cdout_noless100/taxa_plots_Treatment/Treatment_otu_table.biom --output_mapping_fp cdout_noless100/taxa_plots_Treatment/Treatment_map.txt --collapse_fields 'Treatment' Stdout: Stderr: # Sort OTU Table command sort_otu_table.py -i cdout_noless100/taxa_plots_Treatment/Treatment_otu_table.biom -o cdout_noless100/taxa_plots_Treatment/Treatment_otu_table_sorted.biom Stdout: Stderr: # Summarize Taxonomy command summarize_taxa.py -i cdout_noless100/taxa_plots_Treatment/Treatment_otu_table_sorted.biom -o cdout_noless100/taxa_plots_Treatment/ Stdout: Stderr: # Plot Taxonomy Summary command plot_taxa_summary.py -i cdout_noless100/taxa_plots_Treatment/Treatment_otu_table_sorted_L2.txt,cdout_noless100/taxa_plots_Treatment/Treatment_otu_table_sorted_L3.txt,cdout_noless100/taxa_plots_Treatment/Treatment_otu_table_sorted_L4.txt,cdout_noless100/taxa_plots_Treatment/Treatment_otu_table_sorted_L5.txt,cdout_noless100/taxa_plots_Treatment/Treatment_otu_table_sorted_L6.txt -o cdout_noless100/taxa_plots_Treatment//taxa_summary_plots/ Stdout: Stderr: /usr/lib64/python2.7/site-packages/matplotlib/cbook.py:136: MatplotlibDeprecationWarning: The axisbg attribute was deprecated in version 2.0. Use facecolor instead. warnings.warn(message, mplDeprecation, stacklevel=1) Executing commands. # Collapse samples in OTU table by categories command collapse_samples.py -m Diet_mapping.tsv -b cdout_noless100/table_mc2000.biom --output_biom_fp cdout_noless100/taxa_plots_Time/Time_otu_table.biom --output_mapping_fp cdout_noless100/taxa_plots_Time/Time_map.txt --collapse_fields 'Time' Stdout: Stderr: # Sort OTU Table command sort_otu_table.py -i cdout_noless100/taxa_plots_Time/Time_otu_table.biom -o cdout_noless100/taxa_plots_Time/Time_otu_table_sorted.biom Stdout: Stderr: # Summarize Taxonomy command summarize_taxa.py -i cdout_noless100/taxa_plots_Time/Time_otu_table_sorted.biom -o cdout_noless100/taxa_plots_Time/ Stdout: Stderr: # Plot Taxonomy Summary command plot_taxa_summary.py -i cdout_noless100/taxa_plots_Time/Time_otu_table_sorted_L2.txt,cdout_noless100/taxa_plots_Time/Time_otu_table_sorted_L3.txt,cdout_noless100/taxa_plots_Time/Time_otu_table_sorted_L4.txt,cdout_noless100/taxa_plots_Time/Time_otu_table_sorted_L5.txt,cdout_noless100/taxa_plots_Time/Time_otu_table_sorted_L6.txt -o cdout_noless100/taxa_plots_Time//taxa_summary_plots/ Stdout: Stderr: /usr/lib64/python2.7/site-packages/matplotlib/cbook.py:136: MatplotlibDeprecationWarning: The axisbg attribute was deprecated in version 2.0. Use facecolor instead. warnings.warn(message, mplDeprecation, stacklevel=1) Executing commands. # Boxplots (Diet) command make_distance_boxplots.py -d cdout_noless100/bdiv_even2000//weighted_unifrac_dm.txt -f Diet -o cdout_noless100/bdiv_even2000//weighted_unifrac_boxplots/ -m Diet_mapping.tsv -n 999 Stdout: Stderr: # Boxplots (DayofDiet) command make_distance_boxplots.py -d cdout_noless100/bdiv_even2000//weighted_unifrac_dm.txt -f DayofDiet -o cdout_noless100/bdiv_even2000//weighted_unifrac_boxplots/ -m Diet_mapping.tsv -n 999 Stdout: Stderr: # Boxplots (Group) command make_distance_boxplots.py -d cdout_noless100/bdiv_even2000//weighted_unifrac_dm.txt -f Group -o cdout_noless100/bdiv_even2000//weighted_unifrac_boxplots/ -m Diet_mapping.tsv -n 999 Stdout: Stderr: # Boxplots (Treatment) command make_distance_boxplots.py -d cdout_noless100/bdiv_even2000//weighted_unifrac_dm.txt -f Treatment -o cdout_noless100/bdiv_even2000//weighted_unifrac_boxplots/ -m Diet_mapping.tsv -n 999 Stdout: Stderr: # Boxplots (Time) command make_distance_boxplots.py -d cdout_noless100/bdiv_even2000//weighted_unifrac_dm.txt -f Time -o cdout_noless100/bdiv_even2000//weighted_unifrac_boxplots/ -m Diet_mapping.tsv -n 999 Stdout: Stderr: # Boxplots (Diet) command make_distance_boxplots.py -d cdout_noless100/bdiv_even2000//unweighted_unifrac_dm.txt -f Diet -o cdout_noless100/bdiv_even2000//unweighted_unifrac_boxplots/ -m Diet_mapping.tsv -n 999 Stdout: Stderr: # Boxplots (DayofDiet) command make_distance_boxplots.py -d cdout_noless100/bdiv_even2000//unweighted_unifrac_dm.txt -f DayofDiet -o cdout_noless100/bdiv_even2000//unweighted_unifrac_boxplots/ -m Diet_mapping.tsv -n 999 Stdout: Stderr: # Boxplots (Group) command make_distance_boxplots.py -d cdout_noless100/bdiv_even2000//unweighted_unifrac_dm.txt -f Group -o cdout_noless100/bdiv_even2000//unweighted_unifrac_boxplots/ -m Diet_mapping.tsv -n 999 Stdout: Stderr: # Boxplots (Treatment) command make_distance_boxplots.py -d cdout_noless100/bdiv_even2000//unweighted_unifrac_dm.txt -f Treatment -o cdout_noless100/bdiv_even2000//unweighted_unifrac_boxplots/ -m Diet_mapping.tsv -n 999 Stdout: Stderr: # Boxplots (Time) command make_distance_boxplots.py -d cdout_noless100/bdiv_even2000//unweighted_unifrac_dm.txt -f Time -o cdout_noless100/bdiv_even2000//unweighted_unifrac_boxplots/ -m Diet_mapping.tsv -n 999 Stdout: Stderr: # Compare alpha diversity (PD_whole_tree) command compare_alpha_diversity.py -i cdout_noless100/arare_max2000//alpha_div_collated/PD_whole_tree.txt -m Diet_mapping.tsv -c Diet,DayofDiet,Group,Treatment,Time -o cdout_noless100/arare_max2000//compare_PD_whole_tree -n 999 Stdout: Stderr: # Compare alpha diversity (chao1) command compare_alpha_diversity.py -i cdout_noless100/arare_max2000//alpha_div_collated/chao1.txt -m Diet_mapping.tsv -c Diet,DayofDiet,Group,Treatment,Time -o cdout_noless100/arare_max2000//compare_chao1 -n 999 Stdout: Stderr: # Compare alpha diversity (observed_otus) command compare_alpha_diversity.py -i cdout_noless100/arare_max2000//alpha_div_collated/observed_otus.txt -m Diet_mapping.tsv -c Diet,DayofDiet,Group,Treatment,Time -o cdout_noless100/arare_max2000//compare_observed_otus -n 999 Stdout: Stderr: # Group significance (Diet) command group_significance.py -i cdout_noless100/table_even2000.biom -m Diet_mapping.tsv -c Diet -o cdout_noless100/group_significance_Diet.txt Stdout: Stderr: # Group significance (DayofDiet) command group_significance.py -i cdout_noless100/table_even2000.biom -m Diet_mapping.tsv -c DayofDiet -o cdout_noless100/group_significance_DayofDiet.txt Stdout: Stderr: # Group significance (Group) command group_significance.py -i cdout_noless100/table_even2000.biom -m Diet_mapping.tsv -c Group -o cdout_noless100/group_significance_Group.txt Stdout: Stderr: # Group significance (Treatment) command group_significance.py -i cdout_noless100/table_even2000.biom -m Diet_mapping.tsv -c Treatment -o cdout_noless100/group_significance_Treatment.txt Stdout: Stderr: # Group significance (Time) command group_significance.py -i cdout_noless100/table_even2000.biom -m Diet_mapping.tsv -c Time -o cdout_noless100/group_significance_Time.txt Stdout: Stderr: # Compress the filtered BIOM table command gzip cdout_noless100/table_mc2000.biom Stdout: Stderr: # Compress the rarefied BIOM table command gzip cdout_noless100/table_even2000.biom Stdout: Stderr: Logging stopped at 12:53:10 on 03 Aug 2017