I am trying to analyze PAR-CLIP data to identify interactions between cellular RNAs and proteins of interest. I have trimmed adapters from my sequence data and aligned the files using Bowtie. What is the best way to:
(a) Filter aligned reads so that only those that contain diagnostic T>C mutations are retained?
(b) Identify protein binding motifs using the filtered reads?