Question: large insertion/deletion detection from NGS data
2
gravatar for illinois.ks
3.7 years ago by
illinois.ks140
Korea, Republic Of
illinois.ks140 wrote:

Dear all,

I am struggling with identification of large insertion/deletion as well as variants calls. I have used the GATK for snps call s and realized that GATK tools identified up to 50bp(?) insertion/deletion as well as SNPs. 

However  I am wondering whether I could also detect the large size insertion/deletion from paired-end NGS data.

I know that insertion is much more difficult to detect .. And especially, insertions which are larger than the distance of paired ends reads could not be detected (impossible to detect)..  Then, still can I detect the large size deletion ?  

 

Is there any size limit for Indels we can detect ? In my case, my fragment size  is 52bp ~ 546bp.. And median size is 216bp.  In this case,  how long Indel can I detect??  What kind of tool I have to use? 

 

KS.

 

deletion insertion gatk indel • 5.1k views
ADD COMMENTlink modified 3.7 years ago by Brian Bushnell16k • written 3.7 years ago by illinois.ks140
3
gravatar for Dan D
3.7 years ago by
Dan D6.7k
Tennessee
Dan D6.7k wrote:

BreakDancer is the first tool that comes to mind, though its documentation is shoddy. 

I recommend running multiple tools in parallel. A good second choice in my opinion is Scalpel, which uses a microassembly approach.

ADD COMMENTlink modified 3.7 years ago • written 3.7 years ago by Dan D6.7k
3
gravatar for Zev.Kronenberg
3.7 years ago by
United States
Zev.Kronenberg11k wrote:

 

 

There are many good options for structural variation detection using paired end data.  Some tools that are maintained, easy to use, and have relatively good accuracy:

Lumpy:

https://github.com/arq5x/lumpy-sv

Delly:

https://github.com/tobiasrausch/delly

Genome Strip:

http://www.broadinstitute.org/software/genomestrip/

Wham: 

https://github.com/jewmanchue/wham (my tool) 

 

ADD COMMENTlink modified 3.7 years ago • written 3.7 years ago by Zev.Kronenberg11k
1
gravatar for Brian Bushnell
3.7 years ago by
Walnut Creek, USA
Brian Bushnell16k wrote:

BBMap, with ~216bp reads, can detect insertions up to ~60bp and deletions up to 32000bp at default settings.  With non-default settings (increasing the maxindel flag and reducing the minidentity flag) it can detect much longer deletions of hundreds of thousands of bases, and longer insertions as well, but less than read length.  Maybe insertions up to 60% of read length if you set minidentity at the minimum.  This is purely based on mapping, without realignment or reassembly - I do not recommend realignment of BBMap's output during variant calling; calling directly from consensus should be more accurate.  Though if you want insertions longer than a fraction of read length, some kind of assembly is required.

ADD COMMENTlink modified 3.7 years ago • written 3.7 years ago by Brian Bushnell16k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1134 users visited in the last hour