How to find positions with higher depth relative to their surroundings
2
0
Entering edit mode
10 weeks ago
hyperdx1 • 0

Hello.

I try to find positions that have higher depth significantly than their surroundings.

Is there a good way to do this?

Currently I have depth per position data from samtools depth.

depth python samtools • 373 views
ADD COMMENT
0
Entering edit mode
10 weeks ago
d-cameron ★ 2.9k

The usual way of doing this is to use one of the many existing copy number variant caller. Doing this yourself is both a) complicated to do correctly (the stats for depth per position are complicated by the fact that adjacent position depths are highly correlated since they share reads), and b) already done by someone else.

TLDR: save yourself time and effort by using an existing CNV caller.

ADD COMMENT
0
Entering edit mode
10 weeks ago
xiaoguang ▴ 140

I saw you want to get regions that have significantly higher depth than their surroundings; you must get a p-value to show significant as well.

In my opinion, the mostly easy and correct way is using the tool called Genrich, It will calculate a p-value for each peak you got.

ADD COMMENT

Login before adding your answer.

Traffic: 1721 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6