Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Newbie Question: Calculating physical coverage from genome coverage tristanstoeber Illumina/Solexa 4 06-24-2013 11:53 AM
Calculating sequence coverage tahamasoodi Bioinformatics 2 12-08-2012 10:28 AM
calculating the coverage over a gene kjaja Bioinformatics 3 06-01-2012 07:01 AM
calculating the coverage lyz1030 Bioinformatics 3 04-03-2012 06:03 PM
what program can i use for calculating coverage? dkrtndhkd Bioinformatics 2 02-03-2012 06:39 PM

Thread Tools
Old 08-27-2015, 06:33 AM   #1
Junior Member
Location: Princeton, NJ

Join Date: May 2012
Posts: 8
Question Calculating coverage of SNPs

Hello all,

I'm having a problem calculating the depth of coverage of a set of known variants (SNPs). What I want to know is 1) what percentage of known variant sites are captured in my sequencing, 2) to what depth are those variant sites sequenced to (ie. average depth, >5 reads, >10 reads, ect).

The problem is I have hundreds of large bam files from separate experiments, and a very large set of SNPs, so I can't spare the space to output a samtools depth file and just iterate through. I was wondering if there is a tool that alone or in combination with basic unix commands will give me what I want. I have tried using coverageBed (bedtools v. 2.17.0), but I can't seem to get it to give me what I want.

I have the SNPs in a vcf file (that bedtools errors out on), and a bed file that looks like this:

M 647 648 A/C +
M 934 935 G/A +
M 941 942 T/C +
M 1098 1099 T/C +
M 1222 1223 G/A +
M 1247 1248 T/C +
M 1520 1521 G/A +
M 1692 1693 C/T +
M 1831 1832 A/G +
M 1855 1856 G/A +

But bedtools doesn't find any coverage or overlap of these SNPs (which is not true).

Thanks in advance.

Last edited by aprilw; 08-28-2015 at 05:31 AM.
aprilw is offline   Reply With Quote

bam file, bedtools, coverage, samtools, snps

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 02:11 AM.

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