美文网首页
featureCounts 统计reads数量

featureCounts 统计reads数量

作者: 斩毛毛 | 来源:发表于2020-07-20 10:28 被阅读0次

此前,介绍过在转录组分析过程中使用HTseq统计reads数量;HTseq统计reads数量。本次介绍一款速度较于HTseq快很多的reads 数量统计软件。

安装

wget https://sourceforge.net/projects/subread/files/subread-2.0.1/subread-2.0.1-Linux-x86_64.tar.gz
tar zxf subread-2.0.1-Linux-x86_64.tar.gz

添加环境变量即可

简单说明

featureCounts -T 6 -p -t exon -g gene_id -a test.annotation.gtf -o test_featureCounts.txt test.sorted.bam

参数说明
-a: 输入GTF/GFF注释文件
-p: 针对paired-end reads
-F: -a参数的是个,默认为GTF
-g: 选择哪个部分作为reads count;默认为gene_id
-t:  与-g类似,默认为exon
-o:输出文件
-T:线程数

后续则可进行差异分析

相关文章

网友评论

      本文标题:featureCounts 统计reads数量

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