Tweet
Is there a way to control the seed length in the TopHat command?
I know in Bowtie I can use -l/--seedlen or -3/--trim3.
Thanks
Joseph
I know in Bowtie I can use -l/--seedlen or -3/--trim3.
Thanks
Joseph
-
-
i believe you'd use the --segment-length <int> option. my data is 87bp and when left up to it's own devices Tophat splits those into 3 chunks of 25, 25 and 37. I tried a few times to set the --segment-length value to the length of my reads, 87, but it crashed Tophat each time (buffer overflow). I also set --segment-length to 40 which led Tophat to split my reads into 2 parts. that didn't crash it./* Shawn Driscoll, Gene Expression Laboratory, Pfaff
Salk Institute for Biological Studies, La Jolla, CA, USA */ -
Actually, I heard from Cole that this parameter is controlled by TopHat, users cannot set it.Comment
-
users are not allowed to control the parameters sent to Bowtie because Tophat relies on a certain configuration and it has to generate a certain amount of stuff or else it's algorithms don't work. I thought about tweaking the Bowtie settings myself until I realised this from reading about Bowtie and then reading more about Tophat.
however, if you were to set --segment-length to the same length as your reads (in my case 87) then when Tophat passes the reads to Bowtie for alignment in the first stage it forces Bowtie to align the 87bp reads end-to-end. If you allow Tophat to split the reads up into smaller pieces it does so BEFORE Bowtie runs its initial alignments so it'll align reads of 25bp instead of the full 87 in my case. so in a way it is similar to the options you'd pass to Bowtie for seed length but in Tophat's case it applies to the alignment and the junction search./* Shawn Driscoll, Gene Expression Laboratory, Pfaff
Salk Institute for Biological Studies, La Jolla, CA, USA */Comment
-
What would be the difference between the two options (Tophat forcing Bowtie to do end-to-end mapping versus the default option of TopHat where the reads are split)? I am specifically interested in the effect on splice junction mapping.Comment
-
Without an answer directly from the author I'd say your best bet is to just run your data with those different options and compare the results. In fact I'm going to do this myself because I'd like to know more about how these settings alter the output as well./* Shawn Driscoll, Gene Expression Laboratory, Pfaff
Salk Institute for Biological Studies, La Jolla, CA, USA */Comment
-
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...Yesterday, 07:01 AM -
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...04-04-2024, 04:25 PM
Comment