Fdtxjr

Stack Interview Code methods made from class Node and Smart Pointers

Why does AES have exactly 10 rounds for a 128-bit key, 12 for 192 bits and 14 for a 256-bit key size?

How do you make your own symbol when Detexify fails?

What features enable the Su-25 Frogfoot to operate with such a wide variety of fuels?

C++ check if statement can be evaluated constexpr

Non-trope happy ending?

Why should universal income be universal?

What is the English pronunciation of "pain au chocolat"?

Does an advisor owe his/her student anything? Will an advisor keep a PhD student only out of pity?

How much of a Devil Fruit must be consumed to gain the power?

Is there a RAID 0 Equivalent for RAM?

I found an audio circuit and I built it just fine, but I find it a bit too quiet. How do I amplify the output so that it is a bit louder?

Delete multiple columns using awk or sed

Change the color of a single dot in `ddot` symbol

How many arrows is an archer expected to fire by the end of the Tyranny of Dragons pair of adventures?

Is there a nicer/politer/more positive alternative for "negates"?

Creating two special characters

What fields between the rationals and the reals allow a good notion of 2D distance?

Are Captain Marvel's powers affected by Thanos breaking the Tesseract and claiming the stone?

C++ copy constructor called at return

Does the reader need to like the PoV character?

The IT department bottlenecks progress, how should I handle this?

Which was the first story featuring espers?

What does "Scientists rise up against statistical significance" mean? (Comment in Nature)

Variant calling without matched normal sample

Single-sample vs. joint genotypingAre soft-clipped bases used for variant calling in samtools + bcftools?BAM to BigWig without intermediary BedGraphCreating a new reference genome B by changing genome A with mismatches from BAM file without long readsIs it possible to use SNP heterozygosity as a proxy for Indel heterozygosity?Merge 2 VCFs from different variant callersIssues performing variant calling with GATKHuman Genome Alignment and Variant Calling BenchmarksNo variant found using GATK 4.0 HaplotypeCallercalling diploid SNVs from long reads

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?

bam snp gatk wgs vcftools

share|improve this question

edited Mar 14 at 17:20

Feresh Teh

asked Mar 14 at 14:25

Feresh TehFeresh Teh

44912

$endgroup$

add a comment | 

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?

bam snp gatk wgs vcftools

share|improve this question

edited Mar 14 at 17:20

Feresh Teh

asked Mar 14 at 14:25

Feresh TehFeresh Teh

44912

$endgroup$

add a comment | 

$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?

bam snp gatk wgs vcftools

share|improve this question

edited Mar 14 at 17:20

Feresh Teh

asked Mar 14 at 14:25

Feresh TehFeresh Teh

44912

$endgroup$

share|improve this question

edited Mar 14 at 17:20

Feresh Teh

asked Mar 14 at 14:25

Feresh TehFeresh Teh

44912

share|improve this question

edited Mar 14 at 17:20

Feresh Teh

asked Mar 14 at 14:25

Feresh TehFeresh Teh

44912

asked Mar 14 at 14:25

Feresh TehFeresh Teh

44912

asked Mar 14 at 14:25

Feresh TehFeresh Teh

44912

1 Answer
1

active

oldest

votes

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.

share|improve this answer

edited Mar 14 at 14:56

answered Mar 14 at 14:42

Devon Ryan♦Devon Ryan

14.2k21541

$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
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  software.broadinstitute.org/gatk/documentation/article?id=11682
  $endgroup$
  – Feresh Teh
  Mar 18 at 17:54
  

  
  
  
  

add a comment | 

Your Answer

StackExchange.ifUsing("editor", function ()
return StackExchange.using("mathjaxEditing", function ()
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix)
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
);
);
, "mathjax-editing");

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "676"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);



);

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);

Sign up using Google

Sign up using Facebook

Sign up using Email and Password

Post as a guest

Name

Email

Required, but never shown

StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fbioinformatics.stackexchange.com%2fquestions%2f7250%2fvariant-calling-without-matched-normal-sample%23new-answer', 'question_page');

);

Post as a guest

Name

Email

Required, but never shown

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.

share|improve this answer

edited Mar 14 at 14:56

answered Mar 14 at 14:42

Devon Ryan♦Devon Ryan

14.2k21541

$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
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  software.broadinstitute.org/gatk/documentation/article?id=11682
  $endgroup$
  – Feresh Teh
  Mar 18 at 17:54
  

  
  
  
  

add a comment | 

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.

share|improve this answer

edited Mar 14 at 14:56

answered Mar 14 at 14:42

Devon Ryan♦Devon Ryan

14.2k21541

$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
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  software.broadinstitute.org/gatk/documentation/article?id=11682
  $endgroup$
  – Feresh Teh
  Mar 18 at 17:54
  

  
  
  
  

add a comment | 

$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.

share|improve this answer

edited Mar 14 at 14:56

answered Mar 14 at 14:42

Devon Ryan♦Devon Ryan

14.2k21541

$endgroup$

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

share|improve this answer

edited Mar 14 at 14:56

answered Mar 14 at 14:42

Devon Ryan♦Devon Ryan

14.2k21541

answered Mar 14 at 14:42

Devon Ryan♦Devon Ryan

14.2k21541

answered Mar 14 at 14:42

Devon Ryan♦Devon Ryan

14.2k21541