We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

how to use gatk to analysis multi samples at different tumor phase? GenotypeGVCFs or Mutect2

Hello, I am a fresh in gatk, I got some samples from different patient at different phase. for exmaple:
normal phase : 1 patient sample
tumor phase1 : 9 patient samples
tumor phase2 : 7 patient samples
tumor phase3 : 10 patient samples

I went to attain the specific variants at different phases, by reading the gatk best practice, I have finished joint calling to each phase (CombineGVCFs+GenotypeGVCFs) .

#samples is different phase sample
for sample in $samples ; do 
    sample_gvcfs=${sample_gvcfs}"-V $outdir/${sample}/gatk/${sample}.HC.g.vcf.gz \\"\n
time $gatk CombineGVCFs \
    -R $reference/Homo_sapiens_assembly38.fasta \
    ${sample_gvcfs} \
    -O $outdir/population/${outname}.HC.g.vcf.gz && echo "** ${outname}.HC.g.vcf.gz done ** " && \
time $gatk GenotypeGVCFs \
    -R $reference/Homo_sapiens_assembly38.fasta \
    -V $outdir/population/${outname}.HC.g.vcf.gz \
    -O $outdir/population/${outname}.HC.vcf.gz && echo "** ${outname}.HC.vcf.gz done ** "

I get two questions following:

1.After getting vcf and using hard filter, I found there is still too much varients, how can I remove more unrelated variants ? (may be I can use normal sample to filter more varient, is there any tools?)

2.Mutect2 in gatk is the better choice in this situation ?

Any help will be appreciated !

Best Answers


Sign In or Register to comment.