This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Generating Depth of Coverage files from vcf.

Does anyone know if it's possible to obtain a depth of caverage file from a vcf format file? I'm looking to create new files with depth of coverage information in the following formats: VDR, EGFR, RYR1 or GDF.

I ned this depth of coverage information file to run Stargazer software. Also I have to use previously generated vcf files, so no chance of using BAM files to call variants again and in the process obtain the depth of coverage file.

Thanks a lot for your time and help.

ryr1 vcf egfr vdr
• 2,784 views • link
updated by sbstevenlee • sbstevenlee • written by Francisco • Francisco

in the following formats: VDR, EGFR, RYR1 or GDF.

looks like gene names to me, not "format"

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

I'm not familiar with the formats you mention, but AFAIK, VCF files cannot be used to obtain depth data, since they do not contain that type of information. Depth of coverage is usually derived from BAM files, or other formats describing read mappings.
Can you provide a snippet from your VCF file and an example of the desired output?

• 0 views • link
written by liorglic liorglic

2 answers

If you're lucky, variants might include the DP field in either the INFO field or each sample, or both.

• 0 views • link
written by Divon • Divon

Wow, I'm very late to the party here! As the developer of Stargazer, I'd like to provide some insights:

  1. VDR, EGFR, and RYR1 are indeed gene names. These genes serve as a control when normalizing read depth in the target gene (e.g. CYP2D6) to compute per-base copy number.

  2. In contrast, GDF is a file format. If you look at the original paper: "For convenience, we will refer to this output as a target GDF (GATK-DepthOfCoverage format) file".

  3. If you don't have access to BAM files, then it's not possible to perform structural variation detection.

  4. However, it doesn't necessarily mean that you can't perform PGx genotyping at all. I've developed a new PGx genotyping tool called PyPGx. You can use PyPGx to perform PGx genotyping with VCF only. For more details, please see the docs.

• 0 views • link
written by sbstevenlee • sbstevenlee

Log in to answer this question.

More posts like this