SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Pipeline to find somatic mutations david.tamborero Bioinformatics 6 08-09-2013 03:05 AM
Find large indels from PacBio reads? metheuse Pacific Biosciences 1 07-01-2013 03:08 PM
Can VarScan be used to find out Somatic Mutations in RNA seq Data ? swapnil2188 RNA Sequencing 1 12-10-2012 08:13 AM
find all snps/indels prbndr Bioinformatics 2 09-20-2011 11:43 AM
Using samtools/tablet to find indels agc Bioinformatics 8 07-01-2010 08:48 AM

Reply
 
Thread Tools
Old 08-08-2013, 06:45 AM   #1
ymc
Senior Member
 
Location: Hong Kong

Join Date: Mar 2010
Posts: 498
Default How do I find medium length somatic indels?

I used IndelGenotyperV2 to find somatic indels. But the longest I find is 4bp. My reads are single 75bp exome reads. I am interested in indels around 20bp. Is that possible with my reads?

I ran the mapping step with bwa-0.7.5 on default setting. Is that the cause? I heard that the -o and -e flags for bwa aln supposedly have something to do with indels. If so, how do I set it up to find indels around 20bp?

Thanks a lot in advance!
ymc is offline   Reply With Quote
Old 08-09-2013, 07:42 AM   #2
KaiYe
Senior Member
 
Location: amsterdam

Join Date: Jun 2009
Posts: 133
Default

try Pindel. there is no theoretical maximum deletion size. currently maximum insertion size for insertions is read length - 20.
KaiYe is offline   Reply With Quote
Old 08-11-2013, 07:47 AM   #3
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

For single-end reads, try novoalign. You may also consider bwa-mem with (higher matching score):

bwa mem -A2 -B3 -O5 -E1

bwa-backtrack does not find a gap more than 4bp for single-end reads. Pindel would not work with single-end, as I remember (could be wrong).
lh3 is offline   Reply With Quote
Old 08-12-2013, 10:20 AM   #4
KatsenPlatz
Junior Member
 
Location: Vancouver

Join Date: May 2010
Posts: 3
Default

Strelka does a fantastic job of finding longer somatic indels: https://sites.google.com/site/strelk...variantcaller/

For exome data, set the "isSkipDepthFilters" to 1 in order to skip depth filtration.
KatsenPlatz is offline   Reply With Quote
Old 08-19-2013, 06:20 PM   #5
ymc
Senior Member
 
Location: Hong Kong

Join Date: Mar 2010
Posts: 498
Default

I tried "bwa aln -e 50" and it seems to work. Is bwa_mem better than bwa for this? How so?
ymc is offline   Reply With Quote
Old 08-20-2013, 05:51 PM   #6
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

Quote:
Originally Posted by ymc View Post
I tried "bwa aln -e 50" and it seems to work. Is bwa_mem better than bwa for this? How so?
Yes, that may work, but bwa-mem should be better, at least faster. Bwa and bwa-mem are completely different algorithms.
lh3 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 02:57 AM.


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