本文最后更新于99 天前,其中的信息可能已经过时,如有错误请发送邮件到wuyk@163.com
1.测序数据
甘蓝黄花YL-1,YL-2,芥兰白花11-192-1,11-192-2
2.软件及功能
Fastqc and fastp: 质控和去除测序原始数据的接头
Hisat2: 基于参考基因组建立索引
Samtools: 排序
Stringtie: 计算基因表达水平
featureCounts: 计算Counts值
DESeq2: 基因差异表达分析
3.RNA-seq的详细流程
3.1 Fastqc and Fastp
#fastqc Quality Control 会输出网页的质控报告
for i in {1..2}; do fastqc YL1-${i}_R1.fq.gz;done
for i in {1..2}; do fastqc YL1-${i}_R2.fq.gz;done
for i in {1..2}; do fastqc A192-${i}_R1.fq.gz;done
for i in {1..2}; do fastqc A192-${i}_R2.fq.gz;done
#fastp 进行过滤操作,其中 -i 和 -I 表示测序数据的两个端点,而 -o和 -O 则表示已过滤的测序数据
for i in {1..2}; do fastp -i YL1-${i}_R1.fq.gz -I YL1-${i}_R2.fq.gz -o YL1-${i}.R1.fq.ft.gz -O YL1-${i}.R2.fq.ft.gz
for i in {1..2}; do fastp -i A192-${i}_R1.fq.gz -I A192-${i}_R2.fq.gz -o A192-${i}.R1.fq.ft.gz -O A192-${i}.R2.fq.ft.gz