SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
using bwt to align 454 reads haonmada Bioinformatics 7 11-17-2011 01:46 AM
Align reads to contigs ojy Bioinformatics 3 07-25-2011 10:16 AM
Align multiple sequences in tabular or fasta format pchiang Bioinformatics 7 07-01-2011 11:18 PM
Looking for assembler to align very short reads andylai Bioinformatics 0 06-02-2011 09:58 PM
How to align reads to known ncRNA database skblazer Bioinformatics 8 01-07-2011 05:32 PM

Reply
 
Thread Tools
Old 09-19-2011, 11:17 AM   #1
ashwatha
Member
 
Location: India

Join Date: Jul 2011
Posts: 14
Default Reads that can align to multiple places

Hi,

is there a way to parse a BAM file and look for reads that can align to two different places on the reference genome equally well? For example, a read might be aligned to some coordinate in chr1, but maybe there is an equally good alignment elsewhere in the reference.

I believe a BAM file has provision for this kind of info to be set in one of the optional columns (XA:Z ?), but BWA doesn't seem to output this column by default. Is there a way to set an option to BWA to output this? Or is there a more obvious way to get this info?

thanks,
Ashwatha.
ashwatha is offline   Reply With Quote
Old 09-20-2011, 12:35 AM   #2
simonandrews
Simon Andrews
 
Location: Babraham Inst, Cambridge, UK

Join Date: May 2009
Posts: 871
Default

I don't know about BWA, but bowtie can be made to find these kinds of reads by using the --best and --strata options which will report only hits which match equally well. It should then be fairly simple to identify reads with more than one hit.
simonandrews is offline   Reply With Quote
Old 09-21-2011, 08:21 PM   #3
ashwatha
Member
 
Location: India

Join Date: Jul 2011
Posts: 14
Default

Thanks Simon. I will look at Bowtie. If anyone has a solution for this with BWA, I would be interested in that as well.
ashwatha 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 03:54 PM.


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