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 && 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