Tag: VQSR

Do VQSR for HaplotypeCaller calls – Sarek

Expected Behavior Filter the calls from HaplotypeCaller with Variant Quality Score Recalibration according to GATK best practise (Tools VariantRecalibrator, ApplyRecalibration, see gatkforums.broadinstitute.org/gatk/discussion/39/variant-quality-score-recalibration-vqsr or a more recent version) Current Behavior Variant quality score recalibration currently not included. Asked Jan 26 ’18 at 08:25 malinlarsson 1 Answer: Keep in mind, that you’d…

Continue Reading Do VQSR for HaplotypeCaller calls – Sarek

gnomAD v3 link not working for download

gnomAD v3 link not working for download 1 your URL is wrong. $ wget -q -O – “https://storage.googleapis.com/gcp-public-data–gnomad/release/3.1.1/vcf/genomes/gnomad.genomes.v3.1.1.sites.chr1.vcf.bgz” | gunzip -c | head ##fileformat=VCFv4.2 ##hailversion=0.2.64-1ef70187dc78 ##FILTER=<ID=AC0,Description=”Allele count is zero after filtering out low-confidence genotypes (GQ < 20; DP < 10; and AB < 0.2 for het calls)”> ##FILTER=<ID=AS_VQSR,Description=”Failed VQSR filtering…

Continue Reading gnomAD v3 link not working for download