Mutect2 failing to produce F1R2.tar.gz and stat files on TCGA bams (tumor-normal pairs)
I am using GATK version 4.1.9.0 to do somatic variant calling on TCGA BAMs.
Mutect2 does not produce any f1r2.tar.gz or stat files (--f1r2-tar-gz argument is included as below) when I ran the following command on TCGA BAMs.
gatk --java-options "-Xmx6G" Mutect2 -R file.fasta \
-I NORMAL_DNA.bam \
-I TUMOR_DNA.bam \
-normal NORMAL \
--germline-resource $RESOURCE \
--f1r2-tar-gz file.tar.gz \
-O file.vcf.gz
I see no error messages.
How can I resolve this issue? When I ran the same script on a non-TCGA BAMs, it worked just fine. Is this a TCGA specific issue? I would appreciate any inputs/suggestions. Thank you.
• 1,447 views
•
link
0 answers
No answers yet.
Log in to answer this question.