Question: Minimum Length Of Read Segments
gravatar for Du, Jianguang
6.3 years ago by
Du, Jianguang380
Du, Jianguang380 wrote:
Dear All, I am going to run Tophat with RNA-seq dataset to observe alternative splicing events. There is a parameter for Tophat: "Minimum length of read segment". According to "implemented Tophat options", the description for "Minimum length of read segment" is "Each read is cut up into segments, each at least this long. These segments are mapped independently. The default is 25". The length of my reads is 36bps, should I change this parameter based on the length of my reads? How long should I input? Thanks. Jianguang Du
rna-seq tophat • 2.7k views
ADD COMMENTlink modified 6.3 years ago by Jennifer Hillman Jackson25k • written 6.3 years ago by Du, Jianguang380
gravatar for Jennifer Hillman Jackson
6.3 years ago by
United States
Jennifer Hillman Jackson25k wrote:
Hello Jianguang, Unless you know that your genome has exons that are shorter in length than 25 bases long, then it is not necessary to modify this parameter. --> How does TopHat find junctions? If shorter, then set the parameter to be that shortest length. Please keep in mind that lowering this value makes the job more sensitive at the cost of making it also more computationally intensive. If the job becomes too large for the public Main Galaxy instance, a cloud instance would be an alternative: Best, Jen Galaxy team -- Jennifer Jackson
ADD COMMENTlink written 6.3 years ago by Jennifer Hillman Jackson25k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 169 users visited in the last hour