美文网首页
samtools sort

samtools sort

作者: 超级无敌大蜗牛 | 来源:发表于2020-02-23 15:31 被阅读0次

sort 按照默认参数,进行htseq-count统计reads数的时候会报错:

image.png
image.png
image.png

看了 htseq-count的参数后发现htseq-count默认是经name sorted的文件,而samtools -sort的时候默认参数是按照position,所以出现报错很正常。
我拿按照-name sorted过的bam文件去做htseq-count的时候虽然没有报错啦,但是建index的时候出错,如下图

image.png

三脸懵逼中,还是不是很明白。

不管报错,继续跑完程序,发现按照name sort和按照position sort,统计出来的read counts 差了近一倍。


image.png

找到了类似的报错:Samtools sorting and indexing

相关文章

网友评论

      本文标题:samtools sort

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