软件使用方法如下:
#hisat2-build 对参考基因组构建索引
hisat2-build genome.fa genome.fa
#双端read比对
hisat2 --new-summary -p 10 -x genome.fa -1 1.fastq.gz -2 2.fastq.gz -S rnaseq.sam
#单端read比对
hisat2 --new-summary -p 10 -x genome.fa -U rnaseq.fastq.gz -S rnaseq.sam
#排序并生成bam文件
samtools sort -o rnaseq.bam rnaseq.sam
有比对的bam文件,就可以进行后续分析了。
此外,我们在网易云课堂上有各种教学视频,有兴趣可以了解一下:
1. 文章越来越难发?是你没发现新思路,基因家族分析发2-4分文章简单快速,学习链接:基因家族分析实操课程
2. 转录组数据理解不深入?图表看不懂?点击链接学习深入解读数据结果文件,学习链接:转录组(有参)结果解读;转录组(无参)结果解读
3. 转录组数据深入挖掘技能-WGCNA,提升你的文章档次,学习链接:WGCNA-加权基因共表达网络分析
4. 转录组数据怎么挖掘?学习链接:转录组标准分析后的数据挖掘
6. 更多学习内容:linux、perl、R语言画图,更多免费课程请点击以下链接:
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!