kraken注释物种，eggnog注释基因一条龙服务

cd /mnt/10t/database/ kraken2-build –download-library bacteria –db krakenbac –use-ftp kraken2-build –download-library archaea –db krakenarch –use-ftp mkdir /mnt/10t/mzy/24samples/tax cd /mnt/10t/mzy/24samples/tax for i in 123 126 134 131 125 140 132 149 148 122 143 133 145 152 118 147 121 151 135 136 150 146 141 137 do kraken2 –db /mnt/10t/database/krakenbac –report i.DNAbacreport –use-mpa-style –report-zero-counts –threads 80 –gzip-compressed –paired /mnt/10t/mzy/72sample/DNAclip/i.clip.1.fq.gz /mnt/10t/mzy/72sample/DNAclip/i.clip.2.fq.gz > i.DNA.sam kraken2 –db /mnt/10t/database/krakenarch –report i.DNAmetreport –use-mpa-style –report-zero-counts –threads 80 –gzip-compressed –paired /mnt/10t/mzy/72sample/DNAclip/i.clip.1.fq.gz /mnt/10t/mzy/72sample/DNAclip/i.clip.2.fq.gz > i.DNA.sam && rm -f

kraken2 –db /mnt/10t/database/krakenbac –report i.RNAbacreport –use-mpa-style –report-zero-counts –threads 80 –gzip-compressed –paired /mnt/10t/mzy/rawdata/zhuanluzu/rawzhuanlu/i\_R1.fq.gz /mnt/10t/mzy/rawdata/zhuanluzu/rawzhuanlu/i\_R2.fq.gz > i.RNA.sam kraken2 –db /mnt/10t/database/krakenarch –report i.RNAmetreport –use-mpa-style –report-zero-counts –threads 80 –gzip-compressed –paired /mnt/10t/mzy/rawdata/zhuanluzu/rawzhuanlu/i\_R1.fq.gz /mnt/10t/mzy/rawdata/zhuanluzu/rawzhuanlu/i\_R2.fq.gz > i.RNA.sam && rm -f

版权声明：本文内容由互联网用户自发贡献，该文观点仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 举报，一经查实，本站将立刻删除。

发布者：全栈程序员栈长，转载请注明出处：https://javaforall.cn/227655.html原文链接：https://javaforall.cn