SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Trimming left end (5') of reads?? blindtiger454 RNA Sequencing 23 01-11-2018 12:34 AM
Linux Command for joining left & right hand reads ddaneels Bioinformatics 3 04-11-2012 07:16 AM
Is left end (_1.fq) of normally mapped reads always positive strand? bioiion Illumina/Solexa 2 02-15-2012 07:19 AM
Tophat - Left Kept Reads jkozubek Bioinformatics 2 07-27-2011 05:40 AM
What's left of 5 million reads Manu Illumina/Solexa 3 08-18-2010 03:47 AM

Reply
 
Thread Tools
Old 04-11-2012, 02:53 PM   #1
genec
Member
 
Location: San Francisco, CA, USA

Join Date: Oct 2009
Posts: 12
Default Left-aligning indels

I was under the impression that the GATK IndelRealigner was supposed to left-align indels. I know that GATK has LeftAlignIndels and LeftAlignVariants, but those both output vcf files, and I would like an adjusted left-aligned bam file.

Attached is an example where two different deletion calls should clearly be collapsed into just one deletion. This alignment is after running GATK IndelRealigner.

Are there any other tools that can left-align indels in a bam file?
Attached Images
File Type: png Screen shot 2012-04-11 at 3.46.15 PM.png (11.7 KB, 60 views)
genec is offline   Reply With Quote
Old 04-13-2012, 06:31 PM   #2
adaptivegenome
Super Moderator
 
Location: US

Join Date: Nov 2009
Posts: 437
Default

The only other local realignment tool I know of is SRMA.
adaptivegenome is offline   Reply With Quote
Old 04-13-2012, 09:14 PM   #3
nilshomer
Nils Homer
 
nilshomer's Avatar
 
Location: Boston, MA, USA

Join Date: Nov 2008
Posts: 1,285
Default

How good are you at coding? There are some functions already implemented that do this, but you would need to write the wrapper code or adapt them (they are in C/Java or something that isn't perl/python).
nilshomer is offline   Reply With Quote
Old 04-14-2012, 04:32 AM   #4
gaffa
Member
 
Location: Gothenburg/Uppsala, Sweden

Join Date: Oct 2010
Posts: 82
Default

Check out the 'bamleftalign' utility which comes with the FreeBayes SNP caller, it should do this.
gaffa is offline   Reply With Quote
Old 10-16-2012, 03:05 AM   #5
ekg
Member
 
Location: Boston, MA

Join Date: Apr 2010
Posts: 36
Default

I don't have the input data available, but this is exactly the kind of thing which bamleftalign is designed to resolve. Did it work for you?
ekg is offline   Reply With Quote
Reply

Tags
gatk, indel, realigner

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:29 AM.


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