Question: Depth of Coverage for all SNP in reference genome in a BAM file
1
gravatar for jingning629
17 months ago by
jingning6290 wrote:

Hello!

I'm new to this area. Sorry in advance if I say something inaccurately.

Now I have Raw Sequencing Data for many samples, and also a reference genome. For a given sample, I want to find out the depth of coverage (the number of reads that cover a given nucleotide) for all positions of SNPs in the reference genome. Is there any tools that I can use? I found that I can know the Read Depth using GATK HaplotypeCaller but only for the nucleotide having variants (genotype 0/1 or 1/1), which is lack of information for the 0/0 genotype.

Thank you in advance!

ADD COMMENTlink modified 17 months ago by Bastien Hervé4.4k • written 17 months ago by jingning6290

Hello,

have you already mapped and aligned you reads to your reference genome? What genome is it? What do you mean by "for all positions of SNPs in the reference genome"? If you want to know the read depth on position where your sample has a variant, just do a variant calling. In the resulting vcf file there should be that information.

fin swimmer

ADD REPLYlink written 17 months ago by finswimmer12k
0
gravatar for Bastien Hervé
17 months ago by
Bastien Hervé4.4k
Limoges, CBRS, France
Bastien Hervé4.4k wrote:

Maybe not the fastiest way but you can download the dbSNP of your reference genome.

Create a bed file with chromosome names and positions from your dbSNP.

Use samtools depth with option -b (with a specific bed file) on your bam file.

ADD COMMENTlink written 17 months ago by Bastien Hervé4.4k

Thank you! I didn't find a dbSNP file with my reference genome. The reference genome is just a .fa file. Could I generate the dbSNP using it?

ADD REPLYlink written 17 months ago by jingning6290
1

You need to provide some complementary informations. What is your reference genome, which version. Where did you get it (Gencode, Ensembl...) ?

ADD REPLYlink written 17 months ago by Bastien Hervé4.4k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1093 users visited in the last hour