I'm facing a problem calling peaks from a eCLIP experiment. I have 3 bam files: 2 replicates for one condition and a negative control. I don't know if this is the best option but I'm calling the peaks using MACS (from Galaxy), which is supposed to be for ChiP-seq and, for this reason, it doesn't give me the strand of the signal in the output file but I really need this information. Is there a way to tell MACS that I want the strand information to be kept? Or maybe is there a better way to call peaks for a CLIP experiment?
Once I have the bed files with the peaks I will use bedtools to extract the RNAs were I can find a peak for different factors.
Thank you in advance for your help