The current GATK version is 3.7-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!

Get notifications!


You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

Did you remember to?


1. Search using the upper-right search box, e.g. using the error message.
2. Try the latest version of tools.
3. Include tool and Java versions.
4. Tell us whether you are following GATK Best Practices.
5. Include relevant details, e.g. platform, DNA- or RNA-Seq, WES (+capture kit) or WGS (PCR-free or PCR+), paired- or single-end, read length, expected average coverage, somatic data, etc.
6. For tool errors, include the error stacktrace as well as the exact command.
7. For format issues, include the result of running ValidateSamFile for BAMs or ValidateVariants for VCFs.
8. For weird results, include an illustrative example, e.g. attach IGV screenshots according to Article#5484.
9. For a seeming variant that is uncalled, include results of following Article#1235.

Did we ask for a bug report?


Then follow instructions in Article#1894.

Formatting tip!


Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ``` ) each to make a code block as demonstrated here.
Powered by Vanilla. Made with Bootstrap.
Picard 2.9.0 is now available. Download and read release notes here.
GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.

What's the difference between b37 and hg19 resources?

ihleeihlee Member Posts: 2

Hi, all. I have questions on resource bundles.
Are the 'hg19' bundle files just liftover from 'b37' bundles in UCSC-style? If so, why are there some variants in only one version and not the other? For example, the variant 'rs34872315 (on chr1)' is in b37 version of dbsnp137.excluding_sites_after_129.vcf, but not in hg19 version. At first, I thought it's because of the differences in reference genome (vcf files in the bundle are fit for the accompanying reference sequences). But the reference chromosome 1 was the same in both bundles. Can you help me to understand the difference between b37 and hg19 resource bundles?

Tagged:

Answers

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie Posts: 11,427 admin

    Hi there,

    The two versions of the reference genomes are not exactly the same. There are a few differences, for example some bases that are flipped between strands. That is why we have liftover chain files to convert between the two versions. So there may be a few variants that are filtered out in one version relative to the other. But this should affect only a tiny proportion of variants.

    Geraldine Van der Auwera, PhD

  • ihleeihlee Member Posts: 2

    Yes, I expected a few differences in variant sets between one version and the other. But then, why were the accompanying reference sequence files the same? I diffed chromosome by chromosome between two fasta files, and only differences were chromosome M and Y.

  • trptyrphetrptyrphe huileixu@gmail.comMember Posts: 1

    Hi, I couldn't find the hg19Tob37 chain file in the bundle resources. Could you please let me know its updated location? Thanks.

  • Geraldine_VdAuweraGeraldine_VdAuwera Cambridge, MAMember, Administrator, Broadie Posts: 11,427 admin

    Everything we provide should be in the resource bundle; if it's not in there, we don't have it. Did you look in the bundle on our FTP, or just in the small resources directory in the downloaded package?

    Geraldine Van der Auwera, PhD

  • mr_davemr_dave MarylandMember Posts: 4

    For those searching in 2017, the liftover chain files are not in the bundle, you can find them at
    ftp://ftp.broadinstitute.org/Liftover_Chain_Files

Sign In or Register to comment.