sort
按照默认参数,进行htseq-count
统计reads数的时候会报错:
image.png
image.png
看了 htseq-count
的参数后发现htseq-count
默认是经name sorted
的文件,而samtools -sort
的时候默认参数是按照position
,所以出现报错很正常。
我拿按照-name sorted
过的bam文件去做htseq-count
的时候虽然没有报错啦,但是建index的时候出错,如下图
三脸懵逼中,还是不是很明白。
不管报错,继续跑完程序,发现按照name sort和按照position sort,统计出来的read counts 差了近一倍。
image.png
找到了类似的报错:Samtools sorting and indexing
网友评论