BaseRecalibrator as CountCovariates (GATK)

The CountCovariates from GATK]( has been replaced by BaseRecalibrator.

The way we can use this new BaseRecalibrator is:

java -Xmx4g -jar GenomeAnalysisTK.jar 
    -l INFO -T BaseRecalibrator 
    -I chr22.marked.realigned.fixed.bam 
    -R GRCh37.73.dna.22.fa 
    -o recalibration_report.grp 
    --knownSites dbsnp_137.hg19.vcf

This call generates the file we will use with PrintReads: recalibration_report.grp

As for CountCovariates, the new BaseRecalibrator needs the dbSNP table in VCF format. This file can be downloaded from the GATK resource bundle, available via FTP:

  • location: ``
  • username: `gsapubftp-anonymous`
  • password:

In the FTP, the dbSNP VCF files (like snp132.vcf< or dbsnp_137.hg19.vcf) are placed in /bundle/2.5/hg19, with their homologous IDX files.


Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )


Connecting to %s

%d bloggers like this: