SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   plotting small RNA using window size - step size (http://seqanswers.com/forums/showthread.php?t=34202)

StephaniePi83 10-02-2013 08:55 AM

plotting small RNA using window size - step size
 
Dear all,
I'm analyzing small RNA seq data.
I have the genomic coordinates of some elements (called cluster) on the genome.
I have to make some graph, representing the density of small RNA mapping to these cluster using a windows size of 250nt and step size of 25 nt ( i have to do this because it has already been done in a paper).
I'm not sure to well understand :
does it mean that i have to created a BED file called windows.bed for 250nt windows across the element of interest and that each windows overlap by 25 nt.
For example :

Code:

chr1 0 250
chr1 225 475
....

and i create a 2nd bed file of sequence reads called reads.bed

And i can use coverageBed from Bedtools to calculate for each window in windows.bed the number of reads in reads.bed that overlap the window.

Do i well understand the notion of window size / step size ?

Thanks in advance for your help,
Stéph

and then, i can use Bedtools

Kennels 10-02-2013 05:58 PM

Quote:

Originally Posted by StephaniePi83 (Post 117733)
Dear all,
I'm analyzing small RNA seq data.
I have the genomic coordinates of some elements (called cluster) on the genome.
I have to make some graph, representing the density of small RNA mapping to these cluster using a windows size of 250nt and step size of 25 nt ( i have to do this because it has already been done in a paper).
I'm not sure to well understand :
does it mean that i have to created a BED file called windows.bed for 250nt windows across the element of interest and that each windows overlap by 25 nt.
For example :

Code:

chr1 0 250
chr1 225 475
....

and i create a 2nd bed file of sequence reads called reads.bed

And i can use coverageBed from Bedtools to calculate for each window in windows.bed the number of reads in reads.bed that overlap the window.

Do i well understand the notion of window size / step size ?

Thanks in advance for your help,
Stéph

and then, i can use Bedtools

Hi,

For example, it means you count the number of small RNA that mapped to pos 0-250, then 25-275, then 50-300 and so on, and you plot the total number of mapped smallRNA in these 250nt bins.

StephaniePi83 10-03-2013 01:00 AM

Hi Kennels,
I'm not sure to well understand, what is the x-axis of the plot ? [0,25,50, ...] ?
There is redundancy no ? I mean, small RNA that have been counted in the windows [0,250] will also be counted in the windows [25,275], no ?


All times are GMT -8. The time now is 11:42 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.