Question: Output all allele counts at each position using bcftools
gravatar for DVA
4.8 years ago by
United States
DVA550 wrote:


I have a bam file and would like to know all different allele counts at every position - not only the positions got called for variants, and for each position, I want the counts for all alleles shown up. By reading some previous posts: Vcf - Showing Data For Non-Variant Positions As Well As Snps and Is There A Samtools/Bcftools Setting To Call Variants No Matter How Infrequent? , I figured my code as the following:

samtools mpileup -q 1 -uf $REF $SAMPLE.sort.rudup.bam -t DP,DV,DPR,DP4,SP > $

but then when I try to get a vcf (bcftools view -c $ > $, it gave me error: incorrect number of fields (6 != 5) at 0:0

I figured I might miss something important here as I'm new to samtools and bcftools. I would appreciate any help. Thank you in advance!

samtools bcftools • 2.7k views
ADD COMMENTlink modified 4.8 years ago • written 4.8 years ago by DVA550
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: 1131 users visited in the last hour