SEQanswers

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

Reply
 
Thread Tools
Old 08-27-2015, 06:33 AM   #1
aprilw
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.
April

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

Tags
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