Output all allele counts at each position using bcftools
0
0
Entering edit mode
8.1 years ago
DVA ▴ 630

Hello,

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 > $SAMPLE.vc.bcf

but then when I try to get a vcf (bcftools view -c $SAMPLE.vc.bcf > $SAMPLE.vc.vcf), 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!

bcftools samtools • 3.6k views
ADD COMMENT

Login before adding your answer.

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