Mercurial > repos > artbio > sr_bowtie
changeset 0:de9cd3998571 draft
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sr_bowtie commit 19d31884000f86a5e5e07cd98ecd29fe6c9b1a7e
| author | artbio |
|---|---|
| date | Sat, 02 Sep 2017 10:52:05 -0400 |
| parents | |
| children | 5f4fbba31b6a |
| files | sRbowtie.xml test-data/297_reference.fa test-data/al.fa test-data/input.fa test-data/input.fastq test-data/output.bam test-data/output.tab test-data/output2.bam tool-data/bowtie_indices.loc.sample tool_data_table_conf.xml.sample |
| diffstat | 10 files changed, 491 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/sRbowtie.xml Sat Sep 02 10:52:05 2017 -0400 @@ -0,0 +1,283 @@ +<tool id="bowtieForSmallRNA" name="sR_bowtie" version="2.0.0"> + <description>for small RNA short reads</description> + <requirements> + <requirement type="package" version="1.1.2=py27_2">bowtie</requirement> + <requirement type="package" version="1.2">samtools</requirement> + </requirements> + <command detect_errors="exit_code"><![CDATA[ + #if $refGenomeSource.genomeSource == "history": + bowtie-build -f $refGenomeSource.ownFile genome && + ln -s -f '$refGenomeSource.ownFile' genome.fa && + #set index_path = 'genome' + #else: + #set index_path = $refGenomeSource.index.fields.path + #end if + #if $input.extension == "fasta": + #set format = "-f" + #elif $input.extension == "fastq": + #set format = "-q" + #end if + + ## set the method_prefix + #if $method == "RNA": + #set method_prefix = "-v %s -M 1 --best --strata --norc" % str($v_mismatches) + #elif $method == "unique": + #set method_prefix = "-v %s -m 1" % str($v_mismatches) + #elif $method == "multiple": + #set method_prefix = "-v %s -M 1 --best --strata" % str($v_mismatches) + #elif $method == "k_option": + #set method_prefix = "-v %s -k 1 --best" % str($v_mismatches) + #elif $method == "n_option": + #set method_prefix = "-n %s -M 1 --best" % str($v_mismatches) + #elif $method == "a_option": + #set method_prefix = "-v %s -a --best" % str($v_mismatches) + #end if + + ## set the extra_output + #if $additional_fasta == "No": + #set extra_output = "" + #elif $additional_fasta == "al": + #set extra_output = " --al %s " % str($aligned) + #elif $additional_fasta == "unal": + #set extra_output = " --un %s " % str($unaligned) + #else: + #set extra_output = " --al %s --un %s " % (str($aligned), str($unaligned)) + #end if + + #set $method_postfix = " %s %s " % ($method_prefix, $extra_output) + + ## run the bowtie alignement + #if $output_format == "tabular": + bowtie -p \${GALAXY_SLOTS:-4} $method_postfix --suppress 6,7,8 $index_path $format '$input' > $output + #elif $output_format == "sam": + bowtie -p \${GALAXY_SLOTS:-4} $method_postfix -S $index_path $format '$input' > '$output' + #elif $output_format == "bam": + bowtie -p \${GALAXY_SLOTS:-4} $method_postfix -S $index_path $format '$input' | samtools view -u - | samtools sort -@ "\${GALAXY_SLOTS:-4}" -T tmp -O bam -o $output + #end if + ##### | samtools view -uS + ]]></command> + <inputs> + <param format="fasta, fastq" help="Only with clipped, fasta or fastq read files" label="Input fasta or fastq file: reads clipped from their adapter" name="input" type="data" /> + <param help="bowtie parameters adjusted to the type of matching. RNA option match to only one strand" label="What kind of matching do you want to do?" name="method" type="select"> + <option value="RNA">Match on sense strand RNA reference index, multiple mappers randomly matched at a single position</option> + <option value="unique">Match unique mappers on DNA reference index</option> + <option selected="true" value="multiple">Match on DNA, multiple mappers randomly matched at a single position</option> + <option value="k_option">Match on DNA as fast as possible, without taking care of mapping issues (for raw annotation of reads)</option> + <option value="n_option">Match on DNA - RNAseq mode (-n bowtie option)</option> + <option value="a_option">Match and report all valid alignments</option> + </param> + <param help="specify the -v bowtie option" label="Number of mismatches allowed" name="v_mismatches" type="select"> + <option value="0">0</option> + <option selected="true" value="1">1</option> + <option value="2">2</option> + <option value="3">3</option> + </param> + <conditional name="refGenomeSource"> + <param help="Built-ins were indexed using default options" label="Will you select a reference genome from your history or use a built-in index?" name="genomeSource" type="select"> + <option value="indexed">Use a built-in index</option> + <option value="history">Use one from the history</option> + </param> + <when value="indexed"> + <param help="if your genome of interest is not listed - contact instance administrator" label="Select a DNA reference index" name="index" type="select"> + <options from_data_table="bowtie_indexes"> + + </options> + </param> + </when> + <when value="history"> + <param format="fasta" label="Select a fasta file, to serve as index reference" name="ownFile" type="data" /> + </when> + </conditional> + <param help="Note that the BAM will be viewable in trackster only if you choose a full genome referenced for Trackster usage. see the doc below" label="Select output format" name="output_format" type="select"> + <option selected="true" value="tabular">tabular</option> + <option value="sam">sam</option> + <option value="bam">bam</option> + </param> + <param help="to get aligned and unaligned reads in fasta format" label="additional fasta output" name="additional_fasta" type="select"> + <option selected="true" value="No">No</option> + <option value="al">aligned</option> + <option value="unal">unaligned</option> + <option value="al_and_unal">both aligned and unaligned</option> + </param> + </inputs> + <outputs> + <data format="tabular" label="Bowtie Output" name="output"> + <change_format> + <when format="sam" input="output_format" value="sam" /> + <when format="bam" input="output_format" value="bam" /> + </change_format> + <actions> + <conditional name="refGenomeSource.genomeSource"> + <when value="indexed"> + <action name="dbkey" type="metadata"> + <option column="1" name="bowtie_indexes" offset="0" type="from_data_table"> + <filter column="0" compare="startswith" keep="False" type="param_value" value="#" /> + <filter column="0" ref="refGenomeSource.index" type="param_value" /> + </option> + </action> + </when> + <when value="history"> + <action name="dbkey" type="metadata"> + <option name="refGenomeSource.ownFile" param_attribute="dbkey" type="from_param" /> + </action> + </when> + </conditional> + </actions> + </data> + <data format="fasta" label="Matched reads" name="aligned"> + <filter>additional_fasta == "al" or additional_fasta == "al_and_unal"</filter> + </data> + <data format="fasta" label="Unmatched reads" name="unaligned"> + <filter>additional_fasta == "unal" or additional_fasta == "al_and_unal"</filter> + </data> + </outputs> + <tests> + <test> + <param name="genomeSource" value="history" /> + <param ftype="fasta" name="ownFile" value="297_reference.fa" /> + <param name="method" value="unique" /> + <param ftype="fasta" name="input" value="input.fa" /> + <param name="v_mismatches" value="1" /> + <param name="output_format" value="bam" /> + <output file="output.bam" ftype="bam" compare="sim_size" delta="1000" name="output" /> + </test> + <test> + <param name="genomeSource" value="history" /> + <param ftype="fasta" name="ownFile" value="297_reference.fa" /> + <param name="method" value="unique" /> + <param ftype="fastq" name="input" value="input.fastq" /> + <param name="v_mismatches" value="1" /> + <param name="output_format" value="bam" /> + <output file="output2.bam" ftype="bam" compare="sim_size" delta="1000" name="output" /> + </test> + <test> + <param name="genomeSource" value="history" /> + <param ftype="fasta" name="ownFile" value="297_reference.fa" /> + <param name="method" value="multiple" /> + <param ftype="fasta" name="input" value="input.fa" /> + <param name="v_mismatches" value="1" /> + <param name="output_format" value="tabular" /> + <output file="output.tab" ftype="tabular" name="output" /> + </test> + <test> + <param name="genomeSource" value="history" /> + <param ftype="fasta" name="ownFile" value="297_reference.fa" /> + <param name="method" value="multiple" /> + <param ftype="fasta" name="input" value="input.fa" /> + <param name="v_mismatches" value="1" /> + <param name="additional_fasta" value="al" /> + <param name="output_format" value="tabular" /> + <output file="output.tab" ftype="tabular" name="output" /> + <output file="al.fa" ftype="fasta" name="aligned" /> + </test> + </tests> + <help> + +**What it does** + +Bowtie_ is a short read aligner designed to be ultrafast and memory-efficient. It is developed by Ben Langmead and Cole Trapnell. Please cite: Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biology 10:R25. + +.. _Bowtie: http://bowtie-bio.sourceforge.net/index.shtml + +A generic "Map with Bowtie for Illumina" Galaxy tool is available in the main Galaxy distribution. + +However, this Bowtie wrapper tool only takes FASTQ files as inputs. + +The sRbowtie wrapper works with short (-v bowtie mode) reads inputs, in fasta or fastq format, and proposes a simplified set of configurations suited to small RNA analysis. + +------ + +**OPTIONS** + +.. class:: infomark + +This script uses Bowtie to match reads on a reference index. + +Depending on the type of matching, different bowtie options are used: + +**Match on sense strand RNA reference index, multiple mappers randomly matched at a single position** + +Align on RNA reference, SENSE strand, randomly attributing multiple mapper to target with least mismatches: + +*-v [0,1,2,3] -M 1 --best --strata -p 12 --norc* + +**Match unique mappers on DNA reference index** + +Align ONLY unique mappers on DNA reference index + +*-v [0,1,2,3] -m 1 -p 12* + +Note that using this option with -v values other than 0 is questionnable... + +**Match on DNA, multiple mappers randomly matched at a single position** + +Align multiple mappers, randomly attributing multiple mapper to target with least mismatches, number of mismatch allowed specified by -v option: + +*-v [0,1,2,3] -M 1 --best --strata -p 12* + +**Match on DNA as fast as possible, without taking care of mapping issues (for raw annotation of reads)** + +Align with highest speed, not guaranteeing best hit for speed gain: + +*-v [0,1,2,3] -k 1 --best -p 12* + +**Match on DNA - RNAseq mode (-n bowtie option)** + +Align reads in as for RNAseq data alignment + +*-n [0,1,2,3] -M 1 --best -p 12* + +**Match and report all valid alignments** + +Align reads and report all valid alignments + +*-v [0,1,2,3] -a --best -p 12* + + + +----- + +**Input formats** + +.. class:: warningmark + +*Lists of reads, in fasta or fastq format, clipped from their adapter sequence* + +----- + +**OUTPUTS** + +If you choose tabular as the output format, you will obtain the matched reads in tabular bowtie output format (--suppress 6,7,8), having the following columns:: + + Column Description + -------- -------------------------------------------------------- + 1 FastaID fasta identifier + 2 polarity + or - depending whether the match was reported on the forward or reverse strand + 3 target name of the matched target + 4 Offset O-based coordinate of the miR on the miRBase pre-miR sequence + 5 Seq sequence of the matched Read + +If you choose SAM, you will get the output in unordered SAM format. + +.. class:: warningmark + +if you choose BAM, the output will be in sorted BAM format. +To be viewable in Trackster, several condition must be fulfilled: + +.. class:: infomark + +Reads must have been matched to a genome whose chromosome names are compatible with Trackster genome indexes + +.. class:: infomark + +the database/Build (dbkey) which is indicated for the dataset (Pencil - Database/Build field) must match a Trackster genome index. + +Please contact the Galaxy instance administrator if your genome is not referenced + +**Matched and unmatched fasta reads can be retrieved, for further analyses** + + </help> + <citations> + <citation type="doi">10.1186/gb-2009-10-3-r25</citation> + </citations> +</tool>
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/297_reference.fa Sat Sep 02 10:52:05 2017 -0400 @@ -0,0 +1,118 @@ +>FBgn0000005_297 +AGTGACGTATTTGGGTGGTCCAAACCAGCCACTTCCATTATTTCAAAGAAATCAGTAATG +CACTCTAGTAATTTTCCATAACTGTATCCCAGCTGCGCAGACTCGTTTATCTTTTGCAGC +GCAGCGTTCTTTGTAAACATCCTAAAGACCTGCCTAAGCAGATTTGACTGCCCTCTTTCA +ACGCTACCTAATCTTAAGAACCCAAGAGCGAGGCTCTCCCGAAATACAAATATTGTTCAA +ATACTGAGGCTTCTCCTCAATCCAATTTGCATTTGATTTTTAGTCTTAAGCTGAGATCCA +AAGAATAAAGTCGTGAAACTATTTCTCCTAAAAACTATTTTTTATTTCTTGGCGTTGTCC +TTAGTCAACTGACGGGACATTAGTTCGACTCATAAATAAAACAACAATTTTACTGGCGCA +GTCGGTAGGATACAAAAGTATCCGAAAAAAAAGAACCTTCGAATGGAAAATAAGTTAAAT +TTTATAGTCCTGTGCTCGAAACATCTCCCAAAATAAATTCGTGAAAACTCTTCAACTTCA +ATTATAATTCCAATTCGGTTATCCAATAATAAGTGGAAGTGAAATACGAAACAAAAATAT +TAAGTCCAAAGGCAACTAAGTTTTAAAACCAACATATAAAAATAAAAAATTAAAACAATA +TAGAATTTTAATAATACAACACAAAAATTTACAAAACAAAAAAACAAACAAGTGAAACTA +GAAAGCTTAAAAATAATAATAACATTGAATCCGAAACAAAACAAAAAAATAAAACACAAA +AGTTAAAAATTTTACAATAAAAATGTCACAACCAATTATTGCGCTGAGCGACATAAACCT +TGCCGAAGCCCGTCGGCAGCTTAAAGACATTATGCCATTCAAGGGTGATCCAGAAACCCT +TCACACCTTTATCAGCAGAGTGGATTACGTAATTTCGCTCTACCAAACAAATGATGTCCG +ACAACAGAGGATTCTACTGGGAGCCATCGAAAGGAACTTGGACGGACAAATTACACGATC +TTTGGGACTTCCGAACGTCGAAGATTGGCCTACCCTTAAAGCAAGACTCATCGCGGAATT +TAAAATTCAAACACCAAACTACAAACTTCTGGAGAACTTCAGGGAGACACCATACAGAGG +AAGCCTAAGAGCATTCTGCGAAGAAGCGGAGAGACGACGTCAATTACTAATTTCGAAACT +ACACCTGGAAGGTAACCAATCGGATTTTCTTATTTATATTCAGGGTATTAAAGAATCTAT +TAAGATACTGATAAGGAAACTACCAATACAATTATTCACTATTTTAGCCCATCACGATAT +TACAGACTTAAGATCCTTAATTACCATTGCACAAAATGAGGGAATTTATGAAGAACACAT +TAATTTTGAATTTTATGAAAAACCAGAATATCGTAATAAAAATTCAAATTCTAACCAGAA +TTCGAAAACACAAAAATTCAATACAAATGTTCAAACTCAAAATCGACCAAGTTACTCACA +ATATTCCCAACCCTTCCAACCTAATTTTAATCAATACATTCAACCATTTAGACCTAGCTA +TACACAGCAGATAACTAACAACCCACCCATGTGGCACGCACCTAATTATTTCAGACCCAA +CCAATACATAAACCCACAACCCATTATTCAAAAAAATCATTTCCAACAATATCCCAACAA +AGCCCAATTTCCCCAAACAACGCATTTTAGAGGAAATACATACCCTCGACTACAACAACC +CTCTACATATAAAAATACTAACTTCCCGATTACTAAACGACTAAGACCATCGGACAGTGA +ACAAACTAAAATGTCTATTGACGAAATTAGATTCCAAGACGCGCATGAATTCGAACAAGT +CCAACCTAATTATTACGAGCAACAGTATTTTAACCAAAATCAATACAATCCGTATCAAAA +TCATAGCTTCATTAATGAAGGGCAACAACAAGTTCAATTTGTACAAATTAATAACAAACA +AAACCAAAATAATTCTGAACTAAACGAAAATTTTCGGTTAACAGTTCCGGAAAATACGAA +TACATAAAAATAGTATACAAAGGGCGTTCATACAAATGCCTTCTAGACACAGGATCAACA +ATTAATATGATCAATGAAAATATATTTTGTCTTCCCATTCAAAATAGTAGATGTGAAGTT +TTAACATCAAATGGCCCTATTACCTTGAACGACTTGATTATGTTACCCAGAAATAGTATT +TTCAAAAAAACCGAACCATTTTATGTGCACAGATTTTCTAATAATTACGATATGCTAATT +GGCAGAAAATTGTTGAAAAATGCTCAATCAGTTATTAATTACAAAAATGATACAGTTACC +CTTTTTGATCAAACATACAAATTAATTACTTCAGAATCCGAAAGAAACCAAAATTTGTAT +ATCCAAAGGACACCAGAATCAATTGCAAGCTCAGATCAGGAATCAATAAAAAAATTAGAT +TTTTCACAGTTTCGATTAGATCACCTAAATCAGGAGGAAACTTTTAAGTTAAAAGGCTTG +TTAAATAAATTTAGAAATCTTGAATATAAGGAGGGAGAGAAATTAACATTTACAAATACA +ATTAAACACGTACTAAATACAACACATAACTCCCCAATTTATTCGAAACAATACCCACTT +GCGCAAACACACGAAATCGAAGTAGAAAACCAAGTACAGGAAATGCTGAATCAGGGATTA +ATTAGGGAAAGTAATTCTCCATACAATAGTCCTACTTGGGTCGTACCAAAGAAACCGGAT +GCTTCTGGTGCAAATAAGTACAGGGTAGTAATTGATTATAGAAAGCTAAATGAAATAACC +ATACCTGACAGATATCCAATTCCAAATATGGACGAAATTCTTGGCAAACTGGGTAAATGC +CAATATTTTACAACGATCGATCTGGCAAAGGGATTTCATCAAATAGAAATGGACGAAGAA +TCAATTTCTAAAACTGCATTCTCCACAAAAAGCGGTCATTACGAATACCTTCGAATGCCA +TTTGGCCTTAGGAATGCACCCGCTACTTTTCAAAGGTGCATGAATAATATCCTTCGACCG +TTGCTTAACAAACACTGTTTGGTGTATCTGGATGATATTATAATTTTTTCAACATCCCTT +ACAGAACATTTAAATTCAATACAATTAGTTTTTACAAAGCTTGCAGATGCAAATTTAAAA +TTGCAACTAGACAAATGTGAGTTCTTAAAAAAGGAAGCTAACTTTCTTGGTCACATAGTT +ACCCCTGATGGTATTAAACCAAATCCTATTAAAGTTAAAGCCATAGTTTCATACCCAATT +CCGACAAAAGATAAAGAGATAAGAGCTTTCCTTGGATTAACAGGTTATTATCGCAAATTT +ATTCCAAATTACGCAGACATAGCAAAACCCATGACCAGCTGCTTAAAAAAAAGGACAAAG +ATAGATACACAAAAACTTGAGTACATAGAGGCATTCGAAAAACTTAAGGCTTTGATAATT +CGTGACCCAATTTTACAATTACCTGATTTTGAAAAGAAATTTGTTTTAACCACAGATGCA +AGTAACTTGGCCCTCGGGGCTGTCCTTTCTCAAAACGGTCATCCTATATCTTTTATTAGT +AGAACACTTAACGATCACGAATTAAATTACAGTGCTATCGAAAAAGAATTACTTGCCATA +GTTTGGGCCACAAAAACTTTTCGACATTATTTACTAGGACGACAATTTCTCATTGCCAGT +GACCATCAACCTCTTAGATGGCTTCATAACTTAAAGGAACCAGGTGCTAAGTTAGAAAGA +TGGAGAGTTAGATTAAGCGAATACCAATTTAAAATAGATTATATTAAAGGGAAAGAAAAT +TCAGTTGCCGATGCATTATCAAGAATTAAAATTGAAGAAAATCATCATAGTGAAGCTACT +CAACATAGTGCAGAAGAGGACAATAGCAACCTTATTCATTTAACAGAAAAACCAATAAAT +TATTTCAAAAAACAAATAATCTTTATTAAATCCGATAAAAATAAAGTAGAGCATTCAAAA +ATATTCGGTAACTCCATTACCACAATTCAATATGACGTAATGACACTTGAAAAGGCCAAA +CAAATTTTACTCGATCACTTTATCCATAGAAACATTACCATTTATATTGAGAGCGATGTA +GATTTTGAAATCGTTCAAAGAGCACACATAGAAATTGTTAATACCACCTACACAAAAGTA +ATTCGCAGTCTTTTCCTATTAAAGAACGTTGGTTCATACGCCGAATTCAAAGAAATCATA +CTTCAATCACATGAAAAACTTTTACACCCTGGTATACAGAAAATGACAAAATTATTTAAA +GAAAATCACTTCTTTCCAAATAGCCAACTATTAATTCAGAATATAATAAACGAATGCAAC +ATATGCAATTTGGCCAAAACAGAACATAGAAACACCAAAATGCCTTTAAAAATCACACCC +AACCCGGAACATTGCCGAGAAAAATTTGTAGTAGATATTTATTCATCTGAGGGAAAACAT +TACATCAGTTGCATTGATATTTATTCTAAATTCGCTACACTTGAGCAAATTAAAACTAAG +GATTGGATAGAATGCAGAAACGCATTAATGCGCATTTTTAATCAACTAGGAAAACCCAAA +TTATTAAAGGCAGACAGAGACGGAGCTTTCTCCAGTTTAGCTTTAAAGCGATGGCTTGAA +GAAGAAGAAGTCGAATTACAGCTCAATACAGCAAAAAACGGAGTAGCAGACGTCGAAAGA +TTACACAAAACAATAAATGAAAAAATTCGTATAATCAATTCATCTGATGATGAAGAAGTA +AAATTAAGCAAGATAGAAACAATCCTCTACACATACAACCAAAAAATTAAACATGACACT +ACTGGACAGAGACCTGCTCAAATTTTCTTATACGCTGGGCATCCCATATTAGACACTCAA +AAAATTAAAGAGAAGAAAATAGAGAAAATAAATGAAGACAGACGGGAATTTAATATTGAC +ACTAATTACAGAAAAGGTCCACTACAGAAAGGCAAATTAGAAAACCCATTTAAACCAACC +AAAAATGTAGAACAGACAGACCCTGACCATTACAAAATCACTAATAGAAATAGAGTTACG +CACTACTACAAAACACAATTCAAAAAACAAAAGAAAAATAATAAACTCTCAATTTCACAG +GCACCTGGTACCCGATAACACTATTGTTTATACTGATCACAGCTGTTCATGGACAACAAA +TTCAAATTAATAATATTGACACCAACCACGGATATCTCCTTTTTTCTGATAAGCCAGTAC +AGATACCATCCTCCTTTGAACATCACTCCTTAAAAATCAATTTAACTGAAATAGACATCG +TGGTTGACTATTTTGAGCAAAGACTACGAACCGATTACCATGCACCCCAGATCAATTTTT +TATACAATAAAATAAAAAGAGAACTAGCCAGAATAACCCTGAAACATAGAAACAAACGGG +GTTTTATTAACATTGTGGGTTCAGGTTTTAAATACCTATTTGGAACACTAGATGAAAATG +ATCGAGTCGAAATACAGAAAAAACTTGAAATCAACGTCCATAACTCAGTAAAATTACATG +AACTCAACGACGCCATACGATTGATAAATGACGGAATGCAAAAAATACAGAATTATGAAA +ATAACCACACCATCATTGACAGTCTTTTGTTCGAACTAATGCAGTTTACGGAATACATAG +AAGATTTGGAAATGGCTATGCAGCTTTCCAGACTTGGACTGTTTAACCCCAAATTACTAA +ACTACGACAAACTTGAAAATGTGAACAGCCAAAACATTTTGAACATTAAAACATCCACTT +GGATTAACTACAATGATAACCAAGTATTAATCATATCCCACATACCCATTTACCTTTCAC +TAATAAGCACAATTAAAATAATTCCTTACCCAGACTCCAACGGCTATCAGCTAGATTACA +CAGACACACAATCATATTTTGAAAAAGAAAATAAAGTTTATAATACCGAAAATAAAGAAG +TAAAAAATGAATGTGTCACCAATATTATTAAACACTTAAATCCAATTTGTAATTTTAAGC +CAGTACACACGAACGAAATAATAAAATACATAGAACCAAACACAATTGTAACTTGGAACT +TAACCCAAACAATTCTTAACCAAAATTGCCAAAATTCAATTAATAAAATAAAAATAGAAG +GAAACAAAATGATAAGAGTAACGCAATGCAAAATAGAAATCAATAATATAAATTTTAGTG +AAACTCTGTTAGAACCAGAAATAGATTTGACACCACTATACACACCACTTAATATAACAA +AAATAAAAATTGTAAAACACAACGACATTATTGAGATGATTTCAGAGAACAATATTACAC +TTTACATACAAATGATCATTGTAATAATCGCACTAATTTTGTTGTACTCATATTTAAGAT +ATGTATCATTTAAACCATTTATGATGTTGTATGCAAAACTTAAAATAAGAAAAAATCAAA +ATCAAAACACACCACAACAAACAGAAATAGAAGAAATTCCATTTCCCACACTATATCCAT +CAATCCCAGCCCAAGTATAGGCTTCTCTTTAAGGGAAGGGGAGTGACGTATTTGGGTGGT +CCAAACCAGCCACTTCCATTATTTCAAAGAAATCAGTAATGCACTCTAGTAATTTTCCAT +AACTGTATCCCAGCTGCGCAGACTCGTTTATCTTTTGCAGCGCAGCGTTCTTTGTAAACA +TCCTAAAGACCTGCCTAAGCAGATTTGACTGCCCTCTTTCAACGCTACCTAATCTTAAGA +ACCCAAGAGCGAGGCTCTCCCGAAATACAAATATTGTTCAAATACTGAGGCTTCTCCTCA +ATCCAATTTGCATTTGATTTTTAGTCTTAAGCTGAGATCCAAAGAATAAAGTCGTGAAAC +TATTTCTCCTAAAAACTATTTTTTATTTCTTGGCGTTGTCCTTAGTCAACTGACGGGACA +TTAGTTCGACTCATAAATAAAACAACAATTTTACT
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/al.fa Sat Sep 02 10:52:05 2017 -0400 @@ -0,0 +1,10 @@ +>1 +CATTCTGCGAAGAAGC +>2 +GGAGAGACGACGTCAATTACT +>3 +AATTTCGAAACTACACCTGGAAGGTAACCA +>4 +ATCGGATTTTCTTATTTATATTCAGGGTATTAAAGAATCTAT +>5 +TAAGATACTGATAAGGAAACTACCAATA
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/input.fa Sat Sep 02 10:52:05 2017 -0400 @@ -0,0 +1,10 @@ +>1 +CATTCTGCGAAGAAGC +>2 +GGAGAGACGACGTCAATTACT +>3 +AATTTCGAAACTACACCTGGAAGGTAACCA +>4 +ATCGGATTTTCTTATTTATATTCAGGGTATTAAAGAATCTAT +>5 +TAAGATACTGATAAGGAAACTACCAATA
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/input.fastq Sat Sep 02 10:52:05 2017 -0400 @@ -0,0 +1,20 @@ +@HTW-1 +CATTCTGCGAAGAAGC ++ +HHHHHHHHHHHHHHHH +@HTW-2 +GGAGAGACGACGTCAATTACT ++ +HHHHHHHHHHHHHHHHHHHHH +@HTW-3 +AATTTCGAAACTACACCTGGAAGGTAACCA ++ +HHHHHHHHHHHHHHHHHHHHHHHHHHHHHH +@HTW-4 +ATCGGATTTTCTTATTTATATTCAGGGTATTAAAGAATCTAT ++ +HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH +@HTW-5 +TAAGATACTGATAAGGAAACTACCAATA ++ +HHHHHHHHHHHHHHHHHHHHHHHHHHHH
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/output.tab Sat Sep 02 10:52:05 2017 -0400 @@ -0,0 +1,5 @@ +1 + FBgn0000005_297 1151 CATTCTGCGAAGAAGC +2 + FBgn0000005_297 1167 GGAGAGACGACGTCAATTACT +3 + FBgn0000005_297 1188 AATTTCGAAACTACACCTGGAAGGTAACCA +4 + FBgn0000005_297 1218 ATCGGATTTTCTTATTTATATTCAGGGTATTAAAGAATCTAT +5 + FBgn0000005_297 1260 TAAGATACTGATAAGGAAACTACCAATA
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool-data/bowtie_indices.loc.sample Sat Sep 02 10:52:05 2017 -0400 @@ -0,0 +1,37 @@ +#This is a sample file distributed with Galaxy that enables tools +#to use a directory of Bowtie indexed sequences data files. You will +#need to create these data files and then create a bowtie_indices.loc +#file similar to this one (store it in this directory) that points to +#the directories in which those files are stored. The bowtie_indices.loc +#file has this format (longer white space characters are TAB characters): +# +#<unique_build_id> <dbkey> <display_name> <file_base_path> +# +#So, for example, if you had hg18 indexed stored in +#/depot/data2/galaxy/bowtie/hg18/, +#then the bowtie_indices.loc entry would look like this: +# +#hg18 hg18 hg18 /depot/data2/galaxy/bowtie/hg18/hg18 +# +#and your /depot/data2/galaxy/bowtie/hg18/ directory +#would contain hg18.*.ebwt files: +# +#-rw-r--r-- 1 james universe 830134 2005-09-13 10:12 hg18.1.ebwt +#-rw-r--r-- 1 james universe 527388 2005-09-13 10:12 hg18.2.ebwt +#-rw-r--r-- 1 james universe 269808 2005-09-13 10:12 hg18.3.ebwt +#...etc... +# +#Your bowtie_indices.loc file should include an entry per line for each +#index set you have stored. The "file" in the path does not actually +#exist, but it is the prefix for the actual index files. For example: +# +#hg18canon hg18 hg18 Canonical /depot/data2/galaxy/bowtie/hg18/hg18canon +#hg18full hg18 hg18 Full /depot/data2/galaxy/bowtie/hg18/hg18full +#/orig/path/hg19 hg19 hg19 /depot/data2/galaxy/bowtie/hg19/hg19 +#...etc... +# +#Note that for backwards compatibility with workflows, the unique ID of +#an entry must be the path that was in the original loc file, because that +#is the value stored in the workflow for that parameter. That is why the +#hg19 entry above looks odd. New genomes can be better-looking. +#
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool_data_table_conf.xml.sample Sat Sep 02 10:52:05 2017 -0400 @@ -0,0 +1,8 @@ +<!-- Use the file tool_data_table_conf.xml.oldlocstyle if you don't want to update your loc files as changed in revision 4550:535d276c92bc--> +<tables> + <!-- Locations of indexes in the Bowtie mapper format --> + <table name="bowtie_indexes" comment_char="#"> + <columns>value, dbkey, name, path</columns> + <file path="tool-data/bowtie_indices.loc" /> + </table> +</tables>