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

大家好，又见面了，我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46，售后保障稳定

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 && rm -f $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 $i.DNA.sam

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 && rm -f $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 $i.RNA.sam
cd /mnt/10t/mzy/24samples/geneset
/usr/bin/python2.7 /mnt/10t/eggnog-mapper/emapper.py -m diamond –no_annot –no_file_comments –cpu 10 -i /mnt/10t/mzy/24samples/geneset/unique24.faa -o eggnog.txt