Tophat options to report unaligned reads and controlling Bowtie options

Siva · 10-15-2010, 08:38 PM

I run Tophat from command line on our shared Linux environment queue. The gist is this:

Code:

$ ~/tophat -o ./tophat_out_s1 -p 4 --solexa1.3-quals ./genome_bwtind/<AGPv1_genome_ind> s_1_sequence.txt

I supply the options -o (for output file) -p (no. of threads) and --solexa1.3-quals (specifying quality). I also specify the path to the bowtie index and finally the actual sequence.txt (fastq) file.

I understand that Tophat runs Bowtie twice. Can I supply an option to Tophat to only report the unaligned reads after its run? I mean something like the bowtie --un option. Similarly on these lines is it possible to control Bowtie paramters from the Tophat command?

thanks
Siva

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Bowtie Paired-end options	days369	Bioinformatics	5	12-12-2012 05:38 PM
Multiple alignment why and what for (Bowtie -k options)	pmcqt	Bioinformatics	0	11-29-2011 06:48 AM
The optimal options for TopHat in mapping 35bp SOLiD reads	holywoool	Bioinformatics	0	09-22-2011 08:10 PM
Setting Bowtie options from the Tophat command line	GiladZil	RNA Sequencing	2	08-02-2011 02:42 PM
tophat - solexa options	sohnic	Bioinformatics	4	03-25-2011 02:30 PM

10-15-2010, 08:38 PM	#1
Siva Member Location: Ames, IA Join Date: Nov 2008 Posts: 57	Tophat options to report unaligned reads and controlling Bowtie options I run Tophat from command line on our shared Linux environment queue. The gist is this: Code: $ ~/tophat -o ./tophat_out_s1 -p 4 --solexa1.3-quals ./genome_bwtind/<AGPv1_genome_ind> s_1_sequence.txt I supply the options -o (for output file) -p (no. of threads) and --solexa1.3-quals (specifying quality). I also specify the path to the bowtie index and finally the actual sequence.txt (fastq) file. I understand that Tophat runs Bowtie twice. Can I supply an option to Tophat to only report the unaligned reads after its run? I mean something like the bowtie --un option. Similarly on these lines is it possible to control Bowtie paramters from the Tophat command? thanks Siva Last edited by Siva; 10-15-2010 at 11:15 PM. Reason: Clarity