利用leafcutter做可变剪切(新手向)

本文参考于生信技能树
利用leafcutter做可变剪切需要的材料是RNA-seq流程中所获得的bam类型的文件
首先是软件的安装:
1.源代码下载
git clone https://github.com/davidaknowles/leafcutter
#这条命令下载该软件的源代码,需要修改里面的一个脚本 scripts/bam2junc.sh 把软件路径增添进去即可,并将该软件路径添加到环境变量中
2.软件安装:
由于该软件是一个R包所以可以直接在R中安装这个包:
if (!require(“devtools”))
install.packages(“devtools”, repos=‘http://cran.us.r-project.org’)
devtools::install_github(“davidaknowles/leafcutter/leafcutter”)

另外也可以使用conda安装:
conda install -y -c davidaknowles r-leafcutter(我没成功)

正式流程
1.将bam文件转换位junc类型的文件:
*ls .bam >bam_dir.txt #构建bam的路径文件
cat bam_dir.txt|while read id; do file=$(basename id);sample=id ); sample=i</

评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值