From: Adam W. <aw...@sg...> - 2013-06-18 08:36:24
|
Hi, I recently ran samtools mpileup from a Makefile and some extra fields appeared, the QS and RPB fields. e.g: samtools mpileup -gf ref.fa seq.bam | bcftools view -cg - > all.vcf NC_000962 233 . T . 81 . DP=17;QS=1.000000,0.000000,0.000000,0.000000;AF1=0;AC1=0;DP4=10,7,0,0;MQ=60;FQ=-78 PL 0 NC_000962 234 . G . 45 . DP=16;QS=0.942609,0.057391,0.000000,0.000000;RPB=-1.518513e+00;AF1=0;AC1=0;DP4=9,6,0,1;MQ=60;FQ=-42;PV4=0.44,0.2,1,0.09 PL 0 When I run this makefile on another machine those fields aren't there. Any ideas what is making these fields appear? and also what does QS mean? I can't seem to find it in the manual. [Samtools Version: 0.1.18 (r982:295); bcftools 0.1.17-dev (r973:277)] Thanks for any help Adam |