宏基因组分析流程

论坛 期权论坛 脚本     
匿名技术用户   2020-12-21 21:13   44   0
1.md5sum+trimmomatic
md5sum SRR1976948_1.fastq.gz SRR1976948_2.fastq.gz

java -jar /data/XXXXX/software/software/Trimmomatic-0.36/trimmomatic-0.36.jar PE \
-phred33 SRR1976948_1.fastq.gz SRR1976948_2.fastq.gz \
/data/XXXXX/test/MGS/01trim/SRR1976948_1_paired.fq.gz \
/data/XXXXX/test/MGS/01trim/SRR1976948_1_unpaired.fq.gz \
/data/XXXXX/test/MGS/01trim/SRR1976948_2_paired.fq.gz \
/data/XXXXX/test/MGS/01trim/SRR1976948_2_unpaired.fq.gz \
ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36

Remove adapters (ILLUMINACLIP:TruSeq3-PE.fa:2:30:10)
Remove leading low quality or N bases (below quality 3) (LEADING:3)
Remove trailing low quality or N bases (below quality 3) (TRAILING:3)
Scan the read with a 4-base wide sliding window, cutting when the average quality per base drops below 15 (SLIDINGWINDOW:4:15)
Drop reads below the 36 bases long (MINLEN:36)

2.fastqc
fastqc SRR1976948_1_paired.fq.gz
fastqc SRR1976948_2_paired.fq.gz
3.MEGAHIT

4.QUEST评估组装效果

5.Prokka注释基因

6.sourmash比较数据集教程

7.基因丰度估计Salmon

8.分箱宏基因组

分享到 :
0 人收藏
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

积分:7942463
帖子:1588486
精华:0
期权论坛 期权论坛
发布
内容

下载期权论坛手机APP