|
|
|
|
|||||||
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| tophat max number of junctions per read? | rcorbett | Bioinformatics | 0 | 11-23-2011 11:06 AM |
| Tophat 1.3.0 max-deletion-length option BUGFIX | andrewmagis | Bioinformatics | 3 | 06-18-2011 12:07 AM |
| Tophat -j option error | cedance | Bioinformatics | 1 | 03-25-2011 06:05 AM |
| Tophat ignoring '--max-multihits' flag? | polarise | Bioinformatics | 3 | 03-04-2011 12:28 PM |
| tophat max-multihits not heeded | brentp | Bioinformatics | 1 | 06-14-2010 10:54 AM |
 |
|
|
Thread Tools |
12-21-2010, 07:35 PM
|
#1 |
|
Junior Member
Location: MA Join Date: Dec 2010
Posts: 1
|
A question about tophat -g/--max-multihits option
If the --max-multihits option of tophat is set to 1, am I going to get the "best" alignment for the read? For example, will tophat pick the alignment for a read based on the best mapping score or the minimun number of mismatches? I am wondering how tophat works if there are multiple hits and the --max-multihits option is set to 1. Thanks.
|
|
|
|
12-22-2010, 08:02 AM
|
#2 |
|
Member
Location: salt lake city, UT Join Date: Apr 2010
Posts: 72
|
i think you need to use -g , not --max-multihits due to a bug in the arg parsing in tophat.
but if you use -g 1, then you will get the best alignment and send -m 1 to bowtie. this: http://bowtie-bio.sourceforge.net/ma...e-options-best describes how the "best" is determined. it looks at mismatches, then quality. |
|
|
|
|
12-29-2010, 05:36 AM
|
#3 | |
|
Member
Location: Stockholm, Sweden Join Date: Oct 2009
Posts: 62
|
I don't believe that is correct that setting --max-multihits to 1 gives you the BEST alignment for the read. According to the tophat manual --max-multihits 1 will trow away all reads that have more than 1 valid alignment. Hence you will get back sam files that only contain the reads that have exactly one alignment if this option is set to 1.
Quote:
|
|
|
|
|
|
| Thread Tools | |
|
|
