kraken2注释+braken格式转换

由于最近遇到一个宏基因项目，第一次分析注释到的物种数目较少，因此更换数据库试试看看！

于是kraken上榜，OK，let's  try 一下子！

前提是已经安装了kraken和braken，运行命令比较简单

1.注释

##基于cleanreads进行注释：
cd /kraken2/test_kraken2
mkdir sample1_3
/kraken2/kraken2-2.1.2/kraken2  \
	--db /kraken2/database/ \
	--paired  --threads 24 \
	/cleandata/sample1_3_R1.clean.fq.gz \
	/cleandata/sample1_3_R2.clean.fq.gz \
	--output sample1_3/GOE1_3.kraken --report sample1_3.report

查看一下sample1_3.report文件格式：

第一列是相似百分比，最后一列物种，倒数第二列是层级

共有8个样本，预期得到的格式为：行为物种，列为样本名，中间值为物种丰度，如图所示

上代码：

2.格式修改<