Service Notice: Due to the blizzard currently hammering the US Northeast, the Broad is shut down and the GATK forum will be mostly unattended while we hunker down and sip hot cocoa with marshmallows. Assuming the power stays on and we're able to dig ourselves out of the snow when it's all over, normal service should resume Wednesday or Thursday.

LiftoverVariants fails to update the header

nilshomernilshomer Cambridge, MAPosts: 3Member
edited May 2013 in Ask the GATK team

When running the LiftoverVariants command, the resulting VCF header has the old contig names (ex. b37 to hg19, chromosome one's name is still "1" rather than "chr1"). This only affects the header, while the variants themselves are successfully changed. I am running "2013.1-2.4.9-3-g512dc3e". The command I used is shown below: <gatk> -R <human_g1k_v37>.fasta -T LiftoverVariants -V <in>.vcf -dict <hg19>.dict -o <out>.vcf -chain <b37tohg19>.chain

Thanks for your help!

Post edited by nilshomer on

Best Answer

  • ebanksebanks Posts: 684GATK Developer mod
    Answer ✓

    Hey Nils,

    LiftoverVariants is intended to be the first of 2 processing steps for the liftover; the second step is FilterLiftedVariants. In fact, LiftoverVariants is not guaranteed to produce a valid VCF file (which is alluded to in the documentation but I will go ahead and make it much more clear). If you run FilterLiftedVariants on the output of LiftoverVariants you'll see that the contig names in the header have all been correctly replaced.

    Eric Banks, PhD -- Senior Group Leader, MPG Analysis, Broad Institute of Harvard and MIT

Answers

  • ebanksebanks Posts: 684GATK Developer mod
    Answer ✓

    Hey Nils,

    LiftoverVariants is intended to be the first of 2 processing steps for the liftover; the second step is FilterLiftedVariants. In fact, LiftoverVariants is not guaranteed to produce a valid VCF file (which is alluded to in the documentation but I will go ahead and make it much more clear). If you run FilterLiftedVariants on the output of LiftoverVariants you'll see that the contig names in the header have all been correctly replaced.

    Eric Banks, PhD -- Senior Group Leader, MPG Analysis, Broad Institute of Harvard and MIT

Sign In or Register to comment.