gatk
gatk copied to clipboard
I need a PON vcf
Thanks a lot
I am trying to run mutect2 in tumour-only mode for which I need a panel of normal (PON)
I have tried somatic-hg38_1000g_pon.hg38.vcf.vcf which gives
./gatk Mutect2 -R resources_broad_hg38_v0_Homo_sapiens_assembly38.fasta -I /data/Continuum/WES/testAlignmentBROADGenome/results/NG-27280_CLTSS_LTS_001A_lib506241_7636_2_MarkedDup_PicMD.bam -O 3.mt2.vcf -tumor NG-27280_CLTSS_LTS_001A_lib506241_7636_2_MarkedDup_PicMD.targeted_sequencing.sample_name --af-of-alleles-not-in-resource 2.5e-06 --germline-resource af-only-gnomad.hg38.vcf.gz -pon somatic-hg38_1000g_pon.hg38.vcf.vcf
A USER ERROR has occurred: Cannot read file:///data/somatic-hg38_1000g_pon.hg38.vcf.vcf because no suitable codecs found
I know gatk4_mutect2_4136_pon.vcf.gz locates here for which I should register in GDC but because I am a postdoctoral researcher, I can not register
Could you please help me to run mutect2 in tumour-only mode using another publicly available PON?
Thanks for any help
If you download the pon file from the GATK's buckets, change .vcf.vcf to .vcf.gz and run it again.
Did it work? I had the same problem and that change didn't do it
No did not
Fereshteh Izadi @.?anonymous&ep=bwmEmailSignature> Book time to meet with @.?anonymous&ep=bwmEmailSignature>
From: mcollodetti @.> Sent: Monday, 4 March 2024 19:51 To: broadinstitute/gatk @.> Cc: Angel Izadi @.>; Author @.> Subject: Re: [broadinstitute/gatk] I need a PON vcf (Issue #8477)
CAUTION: This email originated from outside of Swansea University. Do not click links or open attachments unless you recognise the sender and know the content is safe.
RHYBUDD: Daeth yr e-bost hwn o'r tu allan i Brifysgol Abertawe. Peidiwch â chlicio ar atodiadau neu agor atodiadau oni bai eich bod chi'n adnabod yr anfonwr a'ch bod yn gwybod bod y cynnwys yn ddiogel.
Did it work? I had the same problem and that change didn't do it
— Reply to this email directly, view it on GitHubhttps://github.com/broadinstitute/gatk/issues/8477#issuecomment-1977332718, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AUS3HJIINO56ONF7EHSUANDYWTGFLAVCNFSM6AAAAAA3SP5AQKVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTSNZXGMZTENZRHA. You are receiving this because you authored the thread.Message ID: @.***>