Is removing strand information when obtaining intergenic regions from GRanges necessary?

Is removing strand information when obtaining intergenic regions from GRanges necessary?

0

Hello Biostars Community,

I am attempting to obtain bed file of intergenic regions from a custom GenomicRanges object.

I am following this tutorial:

research.stowers.org/cws/CompGenomics/Tutorial/peak_assignment.html#orgb42e0a8

The tutorials guides me to remove strand information:

# remove strand information
strand(genic) <- '*'

Could someone “explain like I’m five” why this makes sense?

My guess is: Both strands should have the same genomic information for the most part, right? I guess one strand might be sequenced better than the other strand in some locations therefore removing strand information gives information about both strands to sort of the sum together the information that the either may lack?

Thank you in advance!

Pratik


GenomicRanges


strand


GRanges


intergenic


genic

• 17 views

Read more here: Source link