本文小编为大家详细介绍“picrust2如何使用”,内容详细,步骤清晰,细节处理妥当,希望这篇“picrust2如何使用”文章能帮助大家解决疑惑,下面跟着小编的思路慢慢深入,一起来学习新知识吧。
picrust2用法
#去掉序列ID中注释信息,避免错误:no ASV ids overlap between input FASTA and sequence abundance table
get_fa_by_id.pl $workdir/5.pick_otu_qiime/pick_de_novo_otus/otu_table_clean.txt $workdir/5.pick_otu_qiime/pick_de_novo_otus/rep_set/qiime_rep_set.fasta rep_set.fa rm -rf picrust2_out picrust2_pipeline.py -s rep_set.fa -i $workdir/5.pick_otu_qiime/pick_de_novo_otus/otu_table.biom -o picrust2_out --processes 1 --in_traits COG,EC,KO,PFAM,TIGRFAM add_descriptions.py -i picrust2_out/EC_metagenome_out/pred_metagenome_unstrat.tsv.gz -m EC \ -o picrust2_out/EC_metagenome_out/pred_metagenome_unstrat_descrip.tsv.gz add_descriptions.py -i picrust2_out/KO_metagenome_out/pred_metagenome_unstrat.tsv.gz -m KO \ -o picrust2_out/KO_metagenome_out/pred_metagenome_unstrat_descrip.tsv.gz add_descriptions.py -i picrust2_out/COG_metagenome_out/pred_metagenome_unstrat.tsv.gz -m COG \ -o picrust2_out/COG_metagenome_out/pred_metagenome_unstrat_descrip.tsv.gz add_descriptions.py -i picrust2_out/PFAM_metagenome_out/pred_metagenome_unstrat.tsv.gz -m PFAM \ -o picrust2_out/PFAM_metagenome_out/pred_metagenome_unstrat_descrip.tsv.gz add_descriptions.py -i picrust2_out/TIGRFAM_metagenome_out/pred_metagenome_unstrat.tsv.gz -m TIGRFAM \ -o picrust2_out/TIGRFAM_metagenome_out/pred_metagenome_unstrat_descrip.tsv.gz add_descriptions.py -i picrust2_out/pathways_out/path_abun_unstrat.tsv.gz -m METACYC \ -o picrust2_out/pathways_out/path_abun_unstrat_descrip.tsv.gz
读到这里,这篇“picrust2如何使用”文章已经介绍完毕,想要掌握这篇文章的知识点还需要大家自己动手实践使用过才能领会,如果想了解更多相关内容的文章,欢迎关注亿速云行业资讯频道。
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。