SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Galaxy Bowtie alignment stats report ericguo Bioinformatics 1 12-08-2011 11:57 AM
bowtie mapping report cascoamarillo Bioinformatics 6 09-23-2011 12:53 AM
Bowtie alignment report kevleb Bioinformatics 1 05-04-2011 12:54 AM
Tophat options to report unaligned reads and controlling Bowtie options Siva Bioinformatics 0 10-15-2010 07:38 PM
Novoalign native report - Paired reads sandhya Bioinformatics 0 09-15-2010 06:33 AM

Reply
 
Thread Tools
Old 06-29-2012, 10:16 AM   #1
floydian_slip
Junior Member
 
Location: RTP

Join Date: Jun 2012
Posts: 5
Default Not report promiscuous reads in bowtie

Hi,
I was wondering if there is a way in bowtie to stop further alignment of short reads that have already mapped to more than a certain number of locations. Say I have a read that has aligned perfectly to 20 locations and I want bowtie to not align it further and stop processing it so that it does not waste any more time in trying to align these short promiscuous reads. Is there a way to do this in bowtie: not report/align reads that have aligned to more than say 10 locations?
Or any quick tool afterwards to delete these reads from the sam file?
thanks a lot in advance, any help will be greatly appreciated.
Best
floydian_slip is offline   Reply With Quote
Old 06-30-2012, 12:57 PM   #2
dariober
Senior Member
 
Location: Cambridge, UK

Join Date: May 2010
Posts: 310
Default

Hi- Isn't the -k option what your are looking for?
From bowtie2 manual:

Quote:
-k <int>
In this mode, bowtie2 searches for at most <int> distinct, valid alignments for each read. The search terminates when it can't find more distinct valid alignments, or when it finds <int>, whichever happens first. All alignments found are reported in descending order by alignment score. The alignment score for a paired-end alignment equals the sum of the alignment scores of the individual mates. Each reported read or pair alignment beyond the first has the SAM 'secondary' bit (which equals 256) set in its FLAGS field. For reads that have more than <int> distinct, valid alignments, bowtie2 does not gaurantee that the <int> alignments reported are the best possible in terms of alignment score. -k is mutually exclusive with -M and -a, and -M is the default.

Note: Bowtie 2 is not particularly designed with large -k in mind, and when aligning reads to long, repetitive genomes large -k can be very, very slow.
All the best
Dario
dariober is offline   Reply With Quote
Old 07-02-2012, 07:47 AM   #3
floydian_slip
Junior Member
 
Location: RTP

Join Date: Jun 2012
Posts: 5
Default

Thanks dariober, that works! thanks a lot!
Best
floydian_slip is offline   Reply With Quote
Reply

Tags
align, bowtie, reads, sam bam, short read alignment

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 04:43 PM.


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