2
$\begingroup$

I have WGS .bam files for 3 patients (tumour and its matched derived model namely organoid) but I don't matched normal sample. If I call variants of each patients (tumour and its matched organoid), how I can use read counts at germline heterozygous positions estimated by GATK 3.2-2 to compensate for the absence of matched normal sample? I heard people use dbsnp VCF instead of the matched normal small variant VCF. But, I don't know start from where? which GATK function does that?

I should mention, Calling SNV and indel in many tools returns vcf but I called copy number by varscan that did not return .vcf output so I am not sure what to do for CNV

Any helps please?

$\endgroup$
6
$\begingroup$

If you just want to filter out calls present in dbSNP then use:

java -jar GenomeAnalysisTK.jar \
   -T SelectVariants \
   -R reference.fasta \
   -V patient.vcf \
   --discordance dbSNP.vcf \
   -o patient.filtered.vcf

--discordance will produce calls not present in dbSNP.

$\endgroup$
  • $\begingroup$ Thanks a lot @Devon, Calling SNV and indel in many tools returns vcf but I called copy number by varscan that did not return .vcf output, is this command is applicable yet? $\endgroup$ – Feresh Teh Mar 14 at 14:49
  • 3
    $\begingroup$ It's usually best to mention things like that in your post... I presume this will still work, but I can't say I've ever tried. $\endgroup$ – Devon Ryan Mar 14 at 14:56

Your Answer

By clicking "Post Your Answer", you acknowledge that you have read our updated terms of service, privacy policy and cookie policy, and that your continued use of the website is subject to these policies.

Not the answer you're looking for? Browse other questions tagged or ask your own question.