Question: counting the sequencing reads in 10kb regions on a genome-wide scale
gravatar for Bogdan
11 months ago by
Palo Alto, CA, USA
Bogdan700 wrote:

Dear all,

during a step of a CNV analysis on cancer genomes, I would like to be able to COUNT and DISPLAY along the CHROMOSOME AXIS the NUMBER of READS in 10KB REGIONS of GERMLINE and CANCER GENOMES.

I would like to ask you please for your suggestions about :

-- any fast function/algorithm that COUNTS the reads from a BAM file in 10kb windows ?

-- a package/function that DISPLAY the COUNTS along the chromosome axes ?

many thanks,


sequencing cnv genome • 430 views
ADD COMMENTlink modified 8 weeks ago by Biostar ♦♦ 20 • written 11 months ago by Bogdan700

Thank you gentlemen.

Possibly, could I do it also in R, by using Rsamtools or a related package to COUNT the reads in specific genome WINDOWS ?

ADD REPLYlink written 11 months ago by Bogdan700

Yes, you could. But then you wouldn't match the speed:

In our tests, we can estimate depth across 60X genomes for 30 samples in 30 seconds.

ADD REPLYlink modified 11 months ago • written 11 months ago by h.mon23k
gravatar for igor
11 months ago by
United States
igor7.3k wrote:

You should be able to do it with bedtools. First generate a BED file of 10kb regions with bedtools makewindows. Then calculate coverage for those regions with bedtools coverage.

Alternatively, you can try deepTools bamCoverage to generate a coverage track (you can specify the size of bins).

Of course, whatever package you are using for CNV calling should have an option to generate some sort of graphical output.

ADD COMMENTlink modified 11 months ago • written 11 months ago by igor7.3k
gravatar for h.mon
11 months ago by
h.mon23k wrote:

Have a look at indexcov, but its intervals are of 16kb.

ADD COMMENTlink written 11 months ago by h.mon23k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1195 users visited in the last hour