Question: Allele frequency of a gene from a .bam file.
0
gravatar for Kasthuri
3.4 years ago by
Kasthuri270
United States
Kasthuri270 wrote:

It is easy to get allele frequency for a site (mutation/snp etc.) from a bam file. How do we get the allele frequency of a gene (in some metric) from a .bam file, like dominant allele frequency of a gene, for instance? Any advice from population genetics people here? Thanks!

allele frequency • 2.3k views
ADD COMMENTlink modified 2.7 years ago by azmanr0 • written 3.4 years ago by Kasthuri270

dominant allele frequency of a gene

I'm not sure if that even makes sense. Can you elaborate on what you mean?

ADD REPLYlink written 3.4 years ago by WouterDeCoster44k

I think I mean beta allele frequency like they compute in CNVs (FREEC, for instance).

ADD REPLYlink written 3.4 years ago by Kasthuri270
0
gravatar for Pierre Lindenbaum
3.4 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum130k wrote:

call the position of your gene with samtools http://www.htslib.org/workflow/

samtools mpileup -ugf ref.fa -r "chr1:1234-2345" your.bam | bcftools call -vm -o out.vcf

your should get the variants with the AF fields

ADD COMMENTlink written 3.4 years ago by Pierre Lindenbaum130k

Thank you Pierre. Yes, I get that we can get AF from variants. But I am looking at some metric for a region, in general, like beta allele frequency they compute when calling CNVs.

ADD REPLYlink written 3.4 years ago by Kasthuri270
0
gravatar for azmanr
2.7 years ago by
azmanr0
azmanr0 wrote:

Hi,

For one chromosome, I ran the command: samtools mpileup -ugf GRCz10.fasta mutant_chr12.bam | bcftools call -vm -o chr12.mut.raw.vcf

I seem to have used the same options for mpileup and bcftools call, but I do not seem to get the AF field in my output vcf files. I do get the DP4 field nested in the info column, but is there a way to get a seperate column for AF or at least a seperate column for DP4 in the vcf file?

Thanks, Azman

ADD COMMENTlink written 2.7 years ago by azmanr0
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: 1061 users visited in the last hour