测序比对软件的总结-star

这篇博客简要介绍了如何使用STAR软件进行测序比对。首先通过conda安装STAR,接着下载基因组注释的GTF和FASTA数据来生成索引文件。然后进行比对操作,最终得到的比对结果是bam格式,通过samtools将其转换为sam格式。
摘要由CSDN通过智能技术生成

github的说明书看这里:https://github.com/alexdobin/STAR/blob/master/doc/STARmanual.pdf

这里只做实际操作的简单介绍,里面的具体细节需要在github里面看。

1.安装
使用conda进行安装

conda install -c bioconda star

2.下载基因组的注释数据
GTF

wget https://ftp.ebi.ac.uk/pub/databases/gencode/Gencode_human/release_40/gencode.v40.annotation.gtf.gz

FASTA

wget https://ftp.ebi.ac.uk/pub/databases/gencode/Gencode_human/release_40/GRCh38.p13.genome.fa.gz

3.生成索引文件

STAR  --runMode genomeGenerate --genomeDir /star-h38-index\#输出的index的文件夹
 --genomeFastaFiles ./GRCh38.p13.genome.fa \#fasta文件
 --sjdbGTFfile ./gencode.v40.annotation.gtf\#gtf文件

生成结果如下:

(riboseq) [xx@login01 star-h38-index]$ ll -h
total 30G
-rw-r--r-- 1 med-zhouh med-chenh 4.2K Jul 11 18:10 chrLength.txt
-rw-r--r-- 1 med-zhouh med-chenh  11K Jul 11 18:10 chrNameLength.txt
-rw-r--r-- 1 med-zhouh med-chenh 6.8K Jul 11 18:10 chrName.txt
-rw-r--r-- 1 med-zhouh med-chenh 6.9K Jul 11 18:10 chrStart.txt
-rw-r--r-- 1 med-zhouh med-chenh  54M Jul 11 18:10 exonGeTrInfo.tab
-rw-r--r-- 
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值