Thank you. I see what is the problem now, did not notice the -c in bcftools call command: The DPR output is currently supported only with -m, sorry.
Petr On 7 Oct 2014, at 19:39, Koh, Chushin wrote: > Hi Petr, > > Thank you for your reply. Here is the Version info 'grep' from the same VCF > file: > > ##samtoolsVersion=1.1+htslib-1.1 > ##bcftools_callVersion=1.1+htslib-1.1 > > With Regards, > ChuShin. > > ________________________________________ > From: Petr Danecek [[email protected]] > Sent: Tuesday, October 07, 2014 12:28 PM > To: Koh, Chushin > Cc: [email protected] > Subject: Re: [Samtools-help] samtools V1.1 DPR option > > Hi ChuShin, > > only two values should appear on output at biallelic sites. I believe this > was fixed some time ago > http://github.com/samtools/bcftools/commit/91826ee6 > > Are you sure the output is really from bcftools v1.1? It's easy to check - > there is a version string in the VCF header. > > Petr > > > On 7 Oct 2014, at 18:24, Koh, Chushin wrote: > >> Hi all, >> >> I need clarification on the DPR option (see an example output below). There >> are two alleles reported for each locus, however, the DPR field has 3 values >> at pos 6554 and 4 values at pos 6690. What does the 3rd and 4th value >> correspond to (e.g. what nucleotide) in these cases? Thanks! >> >> Cheers, >> ChuShin >> >> >> -- >> command: >> $ samtools mpileup -B -d 10000 -t DP,DPR -uf reference.fa <BAM_FILES>| >> bcftools call -cv -o raw.vcf >> -- >> >> -- >> output: >> seq1 6554 . G A 121.224 . >> DP=102;VDB=0.360777;SGB=17.9188;RPB=0.404703;MQB=2.10899e-12;MQSB=0.980733;BQB=0.730302;MQ0F=0.147059;AF1=0.293951;G3=0.669799,2.37878e-06,0.330198;HWE=6.51903e-05;AC1=10;DP4=43,35,21,3;MQ=42;FQ=122.706;PV4=0.0038602,0.458111,9.36477e-23,0.161162 >> GT:PL:DP:DPR 1/1:61,15,0:5:0,5,0 0/0:0,9,103:3:3,0,0 >> 0/0:0,24,211:8:8,0,0 0/0:0,30,238:10:10,0,0 0/0:0,36,255:12:12,0,0 >> 0/0:0,18,166:6:6,0,0 0/0:0,30,239:10:10,0,0 0/0:0,0,0:0:0,0,0 >> 0/0:0,21,196:7:7,0,0 1/1:99,21,0:7:0,7,0 0/0:0,24,219:8:8,0,0 >> 0/0:0,24,222:8:8,0,0 1/1:18,21,0:7:0,7,0 0/0:0,0,0:0:0,0,0 >> 1/1:28,9,0:3:0,3,0 0/0:0,18,190:6:6,0,0 0/1:22,6,0:2:0,2,0 >> seq1 6690 . C T 999 . >> DP=152;VDB=0.981285;SGB=3.7794;RPB=0.326667;MQB=0;MQSB=8.54057e-07;BQB=1;MQ0F=0.151316;AF1=0.967821;AC1=33;DP4=1,1,86,64;MQ=35;FQ=15.3744;PV4=1,0.39789,0.00644419,0.357468 >> GT:PL:DP:DPR 1/1:145,39,0:13:0,13,0,0 1/1:157,24,0:8:0,8,0,0 >> 1/1:241,39,0:13:0,13,0,0 1/1:238,33,0:12:0,11,1,0 >> 1/1:255,45,0:15:0,15,0,0 1/1:128,15,0:5:0,5,0,0 >> 1/1:191,30,0:10:0,10,0,0 0/1:0,6,63:2:2,0,0,0 >> 1/1:163,27,0:9:0,9,0,0 1/1:161,30,0:10:0,10,0,0 >> 1/1:255,54,0:18:0,18,0,0 1/1:255,48,0:16:0,16,0,0 >> 1/1:88,24,0:8:0,8,0,0 1/1:51,9,0:3:0,3,0,0 1/1:71,9,0:3:0,3,0,0 >> 1/1:98,9,0:3:0,3,0,0 1/1:39,12,0:4:0,4,0,0 >> -- >> ------------------------------------------------------------------------------ >> Meet PCI DSS 3.0 Compliance Requirements with EventLog Analyzer >> Achieve PCI DSS 3.0 Compliant Status with Out-of-the-box PCI DSS Reports >> Are you Audit-Ready for PCI DSS 3.0 Compliance? Download White paper >> Comply to PCI DSS 3.0 Requirement 10 and 11.5 with EventLog Analyzer >> http://pubads.g.doubleclick.net/gampad/clk?id=154622311&iu=/4140/ostg.clktrk >> _______________________________________________ >> Samtools-help mailing list >> [email protected] >> https://lists.sourceforge.net/lists/listinfo/samtools-help > > > > -- > The Wellcome Trust Sanger Institute is operated by Genome Research > Limited, a charity registered in England with number 1021457 and a > company registered in England with number 2742969, whose registered > office is 215 Euston Road, London, NW1 2BE. -- The Wellcome Trust Sanger Institute is operated by Genome Research Limited, a charity registered in England with number 1021457 and a company registered in England with number 2742969, whose registered office is 215 Euston Road, London, NW1 2BE. ------------------------------------------------------------------------------ Meet PCI DSS 3.0 Compliance Requirements with EventLog Analyzer Achieve PCI DSS 3.0 Compliant Status with Out-of-the-box PCI DSS Reports Are you Audit-Ready for PCI DSS 3.0 Compliance? Download White paper Comply to PCI DSS 3.0 Requirement 10 and 11.5 with EventLog Analyzer http://pubads.g.doubleclick.net/gampad/clk?id=154622311&iu=/4140/ostg.clktrk _______________________________________________ Samtools-help mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/samtools-help
