Output AO using bcftools
I calling variants from sorted bam files using samtools mpileup and then bcftools call -cv input > output.vcf. But in the resulting variant file there is no data reported for alternate allele frequency or AO. Is there a way to get bcftools to output AO?
• 2,410 views
•
link
updated
by
Jorge Amigo
•
written
by
L. A. Liggett
•
written
by
L. A. Liggett
1 answer
you're looking for
bcftools plugin fill-AN-AC
although the plugins' installation is not documented in the manual. you'll have to make sure that you compile the plugins along with the bcftools executable, and set the environment variable BCFTOOLS_PLUGINS appropriately.
• 1 views
•
link
written
by
Jorge Amigo
Log in to answer this question.
More posts like this
-
Generating mpileup file using samtools
written by RT •
-
Retain read names in mpileup and VarScan
written by vasily.dobrovolsky •
-
merge individual runs after bcftools mpileup | bcftools call
written by Leonor •
-
alignement and Variant calling with BWA and samtools
written by simon.burgermeister •
-
Merge bam files and get VCF file
written by khatami.mahshid •
-
Set minimum variant frequency while calling variants with samtools
written by kspata
-
samtools mpileup output only contains vcf header
written by s102060125 •
-
Using SNPs and INDELs to identify genetic differences between two subgroups
written by valerie
-
direct call of indels/snps with bcftools?
written by cristina_sabiers
-
Samtools : what is threshold for the mapping quality?
written by Govardhan Anande
