Question: allele frequency bcf tools
gravatar for jh
5 months ago by
jh30 wrote:

Hi, I would like to get the read counts for the ref and alt allele in a vcf file using bcf tools call.
I use this command:

bcftools mpileup -Ou -f genome.fa {input.bam} | bcftools call -m -o {output.vcf}

But I don't see the counts for the alt allele in the output, only total DP.

Any advice would be greatly appreciated.

snp variant calling • 222 views
ADD COMMENTlink modified 5 months ago by Kevin Blighe66k • written 5 months ago by jh30
gravatar for Kevin Blighe
5 months ago by
Kevin Blighe66k
Kevin Blighe66k wrote:

Hi jh,

You need to also output the AD tag. The full command:

bcftools mpileup \
  --redo-BAQ \
  --min-BQ 30 \
  --per-sample-mF \
  -f "${Ref_FASTA}" \
  Aligned_Sorted_PCRDuped_FiltMAPQ.bam |\
    bcftools call \
      --multiallelic-caller \
      --variants-only \
      -Ob > Aligned_Sorted_PCRDuped_FiltMAPQ.bcf ;


ADD COMMENTlink written 5 months ago by Kevin Blighe66k
Please log in to add an answer.


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