This is a test version of Biostars. For the public version, visit https://www.biostars.org.
gatk legacy bundles (where to get Mills_and_1000G_gold_standard.indels.hg19.sites.vcf.gz)

I need the known indels vcf to run gatk BaseRecalibrator. So I need hg19 (not the b37) version of the known indels: Mills_and_1000G_gold_standard.indels.hg19.sites.vcf.gz

However this file is no longer available at: ftp://gsapubftp-anonymous@ftp.broadinstitute.org/bundle/hg19/Mills_and_1000G_gold_standard.indels.hg19.sites.vcf.gz

the broad institute documentation says the ftp site was disabled last year (https://gatk.broadinstitute.org/hc/en-us/articles/360035890811) so I checked the gcp bucket using the following command: gsutil ls gs://gatk-legacy-bundles/hg19/

output:

gs://gatk-legacy-bundles/hg19/
gs://gatk-legacy-bundles/hg19/ucsc.hg19.dict
gs://gatk-legacy-bundles/hg19/ucsc.hg19.fasta
gs://gatk-legacy-bundles/hg19/ucsc.hg19.fasta.fai

In other words, I cant find Mills_and_1000G_gold_standard.indels.hg19.sites.vcf.gz in the gatk-legacy-bundles bucket. Does anyone know where this file is now available for download?

indels gatk
• 8,780 views • link
updated by Ishak • Ishak • written by nkinney06 • nkinney06

Thanks but those are all hg38 not hg19. I know about liftover tools but the file I need should still be available

• 0 views • link
written by nkinney06 • nkinney06

Your best bet may be to ask GATK support. hg19 version does not seem to be readily available via a search.

• 0 views • link
written by GenoMax GenoMax

These files are for hg38, not hg19. Do you know were are the files for hg19 ?

• 0 views • link
written by Edmond • Edmond

1 answer

You can find the legacy bundle at their ftp site, but you have to be logged in as gsapubftp-anonymous (logging in as anonymous will not show the bundle).

ShujiaHuang on GitHub provides some wget commands to retrieve those files:

• 0 views • link
written by ning • ning

thanks soooooo much

• 0 views • link
written by Ishak • Ishak

Log in to answer this question.

More posts like this
  • Hey guys, I'm having a prob when using GATK4 BQSR . This dbsnp vcf file has chromosomes notated a…
    written by ANImaL • ANImaL
    anilkumar@ak-omen-laptop:~/NGStools/gatk-4.4.0.0$ gatk --java-options "-DGATK_STACKTRACE_ON_USER_EXCEPTION=true" BaseRecalibrator -I "/media/anilkumar/My Passport/CRC/fastq/C_4_mkdp.bam" -R "/media/anilkumar/My Passport/CRC/fastq/hg19.fa" --known-sites "/media/anilkumar/My Passport/CRC/fastq/dbsnp_138.b37.vcf" --known-sites "/media/anilkumar/My Passport/CRC/fastq/Mills_and_1000G_gold_standard.indels.b37.vcf" --known-sites "/media/anilkumar/My Passport/CRC/fastq/1000G_phase1.indels.b37.vcf" -O "/media/anilkumar/My Passport/CRC/fastq/C_4_bqsr.table" Using GATK …
  • Find out genome version given fasta file
    written by nhaus • nhaus
    Hello everybody, I just downloaded the primary assembly from the [GATK bundle][1] along with the dbsnp_146.vcf file (I need it for the BaseRecalibrator function). I …
  • GATK BaseRecalibrator Step
    written by Seigfried Seigfried
    Hello I have been struggling with GATK 4 and feeling an acute lack of the command line pipelines available in previous versions. I took a …
  • hg19 interval list
    written by Nicolas Rosewick Nicolas Rosewick
    Hi, In GATK's resource bundle, there is an interval list for hg38 for scatter-gather purposes. ftp://gsapubftp-anonymous@ftp.broadinstitute.org/bundle/hg38/wgs_calling_regions.hg38.interval_list Quick view of the file : chr1 10001 207666 …
  • error using bqsr and incompatible contigs
    written by bioguy24 bioguy24
    Using the below command in GATK 4.0.6.0, I am getting: A USER ERROR has occurred: Input files reference and features have incompatible contigs: No overlapping …
  • multithread support for gatk
    written by bioguy24 bioguy24
    The below command is meant to use bqsr with multithreading: gatk BaseRecalibrator --num_cpu_threads_per_data_thread=8 \ -I /home/cmccabe/Desktop/fastq/NA12878_marked_duplicates.bam \ -R /home/cmccabe/Desktop/fastq/resources/ucsc.hg19.fasta \ --known-sites /home/cmccabe/Desktop/fastq/fastq/resources/dbsnp_138.b37.vcf \ --known-sites /home/cmccabe/Desktop/fastq/resources/Mills_and_1000G_gold_standard.indels.hg19.vcf …
  • Difference in the following reference genomes, both from Broad Institute FTP?
    written by ari.nazarian • ari.nazarian
    I did some alignment and BQSR using the Homo_sapiens_assembly19.fasta that I found at ftp://ftp.broadinstitute.org/pub/seq/references. I only later discovered the GATK resource bundle at ftp://gsapubftp-anonymous@ftp.broadinstitute.org/bundle/b37/ and …
  • COSMIC vcf file compatibility for Mutect2
    written by erwan.scaon erwan.scaon
    I am running a variant calling pipeline for cancer samples. It includes Mutect2. Working on human, i started with the reference & dbsnp files contained …
  • chr20 length problem
    written by waqasnayab waqasnayab
    Hi, I have an Abi SOLiD WGS csfasta reads. After QC and filtering steps, I aligned the reads via SHRiMP2, specific for csfasta reads. I …
  • Gatk - Couldn'T Read File Bundle/1.5/Hg19/Ucsc.Hg19.Fasta Because Null
    written by muenalan • muenalan
    <p>GATK gives an obscure error message when invoked as</p> <pre><code> $ java -Xmx4g -jar GenomeAnalysisTK.jar RealignerTargetCreator \ -I picard-MarkDuplicates-JOB02_014.bam \ -R references/ftp.broadinstitute.org/bundle/1.5/hg19/ucsc.hg19.fasta \ -o test.intervals …