SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
reads mapping to overlapping genes? beliefbio Bioinformatics 3 04-11-2014 01:18 AM
TopHat and uniquely aligned reads bgibb Bioinformatics 13 10-02-2013 10:45 PM
Problem to use Bedtools after filtering uniquely mapped reads with samtools eilosei Bioinformatics 2 12-21-2011 04:51 PM
Overlapping and non-Overlapping pair-end reads with Tophat senpeng Illumina/Solexa 4 10-16-2011 07:43 PM
non-overlapping reads in isotigs (screenshot) simoblimo 454 Pyrosequencing 2 01-30-2011 07:13 PM

Reply
 
Thread Tools
Old 02-10-2011, 12:30 PM   #1
Yilong Li
Member
 
Location: WTSI

Join Date: Dec 2010
Posts: 41
Default Dindel problem: overlapping windows and non-uniquely re-aligned reads

Hi all!

We are trying to use Dindel for generating re-aligned .bams from Illumina PE-read data. I was able to run through the Dindel pipeline, but noticed that in the output .sam files (generated following the instructions given in Dindel manual), many reads were output more than once.

When examining my data, I realized that overlapping windows were created in Dindel Stage 2. (The makeWindows.py step.) This is probably the reason why I was getting multiple reads in the concatenated .sam files.

Has anyone else encountered the same problems while attempting to get re-aligned .bam files? Do you know any solutions? Merging windows would be an obvious option, but could this result in large windows with unmanageable number of distinct haplotypes?

Any comments are welcome! 8)
Yilong
Yilong Li is offline   Reply With Quote
Old 02-10-2011, 01:59 PM   #2
csoong
Member
 
Location: Connecticut

Join Date: Jun 2009
Posts: 74
Default

just curious, are you outputing SAM to look at the indel sequences? I never became aware of that function. Thanks.
csoong is offline   Reply With Quote
Old 02-10-2011, 02:48 PM   #3
gaffa
Member
 
Location: Gothenburg/Uppsala, Sweden

Join Date: Oct 2010
Posts: 82
Default

There is a warning about this in the Dindel manual (though no solution).
gaffa is offline   Reply With Quote
Old 02-10-2011, 11:57 PM   #4
Yilong Li
Member
 
Location: WTSI

Join Date: Dec 2010
Posts: 41
Default

Thanks for the replies csoong and gaffa.

I want to output the re-aligned reads in a new .bam file to visualize them in a in-house program. Also we could use the re-aligned .bams in GATK UnifiedGenotyper. The actual indel sequences can be seen in the VCF file produced in the last step of the Dindel pipeline I think.

gaffa: yeah I've seen the warning in the manual. I am wondering if anybody been able to combine the re-aligned reads into a single .sam file.
Yilong Li is offline   Reply With Quote
Old 02-28-2011, 05:12 AM   #5
Yilong Li
Member
 
Location: WTSI

Join Date: Dec 2010
Posts: 41
Default

Bump. Please share any ideas!
Yilong Li is offline   Reply With Quote
Old 03-07-2011, 03:10 PM   #6
libiyagirl
Member
 
Location: canada

Join Date: Jan 2011
Posts: 27
Default

Hey I merged all the windows from dindel into a single bam file using samtools merge command. But I didn't check if there's any overlapping windows. Do you always have this problem?
libiyagirl is offline   Reply With Quote
Reply

Tags
dindel

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:15 PM.


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