allele frequency bcf tools
1
0
Entering edit mode
17 months ago
jh ▴ 30

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 • 897 views
ADD COMMENT
1
Entering edit mode
17 months ago

Hi jh,

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

bcftools mpileup \
  --redo-BAQ \
  --min-BQ 30 \
  --per-sample-mF \
  --annotate FORMAT/AD,FORMAT/ADF,FORMAT/ADR,FORMAT/DP,FORMAT/SP,INFO/AD,INFO/ADF,INFO/ADR \
  -f "${Ref_FASTA}" \
  Aligned_Sorted_PCRDuped_FiltMAPQ.bam |\
    bcftools call \
      --multiallelic-caller \
      --variants-only \
      -Ob > Aligned_Sorted_PCRDuped_FiltMAPQ.bcf ;

Kevin

ADD COMMENT

Login before adding your answer.

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