Output AO using bcftools

Latest

Open

RNA-Seq

ChIP-Seq

SNP

Assembly

Tutorials

Tools

Jobs

Forum

Home

Welcome to Biostar! about • faq • rss

Log In

Sign Up

Question: Output AO using bcftools

0

3.1 years ago by

angrypigeon • 120

angrypigeon • 120 wrote:

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?

sequencing alignment • 950 views

ADD COMMENT • link •

modified 3.1 years ago by Jorge Amigo ♦ 11k • written 3.1 years ago by angrypigeon • 120

1

3.1 years ago by

Jorge Amigo ♦ 11k

Santiago de Compostela, Spain

Jorge Amigo ♦ 11k wrote:

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.

ADD COMMENT • link written 3.1 years ago by Jorge Amigo ♦ 11k

Please log in to add an answer.

Similar posts • Search »

bcftools split and filter variants in vcf
I use bcftools norm to split and left-normalize variants (like the line in file). After norm the ...
direct call of indels/snps with bcftools?
I am wondering... With vcftools I can call directly my snps and indels from my vcf file, but wit...
Variant Allele Frequency from Ion Torrent data
I have a vcf file with SNPs called using the Torrent Variant Caller and I want to calculate their...
samtools variant calling
Excuse me : After calling variants with samtools and bcftools, i got vcf file. Could i recal...
Short Variants calling error when using bcftools mpileup
Hello, I want to call gvcf with bcftools mpileup containing variant site and reference homologou...
Variant Calling With Bcftools Call
I have been trying to call variants using bcftools call after installing it as instructed here. ...
Variant calling on bacterial population using Freebayes gave wrong allele frequency and extremely low QUAL
I am working on genomic data of 8 evolving bacterial populations. Each one of them is a communi...
how to add QualityOfDepth(QD) annotation to freebayes vcf
Is there a tool already out there that can add the QualityOfDepth(QD) annotation to a VCF file cr...
automatically create index files from output in bcftools
say i am running a command in bcftools like /share/apps/genomics/bcftools-1.9/bin/bcftools v...
Variant calling by using samtools and bcftools
Hi all I want to do variant calling by using samtools and bcftools. As the pipelines that I foun...
Samtools/bcftools force-calling using an input set of variants
I have an input set of variants in VCF format (`input.vcf.gz`), which includes indels and SNPs. M...
How to use bcftools to calculate AF INFO field from AC and AN in VCF?
Hi, I need to add into my VCFs the AF INFO field, which is not output when variant calling with...
Bcftools Call method not detecting single nucleotide deletions?
Hi, I'm new at this and trying to construct a sequence validation pipeline. I'm using this comm...
de novo variant SNP calling from barcoded single-cell data
I have a BAM file containing sorted, aligned reads from a number of single cells. Each read has a...
getting the HOM_REF genotypes using bcftools mpileup+call ?
I've been asked to find the status of some variants (= overlaping a bed file) in a '1000 genomes'...
Calling CNV using bcftools
Hi, I used VarScan2 to call variants from whole exome sequencing data. The vcf file output looks ...
variant annotation issue
Hi, I have recently performed a variant annotation using Snpeff in water buffalo. The variants w...
BCFtools consensus sequence
I am trying to incorporate variants in my ref.fa file using bcftools consensus. But the output co...

Content

Help

About
FAQ

Access

RSS
Stats
API

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by Biostar version 2.3.0

Traffic: 1531 users visited in the last hour