2

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

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

它似乎可以按我的意愿工作,但是当我通过 .bamfile 将 .bamfile 转换为 .bed 时

bedtools bamtobed -i file.bam > file.bed

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

有什么建议有什么问题吗?

谢谢!

4

0 回答 0