Fdtxjr

Creating two special characters

Biological Blimps: Propulsion

What's the name of the logical fallacy where a debater extends a statement far beyond the original statement to make it true?

Can I cause damage to electrical appliances by unplugging them when they are turned on?

Why does the Sun have different day lengths, but not the gas giants?

When were female captains banned from Starfleet?

What is going on with gets(stdin) on the site coderbyte?

How much theory knowledge is actually used while playing?

What are some good ways to treat frozen vegetables such that they behave like fresh vegetables when stir frying them?

How to draw a matrix with arrows in limited space

Is it allowed to activate the ability of multiple planeswalkers in a single turn?

How to make money from a browser who sees 5 seconds into the future of any web page?

It grows, but water kills it

Giving feedback to someone without sounding prejudiced

Why Shazam when there is already Superman?

Strong empirical falsification of quantum mechanics based on vacuum energy density?

What does Apple's new App Store requirement mean

awk assign to multiple variables at once

Will number of steps recorded on FitBit/any fitness tracker add up distance in PokemonGo?

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

Is this toilet slogan correct usage of the English language?

Why should universal income be universal?

Did the UK lift the requirement for registering SIM cards?

Were Persian-Median kings illiterate?

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