Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
get genes located in a genomic region adrian Bioinformatics 2 09-12-2015 09:21 PM
Tools that report read count AND read names that map to genomic features. foolishbrat Bioinformatics 1 02-05-2014 01:21 AM
Pull read sequences out of a genomic region of interest? JenBarb RNA Sequencing 1 02-13-2013 04:06 PM
List all genes belonging to a certain genomic region Marius Bioinformatics 3 02-18-2011 12:33 AM
how to get annotation form genomic region? dwb0211 Bioinformatics 7 01-13-2010 10:35 PM

Thread Tools
Old 05-17-2016, 07:15 AM   #1
Junior Member
Location: Montpellier, France

Join Date: May 2012
Posts: 6
Default How to get the median of read distribution in a genomic region?


I have reads aligned on my genome and a file of genomic regions (peaks from peakcalling).
I have reads falling into each regions according to a certain distribution and I would like to get the median of this distribution for each region, that is the genomic position corresponding to the median.

Do you know if an existing tool can do this ? Or do you have an idea how to compute this?

jajclement is offline   Reply With Quote
Old 05-17-2016, 09:52 AM   #2
Rick Westerman
Location: Purdue University, Indiana, USA

Join Date: Jun 2008
Posts: 1,104

There might be a specific tool. Aside from that, assuming you have your alignment in a BAM file then I see two methods.

1. Use the BedTools package and in particular coverageBed.

2. Use 'samtools depth' to get read depth along the reference. Put this into Excel or other spreadsheet program in order to do the median calculation.
westerman is offline   Reply With Quote
Old 05-17-2016, 01:53 PM   #3
Senior Member
Location: Bay Area

Join Date: Jun 2012
Posts: 120

The GenomicRanges package in R/Bioconductor would be well suited to this as well.
cmbetts is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 05:26 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO