【原创】hisat2 mapping转录组分析笔记2018-09-27

先下载人类基因组index

cd referece 

mkdir index 

cd index

wget -c ftp://ftp.ccb.jhu.edu/pub/infphilo/hisat2/data/hg19.tar.gz

tar -zxvf hg19.tar.gz

下载的gz文件先解压,得到好多文件,要用cat合在一起。

cat genome.1.ht2 genome.2.ht2 genome.3.ht2 genome.4.ht2 genome.5.ht2 genome.6.ht2 genome.7.ht2 genome.8.ht2 > genome

hisat2 -t -x~/rna_seq_analysis/reference/index/hg19/genome -1 ~/rna_seq_analysis/fastq/SRR5894154_1.fastq.gz -2 ~/rna_seq_analysis/fastq/SRR5894154_2.fastq.gz -s ~/rna_seq_analysis/aligned/SRR5894154.sam

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容