使用.bed或.bam进行Peakcalling时MACS2的结果不同

时间:2015-11-10 16:45:47

标签: bioinformatics

我遇到了以下问题: 我使用MACS2(2.1.0.20140616)和以下短命令行:

macs2 callpeak -t file.bam -f bam -g 650000000 -n Test -B --nomodel -q 0.01

它似乎按我的意愿工作,但是当我通过

将.bamfile转换为.bed时
bedtools bamtobed -i file.bam > file.bed

并在此使用MACS2,我得到了更多的峰值。据我所知,.bed文件应该包含与.bam文件相同的信息,所以这有点奇怪。

有什么建议是什么问题?

谢谢!

0 个答案:

没有答案