usage: htseq-count [options] alignment_file gff_file
-f {sam,bam} (default: sam)
-r {pos,name} (default: name)
-s {yes,no,reverse} (default: yes) #此处关于选项-s为我自己的认识，不一定对
#数据是否来源于链特异性测序，链特异性是指在建库测序时，只测mRNA反转录出的cDNA序列，而不测该cDNA序列反向互补的另一条DNA序列；换句话说就是，链特异性能更准确反映出mRNA的序列信息
#我们知道在gff/gtf中第7列是+-信息，+表示来源于参考基因组序列正链，-表示参考基因组序列的反向互补链
#sam/bam文件的第2列是flag信息，也可以看出比对到正链还是负链
#stranded=no，无链特异性，一条reads通过flag列知道比对到+还是-链后，不管是不是和gff/gtf相匹配，都算是这个feature中的
#stranded=yes, 且se测序，要和gff/gtf相匹配，才算是这个feature中的
#stranded=yes, 且pe测序，要和gff/gtf相匹配，才算是这个feature中的
#stranded=reverse，是yes的相反，这时不是和gff/gtf相匹配了，而是恰好相反，可能源于另一种链特异性，只测cDNA序列反向互补的另一条DNA序列
-a MINAQUAL (default: 10)
#忽略比对质量低于此值的比对结果
-t FEATURETYPE
#feature type (3rd column in GFF file) to be used, all features of other type are ignored (default, suitable for Ensembl GTF files: exon)
#没想到这个还能自己设置
-i IDATTR
#GFF attribute to be used as feature ID (default, suitable for Ensembl GTF files: gene_id)
-m {union,intersection-strict,intersection-nonempty} (default: union)
usage: htseq-count
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 参考基因组下载有三大全文网站提供参考基因组下载，它们分别是：1.NCBI （https://www.ncbi.nl...
- 理论基础 在上篇的比对中，我们需要纠结是否真的需要比对，如果你只需要知道已知基因的表达情况，那么可以选择align...
- # Because the id list get from htseq-count was exactly sa...
- htseq-count计数的相关内容前面在不同的学习阶段写过两次，分别是合并htseq-count的结果得到cou...
- 和featurecounts一样，htseq-count也是一款进行raw count定量的软件。该软件采用pyt...
