基于USEARCH 进行Alpha,Beta,进化树等分析操作
上回收到,使用USEARCH处理完双端测序数据,生成OTU丰度表及代表性序列。这里,继续使用USEARCH进行微生物的Alpha、Beta、进化树等分析操作。主要分为两个部分,第一部分为处理双端测序的脚本修正,第二部分为Alpha、Beta、进化树等分析操作分析脚本。
第一部分:
#!/bin/bash
### 使用 vsearch+usearch11+QIIME 处理双端测序数据
### 目标 OTU表和注释信息
## 目录
# temp/: 临时文件目录
# res/: 结果文件目录
# data/:双端测序数据集
# database/:参考数据库
# err/: 错误信息
# log/: 操作日志
## 文件
HOME=$HOME/Lab_2 ### 操作空间
temp=${HOME}/temp2 ### 临时文件存放目录
res=${HOME}/res2 ### 结果文件存放目录
data=${HOME}/seq ### 测序数据存放位置
refdb=${HOME}/database ### 参考测序数据库
err=${HOME}/err2 ### 错误信息
log=${HOME}/log2 ### 日志信息
###########################################################################################
# 实验准备