The current GATK version is 3.6-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Powered by Vanilla. Made with Bootstrap.
Last chance to register for the GATK workshop next week in Basel, Switzerland! http://www.sib.swiss/training/upcoming-training-events/training/gatk-workshop-lecture

Empty output when using VariantsToVCF with a simple BED file

kevyinkevyin Posts: 13Member
edited September 2013 in Ask the GATK team

java -Xmx2g -jar /share/ClusterShare/software/contrib/gi/gatk/2.6-4-g3e5ff60/GenomeAnalysisTK-2.6-4-g3e5ff60/GenomeAnalysisTK.jar -R /share/ClusterShare/biodata/contrib/gi/gatk-resource-bundle/2.3/hg19/ucsc.hg19.fasta -T VariantsToVCF -o bedtovcf.vcf --variant:BED single.bed --dbsnp /share/ClusterShare/biodata/contrib/gi/gatk-resource-bundle/2.3/hg19/dbsnp_137.hg19.vcf -L chr1

single.bed

chr1 10153 10154 1

The output is empty

Tagged:

Best Answers

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 10,469Administrator, Dev admin
    Answer ✓

    You could probably write a script to generate dummy VCF records from the BED intervals; you'd need to include a function to retrieve the ref sequence.

    Or you could make your life much, much easier and just call variants on your reference using the BED file as intervals list (see -L argument) in EMIT_ALL_SITES mode of the UnifiedGenotyper. That will produce a valid VCF that you can then edit to fit your needs.

    Geraldine Van der Auwera, PhD

Answers

  • kevyinkevyin Posts: 13Member

    @Geraldine_VdAuwera
    Hi there, thx for the reply.
    Is there any format I can use to make a set of custom locations and create a "skeleton" vcf from it?
    Cheers

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 10,469Administrator, Dev admin
    Answer ✓

    You could probably write a script to generate dummy VCF records from the BED intervals; you'd need to include a function to retrieve the ref sequence.

    Or you could make your life much, much easier and just call variants on your reference using the BED file as intervals list (see -L argument) in EMIT_ALL_SITES mode of the UnifiedGenotyper. That will produce a valid VCF that you can then edit to fit your needs.

    Geraldine Van der Auwera, PhD

Sign In or Register to comment.