继续前面的学习，前面已经把软件安装完成，数据库准备好，下面就是分析的过程了，基本上按照原文的命令进行的，由于教程中没有给出tara_f135_full_megahit.fasta这个文件，这里我们就把这几个样本的组装提到了前面，自己组装获得这个序列，然后再进行物种注释。
前面提到，已经有家公司通过宏转录组(Metatranscriptomics)测序检测肠道微生物，面向消费者提供检测服务。对宏转录组充满了好奇，有这样的比方说，**宏基因组可以告诉我们这个微生物群落可能有什么样的功能（潜能），宏转录组就是告诉我们群落正在做什么**，相比宏基因组的眉毛胡子一把抓，宏转录组是更针对当下的结果。由于测序的目标序列少了很多，结果不是变态大，对计算机的配置要求也相对降低。苦于想学宏基因组暂时没有服务器的我，就退而求其次试试宏转录组了，相信不会让我失望。之前学习过单转录组数据的分析，一般的笔记本（双核，8g ram）扛了下来。鉴于中文网络上能找到的宏转录组教程基本没有，只在Github上搜索到两个，选其中一个学习下。
本教程将带您完成处理宏转录组数据的流程。
到目前为止，我们已经完成了所有工作，并复制并粘贴了许多命令来完成所需的操作。这可行！但是也可能很耗时，并且更容易出错。接下来，我们将向你展示如何将所有这些命令放入Shell脚本中。