IndelRealigner/RealignerTargetCreator known site bundle files

mikemike Posts: 103Member
edited November 2012 in Ask the GATK team

Hi, For both IndelRealigner/RealignerTargetCreator, there is an option for known indel sites as below:

-known /path/to/indels.vcf

However, from the bundle files collection such as from hg19, there are several vcf files:

1000G_indels_for_realignment.hg19.vcf
1000G_omni2.5.hg19.sites.vcf
1000G_omni2.5.hg19.vcf
dbsnp_132.hg19.excluding_sites_after_129.vcf
dbsnp_132.hg19.vcf
hapmap_3.3.hg19.sites.vcf
hapmap_3.3.hg19.vcf
indels_mills_devine.hg19.sites.vcf
indels_mills_devine.hg19.vcf
NA12878.HiSeq.WGS.bwa.cleaned.raw.subset.hg19.sites.vcf
NA12878.HiSeq.WGS.bwa.cleaned.raw.subset.hg19.vcf

amongst them, just based on the names, 1000G_indels_for_realignment.hg19.vcf and indels_mills_devine.hg19.sites.vcf look like the files supposed to use for IndelRealigner/RealignerTargetCreator, Could you clarify the exact files for this purpose?

Since for old version, I have used 1000G_phase1.indels.hg19.vcf and Mills_and_1000G_gold_standard.indels.hg19.sites.vcf. and I compared the new and old files, quite different now.

Thanks

Mike

Post edited by Geraldine_VdAuwera on

Best Answer

Answers

  • Geraldine_VdAuweraGeraldine_VdAuwera Posts: 6,467Administrator, GATK Developer admin
  • mikemike Posts: 103Member
    edited November 2012

    Hi, Geraldine:

    Thanks for the input! However, the article seems not updated for the new version GATK v2.0 or newer. For example, the article mentioned for realignment, we shall use:

    Mills_and_1000G_gold_standard.indels.b37.sites.vcf
    1000G_phase1.indels.b37.vcf (currently from the 1000 Genomes Phase I indel calls)
    

    which is exactly what I used for the old version I described in my original comments above. But if we look at the bundle of the new version, those files are gone or at least the names somewhat changed more or less, I copied and pasted again the files in the bundle for the new version as below:

    1000G_indels_for_realignment.hg19.vcf
    1000G_omni2.5.hg19.sites.vcf
    1000G_omni2.5.hg19.vcf
    dbsnp_132.hg19.excluding_sites_after_129.vcf
    dbsnp_132.hg19.vcf
    hapmap_3.3.hg19.sites.vcf
    hapmap_3.3.hg19.vcf
    indels_mills_devine.hg19.sites.vcf
    indels_mills_devine.hg19.vcf
    NA12878.HiSeq.WGS.bwa.cleaned.raw.subset.hg19.sites.vcf
    NA12878.HiSeq.WGS.bwa.cleaned.raw.subset.hg19.vcf
    

    I think for realignment, I shall use 1000G_indels_for_realignment.hg19.vcf, but what about indels_mills_devine.hg19.sites.vcf or indels_mills_devine.hg19.vcf, which one to use for realignment?

    Thanks again

    Mike

    Post edited by Geraldine_VdAuwera on
  • mikemike Posts: 103Member

    Thanks a lot for the great detailed info, Geraldine! Appreciated! Mike

  • mikemike Posts: 103Member

    Hi, Geraldine:

    Sorry,. I just realized that your web page is actually the new version. Our own installation has some confusion about the new and old versions, which was caused by our installation staffs. Sorry about confusion. your web page is fine on that.

    Thanks any way for the info!

    Best

    Mike

Sign In or Register to comment.