一、参考基因组:https://hgdownload.soe.ucsc.edu/goldenPath/oviAri4/bigZips/oviAri4.fa.gz
二、GTF:https://hgdownload.soe.ucsc.edu/goldenPath/oviAri4/bigZips/genes/oviAri4.ncbiRefSeq.gtf.gz
三、hisat2 构建INDEX流程:
1.hisat2_extract_exons.py oviAri4.ncbiRefSeq.gtf > oviAri4.ncbiRefSeq.gtf.exon 2.hisat2_extract_splice_sites.py oviAri4.ncbiRefSeq.gtf > oviAri4.ncbiRefSeq.ss &
3.hisat2-build -p 6 --exon oviAri4.ncbiRefSeq.gtf.exon#提取外显子# --ss oviAri4.ncbiRefSeq.ss#提取剪切位点# oviAri4.fa#参考基因组# oviAri4.fa_ex_ss.gtf > hisat2_build.log 2>&1 &
四、DESeq2得到差异表达基因后,log2foldchange的绝对值> 1 同时padj值小于0.05的筛选条件,共获得了5552个差异表达基因。
想请问一下这么多的差异表达基因可以改筛选条件么?
这家伙很懒,还没有设置简介