美文网首页走进转录组简友广场
RNA-seq上游分析(3)-使用featureCounts对b

RNA-seq上游分析(3)-使用featureCounts对b

作者: 灵活胖子的进步之路 | 来源:发表于2020-11-25 08:48 被阅读0次

使用featureCounts对bam文件进行定量。

## 定义输入输出文件夹
gtf=/home/zk/project/reference/annotateGTF/Homo_sapiens.GRCh38.101.gtf.gz
inputdir=/media/zk/crlm/map
outputdir=/media/zk/crlm/count

# featureCounts对bam文件进行计数
nohup featureCounts -T 16 -p -t exon -g gene_id -a $gtf -o all.id.txt $inputdir/*.sorted.bam >count.log 2>&1 &

# 对定量结果质控
nohup multiqc all.id.txt.summary >countqc.log 2>&1 &

# 得到表达矩阵
cat all.id.txt | cut -f1,7- > counts.txt

相关文章

网友评论

    本文标题:RNA-seq上游分析(3)-使用featureCounts对b

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