Using Bedtools coverage to find BAM coverage?
1
0
Entering edit mode
21 months ago
vctrm67 ▴ 20

Hello,

I am trying to use bedtools coverage to compute the depth at particular SNP sites that I list in a BED file. I am running this command:

bedtools coverage -a sample.bed -b sample.bam -mean > test

However, the output seems to list that the average coverage per site in the BED file is 0. I don't think this is true because I also ran GATK's DepthofCoverage, and it differs dramatically from what bedtools shows. Am I using bedtools incorrectly?

bedtools • 603 views
ADD COMMENT
0
Entering edit mode
21 months ago

You may want to switch the BAM to be 'a' and use -abam, and also add the -d command line parameter. Please take a look through the documentation: https://bedtools.readthedocs.io/en/latest/content/tools/coverage.html

Also check that your contig names are the same in both 'a' and 'b'.

Kevin

ADD COMMENT

Login before adding your answer.

Traffic: 1682 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