美文网首页
RSeQC鉴定文库是否为链特异性文库以及其他小功能

RSeQC鉴定文库是否为链特异性文库以及其他小功能

作者: njmujjc | 来源:发表于2019-03-20 20:18 被阅读0次

首先将基因组gtf文件转换为bed:

perl /media/shen/disk3/biosoft/gtf2bed/gtf2bed.pl gtf文件 > hg19.bed

脚本下载地址:https://github.com/ExpressionAnalysis/ea-utils/blob/master/clipper/gtf2bed

conda下载RSeQC:

source activate py2.7

conda install -c bioconda rseqc

infer_experiment.py -r hg19.bed -i align.bam

2. gene body coverage

geneBody_coverage.py -i 'Br_RNA1.bam' -r '/media/shen/disk1/jjc/reference/hg19/GRCh37-ERCC.bed' -o Br_RNA1_genebodycov

相关文章

网友评论

      本文标题:RSeQC鉴定文库是否为链特异性文库以及其他小功能

      本文链接:https://www.haomeiwen.com/subject/fbxxvqtx.html