GATK MUTECT2
0
0
Entering edit mode
2.9 years ago

I am trying to make PON for WES samples (tumor-blood matched sample) by following GATK recommendations, there is another method with 3-step procedures:

step 1. Run Mutect2 in tumor-only mode for each normal sample:

gatk Mutect2 -R reference.fasta -I normal1.bam -max-mnp-distance 0 -O normal1.vcf.gz

-> I can get vcf.gz files

step 2. Create a GenomicsDB from the normal Mutect2 calls

gatk GenomicsDBImport -R reference.fasta -L intervals.interval_list --genomicsdb-workspace-path pon_db -V normal1.vcf.gz -V normal2.vcf.gz -V normal3.vcf.gz -V ...

-> from step 2, which is the output file???, is there any default interval_list files?

step 3. Combine the normal calls using CreateSomaticPanelOfNormals:

gatk CreateSomaticPanelOfNormals -R reference.fasta --germline-resource af-only-gnomad.vcf.gz -V gendb://pon_db -O pon.vcf.gz

-> from step 3, where can I get --germline-resource ????, and -V gendb://pon_db?? can I use it at it is?

I have no idea about this command... for each step, which input file do I need and which is the output??? Let me know...

Thank you very much in advance for any help!

gatk mutect2 • 786 views
ADD COMMENT

Login before adding your answer.

Traffic: 1356 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6