3.1 years ago by
I believe when you extract the list of genes from your data you have the strand specificity right? so then you will be able to understand which genes correspond to which strand be it + or - thus giving you strand specific feature. Then you can grep your output based on strand features.
This will give you two lists of promoters that have either + or - strandedness. Once you have it when you can overlap the genes to see bidirectional genes , since those which will overlap at refeseqIDs or gene symbols should be shared at the level of both strands. I believe this will help.