This is a test version of Biostars. For the public version, visit https://www.biostars.org.
GATK MUTECT2

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

0 answers

No answers yet.

Log in to answer this question.