Entering edit mode
7.2 years ago
rob.costa1234
▴
310
I have a list of 10 genes that i want to find level of enrichment in a public data set. So I downloaded FAStq files and run the alignment with Bowtie, Macs for peaks. Now my question is how should i Specifically check the status of these selected 10 genes in analyzed data. I can take an indirect approach that annotate MACS peaks first and see overlap of the selected genes. I was wondering is there any direct approach of sub-setting Peaks data corresponding to my 10 selected genes. I only have gene ids- hugo gene names.
Thanks for the input.
Assuming that you mean enrichment OF some transcription factor/histone modification ON the preselected genes here is one thing you can do.
Just upload the peak.bed file from MACS output and the alignment file (SAM/BAM) directly onto some genome browser. Again, assuming that the genome browser that you use, also has the genome and gene list annotated (else you may have to do it by yourself), you can just search for your preselected genes and take a look at the peak profile on each genomic locus one-by-one.