# HG changeset patch # User devteam # Date 1444420132 14400 # Node ID 4eb3c3beb9c7d3bb4a371c330c54f9b106f0393e # Parent 6202ec8aab61febdfad88a666bcae3d51e665a87 planemo upload commit 4ec9eb1570ea116d83f5464a786af6e14fb5b57d diff -r 6202ec8aab61 -r 4eb3c3beb9c7 tophat2_wrapper.xml --- a/tophat2_wrapper.xml Tue Jul 21 16:02:24 2015 -0400 +++ b/tophat2_wrapper.xml Fri Oct 09 15:48:52 2015 -0400 @@ -1,4 +1,4 @@ - + Gapped-read mapper for RNA-seq data tophat2 --version @@ -25,7 +25,7 @@ ## tophat2 - + ## Change this to accommodate the number of threads you have available. --num-threads \${GALAXY_SLOTS:-4} @@ -35,7 +35,7 @@ #if str($params.bowtie_n) == "Yes": --bowtie-n #end if - + --read-edit-dist $params.read_edit_dist --read-realign-edit-dist $params.read_realign_edit_dist -a $params.anchor_length @@ -48,7 +48,7 @@ --segment-mismatches $params.seg_mismatches --segment-length $params.seg_length --library-type $params.library_type - + ## Indel search. #if $params.indel_search.allow_indel_search == "Yes": ## --allow-indels @@ -78,11 +78,11 @@ #else: --no-coverage-search #end if - + #if str($params.microexon_search) == "Yes": --microexon-search #end if - + #if $params.fusion_search.do_search == "Yes": --fusion-search --fusion-anchor-length $params.fusion_search.anchor_len @@ -92,13 +92,13 @@ --fusion-multipairs $params.fusion_search.multipairs --fusion-ignore-chromosomes "$params.fusion_search.ignore_chromosomes" #end if - + #if $params.bowtie2_settings.b2_settings == "Yes": #if $params.bowtie2_settings.preset.b2_preset == "Yes": --b2-$params.bowtie2_settings.preset.b2_preset_select - #end if + #end if #end if - + #end if ## Read group information. @@ -113,7 +113,7 @@ #if $singlePaired.sPaired != "single" -r $singlePaired.mate_inner_distance --mate-std-dev=$singlePaired.mate_std_dev - + #if str($singlePaired.report_discordant_pairs) == "No": --no-discordant #end if @@ -127,7 +127,7 @@ ${index_path} "$singlePaired.input1" #end if - + @@ -203,13 +203,13 @@ - + - + @@ -226,7 +226,7 @@ - + @@ -261,10 +261,10 @@ - - - - + + + + @@ -384,7 +384,7 @@ - - + @@ -512,7 +512,7 @@ - + @@ -529,13 +529,12 @@ -**Tophat Overview** +**TopHat Overview** -TopHat_ is a fast splice junction mapper for RNA-Seq reads. It aligns RNA-Seq reads to mammalian-sized genomes using the ultra high-throughput short read aligner Bowtie(2), and then analyzes the mapping results to identify splice junctions between exons. Please cite: Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, and Salzberg SL. TopHat2: accurate alignment -of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 14:R36, 2013. +TopHat_ is a fast splice junction mapper for RNA-Seq reads. It aligns RNA-Seq reads to mammalian-sized genomes using the ultra high-throughput short read aligner Bowtie(2), and then analyzes the mapping results to identify splice junctions between exons. -.. _Tophat: http://ccb.jhu.edu/software/tophat/ - +.. _TopHat: http://ccb.jhu.edu/software/tophat/ + ------ **Know what you are doing** @@ -550,13 +549,13 @@ **Input formats** -Tophat accepts files in Sanger FASTQ format. Use the FASTQ Groomer to prepare your files. +TopHat accepts files in Sanger FASTQ format. Use the FASTQ Groomer to prepare your files. ------ **Outputs** -Tophat produces two output files: +TopHat produces two output files: - junctions -- A UCSC BED_ track of junctions reported by TopHat. Each junction consists of two connected BED blocks, where each block is as long as the maximal overhang of any read spanning the junction. The score is the number of alignments spanning the junction. - accepted_hits -- A list of read alignments in BAM_ format. @@ -568,27 +567,27 @@ ------- -**Tophat settings** +**TopHat settings** -All of the options have a default value. You can change any of them. Some of the options in Tophat have been implemented here. +All of the options have a default value. You can change any of them. Some of the options in TopHat have been implemented here. ------ -**Tophat parameter list** +**TopHat parameter list** -This is a list of implemented Tophat options:: +This is a list of implemented TopHat options:: - -r This is the expected (mean) inner distance between mate pairs. For, example, for paired end runs with fragments - selected at 300bp, where each end is 50bp, you should set -r to be 200. There is no default, and this parameter + -r This is the expected (mean) inner distance between mate pairs. For, example, for paired end runs with fragments + selected at 300bp, where each end is 50bp, you should set -r to be 200. There is no default, and this parameter is required for paired end runs. --mate-std-dev INT The standard deviation for the distribution on inner distances between mate pairs. The default is 20bp. - -a/--min-anchor-length INT The "anchor length". TopHat will report junctions spanned by reads with at least this many bases on each side of the junction. Note that individual spliced - alignments may span a junction with fewer than this many bases on one side. However, every junction involved in spliced alignments is supported by at least one + -a/--min-anchor-length INT The "anchor length". TopHat will report junctions spanned by reads with at least this many bases on each side of the junction. Note that individual spliced + alignments may span a junction with fewer than this many bases on one side. However, every junction involved in spliced alignments is supported by at least one read with this many bases on each side. This must be at least 3 and the default is 8. -m/--splice-mismatches INT The maximum number of mismatches that may appear in the "anchor" region of a spliced alignment. The default is 0. -i/--min-intron-length INT The minimum intron length. TopHat will ignore donor/acceptor pairs closer than this many bases apart. The default is 70. -I/--max-intron-length INT The maximum intron length. When searching for junctions ab initio, TopHat will ignore donor/acceptor pairs farther than this many bases apart, except when such a pair is supported by a split segment alignment of a long read. The default is 500000. - -g/--max-multihits INT 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 + -g/--max-multihits INT 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. -G/--GTF [GTF 2.2 file] Supply TopHat with a list of gene model annotations. TopHat will use the exon records in this file to build a set of known splice junctions for each gene, and will attempt to align reads to these junctions even if they would not normally be covered by the initial mapping. -j/--raw-juncs [juncs file] Supply TopHat with a list of raw junctions. Junctions are specified one per line, in a tab-delimited format. Records look like: [chrom] [left] [right] [+/-], left and right are zero-based coordinates, and specify the last character of the left sequenced to be spliced to the first character of the right sequence, inclusive.