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.
Announcement
Collapse
No announcement yet.
X
-
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.
-
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.
Instructs TopHat to allow up to this many alignments to the reference for a given read, and suppresses all alignments for reads with more than this many alignments. The default is 40.
Comment
Latest Articles
Collapse
-
by seqadmin
The recent pandemic caused worldwide health, economic, and social disruptions with its reverberations still felt today. A key takeaway from this event is the need for accurate and accessible tools for detecting and tracking infectious diseases. Timely identification is essential for early intervention, managing outbreaks, and preventing their spread. This article reviews several valuable tools employed in the detection and surveillance of infectious diseases.
...-
Channel: Articles
11-27-2023, 01:15 PM -
-
by seqadmin
Microbiome research has led to the discovery of important connections to human and environmental health. Sequencing has become a core investigational tool in microbiome research, a subject that we covered during a recent webinar. Our expert speakers shared a number of advancements including improved experimental workflows, research involving transmission dynamics, and invaluable analysis resources. This article recaps their informative presentations, offering insights...-
Channel: Articles
11-09-2023, 07:02 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 12-01-2023, 09:55 AM
|
0 responses
21 views
0 likes
|
Last Post
by seqadmin
12-01-2023, 09:55 AM
|
||
Started by seqadmin, 11-30-2023, 10:48 AM
|
0 responses
20 views
0 likes
|
Last Post
by seqadmin
11-30-2023, 10:48 AM
|
||
Started by seqadmin, 11-29-2023, 08:26 AM
|
0 responses
15 views
0 likes
|
Last Post
by seqadmin
11-29-2023, 08:26 AM
|
||
Started by seqadmin, 11-29-2023, 08:12 AM
|
0 responses
18 views
0 likes
|
Last Post
by seqadmin
11-29-2023, 08:12 AM
|
Comment