Calculating Sequencing Depth Of A Genomic Region
1
0
Entering edit mode
8.5 years ago
User 1933 ▴ 340

In an whole-exome study project, I need to calculate the sequencing depth for a given set of genes. It seems DepthOfCoverage from GATK can do this for me. However, I am curious if I can get the similar information from the VCF file ? Alternatively is there any other tool I can use ?

vcf exome-sequencing • 4.2k views
ADD COMMENT
1
Entering edit mode

As far as I know, VCF file will provide you depth only at the position of the variant. It depends if you are interested in depth at variants, average depth across the gene or in depth for every position in the gene. If the last one is your case, you can use "samtools depth" over the bam file (alignment file of reads over reference).

ADD REPLY
0
Entering edit mode
8.5 years ago

GATK DepthOfCoverage will give you want you want, but it takes bam files.

Also, read this, it's a nice discussion of using bedtools as well.

I wouldn't calculate my read depth across the target just using a vcf, since that file only includes your variant calls.

ADD COMMENT

Login before adding your answer.

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