bcftools add AD column

Hi everyone.

I would like to add to the vcf file that I obtained with samtools mpileup the field AD. I did the following command:

samtools mpileup -A -d 6000 -B -C 50 -g -f $ref  -s $input.bam | bcftools call -mv -Oz -o $output.vcf

In my vcf file I had these informations in the INFO field:

DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,0,1;MQ=39

I would like to have also AD field. Can anyone help me? Thank you!!

ad bcftools vcf

• 1,868 views • link

updated 5.2 years ago by zx8754

• written 5.2 years ago by martina.gallinaro •

1 answer

bcftools mpileup  -a 'INFO/AD' (...)

• 0 views • link

written 5.2 years ago by Pierre Lindenbaum

More posts like this

Filter 2 mapped VCF files so the "POS" are the same?
written 4.2 years ago by JBC •

I have two VCF files from different individuals. They are both mapped and have a number of non-overlapping positions. Is there a way to filter …
losing data with VCFtools
written 4.2 years ago by JBC •

I am running VCFtools. My input file looks like this: ``` NW_0148388.1 404 . A G 225 . DP=15;VDB=0.310396;SGB=-0.686358;MQSB=0.640197;MQ0F=0;AC=2;AN=2;DP4=0,0,7,7;MQ=43 GT:PL 1/1:255,42,0 NW_0148388.1 495 . T …
Bcftools how to add DP to FORMAT field (get per sample read depth for REF vs ALT alleles )
written 5.0 years ago by Linda •

I'm trying to achieve what this post was looking for https://www.biostars.org/p/77012/ Currently this is my command: bcftools mpileup -Ou --max-depth 8000 --min-MQ 30 --min-BQ 30 …
Mpileup Multiple Samples - Individual DP Values
written 6.4 years ago by Ks-seq •

Hi I used mpileup (samtools and bcftools) for multiple samples (.bam files), and found that there is only single INFO field providing values for DP,AD …
bcftools-merge for vcf files
written 7.0 years ago by evelyn

Hello, I have three `vcf` files created with `bwa` and `bcftools`. I have merged these three files using `bcftools-merge` command as following: bcftools merge file1.vcf.gz …
BCF missing sample names
written 7.0 years ago by marcus.hooker •

I called SNPs with bcftools mpileup and I had 120 input sample files, my command looked like this, bcftools mpileup -d 8000 -f Reference.fasta -Ob …
Understanding a VCF - Genotype likelihoods and genotype call
written 8.1 years ago by cetin.m

I calculated genotype likelihoods from a BAM file with the command: "samtools mpileup -g -f filename" and got the BCF file. Then I used "bcftools …
What is " ./.:. " in the GT output of VCF/BCF
written 8.8 years ago by rightmirem •

I have 70 BAM files (1 for each of 70 complete genome samples), against which I ran ```mpileup``` and ```bcftools call``` to get 70 BCF …
Calling SNPs with samtools: how to specify how many supporting observations are required?
written 9.3 years ago by Lina F

Hi all, I am interested in calling SNPs for some alignment data that I have. I am looking at solutions with `samtools mpileup` and `freebayes`. …
Good filtering practices for VCF (Bcftools 1.3)
written 10.5 years ago by umn_bist •

I have a tumor sample without matches. This sample is most likely impure without a certain ploidy, copy #. I am wondering for such samples …