SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
running GATK or Mutect or a somatic variant caller individual chromosomes bogdan Bioinformatics 1 10-08-2016 11:49 PM
VarScan somatic indel caller working incorrectly? ElenaN Bioinformatics 1 03-03-2014 01:31 AM
Somatic variant caller with no normal control linda_kohn Bioinformatics 0 09-11-2013 04:04 AM
variant caller that can accomodate mosaicism? ymc Bioinformatics 0 08-21-2012 04:46 PM

Reply
 
Thread Tools
Old 04-14-2016, 10:18 PM   #1
PKhongthon
Junior Member
 
Location: Thailand

Join Date: Apr 2015
Posts: 2
Default The problem with sensitivity analysis of somatic variant caller

Hello!

I'm working on somatic variant identification from cancer samples. So, I try to simulate Illumina paired-end read with known variants by using VarSim (http://bioinform.github.io/varsim/), which this tools can simulate the human genome with both germline and somatic mutations (use ART tool as a simulator). In addition, I had simulated tumor reads by vary allele frequency in the range of 0.1 - 0.5.

As this paper "Detecting somatic point mutations in cancer genome sequencing data: a comparison of mutation callers" (DOI: 10.1186/gm495) mentioned about the sensitivity of any somatic variant caller including MuTect, SomaticSniper, VarScan, Strelka, and JointSNVMix. I had followed the command of any tools that they provide in supplement to calling somatic variations. However, the sensitivity is quite low when compare with their reports. The maximum is around 0.9 at allele frequency equal to 0.5 (MuTect) while at allele frequency equal to 0.5 (MuTect) I got the sensitivity just 0.7 (my simulated reads are chr22).

So, I have no idea that why the results are very different, even I use only chr22 but I think the results shouldn't be much different.

Additional, I aligned the reads by using bwa-mem then mark duplicates with Picard-tool and realigned, recalibration with GATK
PKhongthon is offline   Reply With Quote
Reply

Tags
illumina, pipeline, sensitivity, simulated reads, somatic variant detection

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:41 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO