Finding Somatic Mutations From Whole Exome
2
1
Entering edit mode
12.5 years ago
Angel ▴ 220

Q1) How do I call and report somatic mutations after getting the vcf files as below?

samtools mpileup -u hg19.fa foo.sorted.bam | bcftools view - > foo.sorted.bcf

bcftools view food.sorted.bcf | vcfutils varFilter -D200 > foo.d200.flt.vcf
exome mutation • 5.0k views
ADD COMMENT
0
Entering edit mode

Do you have a matched normal sample? Without a comparison, you won't be able to tell somatic from germline mutations.

ADD REPLY
0
Entering edit mode

It is not likely that all cells whose DNA was sequenced originated from the same clonal population. It is much more likely, especially for a clinical sample, that they've sequenced a mixture of cells of different origins, and hence different genomes. A matched normal is helpful, but may not be necessary.

ADD REPLY
4
Entering edit mode
12.5 years ago

Take a look at the -T pair option of bcftools. A non-comprehensive list of somatic variant callers includes varscan, somatic sniper, joint-snv-mix, and mutect (licensed from the broad). I suspect that many folks use a custom script for this application, though.

ADD COMMENT
0
Entering edit mode

Thanks so much. Sorry again for my delayed response. I will be trying some of these today.

ADD REPLY
3
Entering edit mode
12.4 years ago
Dan Koboldt ▴ 60

For exome sequencing of tumor-normal pairs, we use VarScan and SomaticSniper.

ADD COMMENT
0
Entering edit mode

Thanks very much and I will begin with varscan.

ADD REPLY
0
Entering edit mode
ADD REPLY

Login before adding your answer.

Traffic: 2664 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6