This is a test version of Biostars. For the public version, visit https://www.biostars.org.
vcftools merge files problem

I'm trying to merge two vcf files, however i receive an error message 

vcf-merge SNPs_AXL_Australian_m.vcf.gz SNPs_AXL_Australian_w.vcf.gz > combined

Using column name 'SS6004477' for SNPs_AXL_Australian_m.vcf.gz:SS6004477
Using column name 'SS6004478' for SNPs_AXL_Australian_w.vcf.gz:SS6004478
Wrong number of values in SS6004477/PL at 19:41725108 .. nAlleles=1, nValues=3.
Expected 1 values for diploid genotypes or 1 for haploid genotypes.

at Vcf.pm line 172, <__ANONIO__> line 2.
Vcf::throw('Vcf4_1=HASH(0x7fe63b238a48)', 'Wrong number of values in SS6004477/PL at 19:41725108 .. nAll...') called at Vcf.pm line 1767
VcfReader::parse_AGtags('Vcf4_1=HASH(0x7fe63b238a48)', 'HASH(0x7fe63b12fd30)') called at /Users/lucashenriquesviscardi/Downloads/vcftools_0.1.13/bin/vcf-merge line 464
main::merge_vcf_files('HASH(0x7fe63b82bea8)') called at /Users/lucashenriquesviscardi/Downloads/vcftools_0.1.13/bin/vcf-merge line 12

I tried also with bcftools merge option but I still stuck in here, does anyone have any idea how deal with it?

vcftools bcftools merge
• 3,459 views • link
written by viscardi.neander • viscardi.neander

"Wrong number of values in SS6004477/PL at 19:41725108 .. nAlleles=1, nValues=3." "Wrong number of values in SS6004477/PL at 19:41725108"

So you have too many values at bp 19:41725108

doing grep "19:41725108" yourvcffile > tempfile.txt

then you check what is wrong

• 0 views • link
written by Floris Brenk Floris Brenk

I check this line already, but I believe that all this information is needed to merge two different populations...

#file 1
  #CHROM    POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  SS6004472
  19    41725108    .   G   .   123.27  .   AC=0;AF=0.00;AN=2;DP=31;MQ=60.00;MQ0=0;pSC=0.378          GT:DP:GQ:PL:A:C:G:T:IR    0/0:31:93.28:0,93,1162:0,0:0,0:13,18:0,0:0

 

#file 2
#CHROM  POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  SS6004477
19  41725108    .   G   .   126.3   .   AC=0;AF=0.00;AN=2;DP=34;MQ=60.00;MQ0=0;pSC=0.378    GT:DP:GQ:PL:A:C:G:T:IR  0/0:34:96.3:0,96,1233:0,0:0,0:16,18:0,0:0
• 0 views • link
updated by Ram Ram • written by viscardi.neander • viscardi.neander

Then try to remove this line and merge again... does this works?

grep -v "19:41725108" yourvcffile > tempfile.vcf

• 0 views • link
written by Floris Brenk Floris Brenk

Still not working:

 Incorrect number of PL fields (3) at 19:41725109, cannot merge.
• 0 views • link
written by viscardi.neander • viscardi.neander

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • 1000G Tabix download: EOF marker is absent
    written by Guilherme • Guilherme
    I want to download the last release of the phased 1000Genomes (high coverage), that it is in the hg38 build but only for a set …
  • Problem with vcf-merge
    written by vncnt.anna • vncnt.anna
    Hi I am trying to merge two VCFs which I got from using platypus joint variant calling on bam files which I had to split …
  • vcf-merge is giving me a few SNPs in the merged vcf file
    written by shinken123 shinken123
    Hi I am using vcf-merge to join to vcf files, one contain the data of 56 individuals and the other of only one. One of …
  • Merge indel VCF files
    written by Randel • Randel
    <p>I am merging tons of indel VCF files using <code>vcf-merge</code>. It does not work and I found these VCF files cannot pass the <code>vcf-validation</code>. Is …
  • Problems with vcf-merge
    written by devenvyas devenvyas
    I have a list of SNPs from an old Illumina array. I've concatenated and filtered two sets of VCF files (http://cdna.eva.mpg.de/neandertal/altai/AltaiNeandertal/VCF/, http://cdna.eva.mpg.de/denisova/VCF/hg19_1000g/) from two diploid …
  • Vcf Tools Error
    written by tpeters1213 • tpeters1213
    <p>I'm working on a project that involves merging four VCF files. The first three datasets ran smoothly through vcftools' vcf-merge, but merging the final dataset …
  • Vcf-Tools : Error "Broken Vcf Header, No Column Names?'"
    written by Pierre Lindenbaum Pierre Lindenbaum
    <p>I've indexed a set of vcf with tabix.</p> <pre><code>$ gunzip -c A.vcf.bgz ##fileformat=VCFv4.1 (...) #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT A chr1 …
  • Vcf-Compare Vcftools Error
    written by Leandro Batista • Leandro Batista
    <p>Hello,</p> <p>I'm using VCFtools to compare 2 different vcf files. When I run vcf-compare with the following command:</p> <p><em>vcf-compare Software/snpEff/chr11snps.vcf.gz Software/snpEff/chr11snps_output.vcf.gz</em></p> <p>I get the following …
  • Vcftools: Vcf-Merge How To Run ?
    written by Deniz • Deniz
    <p>[perl]$ ./vcf-merge 1.flt.vcf.gz 2.flt.vcf.gz 3.flt.vcf.gz 4.flt.vcf.gz 5.flt.vcf.gz 6.flt.vcf.gz 7.flt.vcf.gz 8.flt.vcf.gz 9.flt.vcf.gz 10.flt.vcf.gz 11.flt.vcf.gz | bgzip -c &gt; merged.vcf.gz </p> <p>-bash: bgzip: command not found The …
  • Help On Merging Vcf Files By Using Vcftools
    written by Jianfengmao • Jianfengmao
    <p>Dear BioStarers,</p> <p>I am learning VCFtools by executing VCFtools commands on VCF files in Examples folder of the VCFtools installation path. </p> <p>Please help me …