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
-
Many organizations study rare diseases, but few have a mission as impactful as Rady Children’s Institute for Genomic Medicine (RCIGM). “We are all about changing outcomes for children,” explained Dr. Stephen Kingsmore, President and CEO of the group. The institute’s initial goal was to provide rapid diagnoses for critically ill children and shorten their diagnostic odyssey, a term used to describe the long and arduous process it takes patients to obtain an accurate...12-16-2024, 07:57 AM -
Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.
Long-Read Sequencing
Long-read sequencing has seen remarkable advancements,...12-02-2024, 01:49 PM
Comment