A question to all users that run CLIP/RIP-seq and then try to find peaks in their data.
Since I couldn't find any specific tool for finding peaks in RNA-based data, I have decided to use the classical ChIP-seq tools (the logic is the same I think).
I align reads to the genome (I don't really care about junction reads), and then run one of the ChIP-seq tools such as MACS.
I would like to hear about the parameters people use when they run ChIP-seq tools on CLIP/RIP-seq data. For example many of the ChIP-seq tools use the "genome size" for random expectation calculations, and different genome sizes give different output.
Do you use the transcriptome size or genome size when trying to find significant peaks?
Any other issues that should be taken care of?
Thanks
Mali
Since I couldn't find any specific tool for finding peaks in RNA-based data, I have decided to use the classical ChIP-seq tools (the logic is the same I think).
I align reads to the genome (I don't really care about junction reads), and then run one of the ChIP-seq tools such as MACS.
I would like to hear about the parameters people use when they run ChIP-seq tools on CLIP/RIP-seq data. For example many of the ChIP-seq tools use the "genome size" for random expectation calculations, and different genome sizes give different output.
Do you use the transcriptome size or genome size when trying to find significant peaks?
Any other issues that should be taken care of?
Thanks
Mali
Comment