使用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
网友评论