|
From: Petr D. <pd...@sa...> - 2017-09-08 19:52:14
|
Hi Janet, that is curious. Any chance you could share a small test case to reproduce the behavior? Best wishes, Petr On Fri, 2017-09-08 at 13:30 -0400, Janet Cady wrote: > I generated a VCF file using the command: > > samtools mpileup -f ref.fa -uv sample.bam -t > INFO/AD,INFO/ADF,INFO/ADR | bcftools call -m -A -o sample.vcf > > > When I run vcf-validator on the resulting file I get a message like > this for all positions with no ALT alleles: > > INFO field at chrM:16567 .. INFO tag [AD=154] expected > different number of values (expected -1, found 1),INFO tag [ADF=151] > expected different number of values (expected -1, found 1),INFO tag > [ADR=3] expected different number of values (expected -1, found 1) > > Positions with one or more ALT alleles don't get have any messages. > > The meta-information lines for these fields say: > ##INFO=<ID=AD,Number=R,Type=Integer,Description="Total allelic > depths"> > ##INFO=<ID=ADF,Number=R,Type=Integer,Description="Total allelic > depths on the forward strand"> > ##INFO=<ID=ADR,Number=R,Type=Integer,Description="Total allelic > depths on the reverse strand"> > > Based on the VCF spec (v4.2), Number=R means that these fields should > have a value for each possible allele including the reference. > Doesn't this mean it should have 1 value? Why would it expect -1 > values? > > > > ------------------------------------------------------------------- > ----------- > Check out the vibrant tech community on one of the world's most > engaging tech sites, Slashdot.org! http://sdm.link/slashdot > _______________________________________________ > Vcftools-help mailing list > Vcf...@li... > https://lists.sourceforge.net/lists/listinfo/vcftools-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. |