SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
running GATK or Mutect or a somatic variant caller individual chromosomes bogdan Bioinformatics 1 10-09-2016 12:49 AM
The problem with sensitivity analysis of somatic variant caller PKhongthon Illumina/Solexa 0 04-14-2016 11:18 PM
VarScan somatic indel caller working incorrectly? ElenaN Bioinformatics 1 03-03-2014 02:31 AM
Somatic variant caller with no normal control linda_kohn Bioinformatics 0 09-11-2013 05:04 AM

Reply
 
Thread Tools
Old 11-04-2018, 09:52 PM   #1
Buckethead84
Junior Member
 
Location: New York

Join Date: Aug 2017
Posts: 7
Default Amplicon Somatic Mutation Caller

Hi everyone,

I'm relatively new to bioinformatics, and I am trying to use some deep amplicon sequencing data from the MiSeq to validate somatic mutations identified from either WES or other targeted sequencing methods.

In terms of experimental design, I have 3 amplicons tiled over each of the mosaic variants in question (range from 3% - 22% AAF) and an average sequencing depth of 10,000X/amplicon using 22 cycles of initial PCR & 8 cycles of indexing with Q5 polymerase.

I'm having trouble identifying the right tools to analyze this kind of data, since I am using unpaired samples and all of the reads will have the same start/end site, preventing me from being able to collapse the reads. Any suggestions? I've combed through pre-existing threads without much luck.
Buckethead84 is offline   Reply With Quote
Old 11-05-2018, 02:16 AM   #2
henry.wood
Member
 
Location: Leeds, UK

Join Date: Apr 2010
Posts: 63
Default

I had a problem like this a while back. If your tiled amplicons can be identified and separated form their start-end positions, you can split them into separate bam files. Varscan2 has the option of mutation calling from multiple input bam files, and works nicely on PCR products. You can't use the somatic mode, if you don't have the paired data, but it picks up the mutations fine. You just need to parse the output to see how many of your mutations are called in multiple amplicons.
Hope that helps. Follow up questions would need me to dig around in old code, so might take a bit longer.
henry.wood is offline   Reply With Quote
Reply

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 07:26 PM.


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