# HG changeset patch # User iuc # Date 1438757443 14400 # Node ID b482293b298762f936f7c85ec7969fc63b219543 # Parent a8ac09e937f3f25089a849dea24a94fd850119ef planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/sortmerna commit 844d980a7de5f199487ca9591420a23df63c5246-dirty diff -r a8ac09e937f3 -r b482293b2987 sortmerna.py --- a/sortmerna.py Mon Aug 03 08:18:26 2015 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,72 +0,0 @@ -#!/usr/bin/env python - -""" -Runs SortMeRNA -""" - -import subprocess -import optparse -import shlex - - -def main(): - """Parse the command line, exectutes SortMeRNA and buildtrie if neeeded.""" - #TODO: Put all SortMeRNA options in the command-line parser - parser = optparse.OptionParser() - parser.add_option('--sortmerna', dest='sortmerna_cmd', help='') - parser.add_option('--buildtrie', dest='buildtrie', - default=False, action='store_true', help='') - (options, args) = parser.parse_args() - if not args: - raise Exception('Please provide at least one database') - - if options.buildtrie: - buildtrie = 'buildtrie' - for database in args: - run_buildtrie([buildtrie, '--db', database]) - - if options.sortmerna_cmd: - sortmerna = 'sortmerna' - run_sortmerna([sortmerna] + - shlex.split(options.sortmerna_cmd) + - ['-m', '262144', '-n', str(len(args)), '--db'] + - args) - - -def run_buildtrie(cmd): - """Run the BuildTrie program.""" - try: - stdout_arg = subprocess.PIPE - stderr_arg = subprocess.PIPE - child_process = subprocess.Popen(args=" ".join(cmd), shell=True, - stdin=None, stdout=stdout_arg, - stderr=stderr_arg) - stdout_str, stderr_str = child_process.communicate() - return_code = child_process.returncode - if return_code is not 0: - raise Exception(stderr_str) - - except Exception, error: - raise Exception('Error while running Buildtrie:\n' + - '\n'.join([str(error), stdout_str, stderr_str])) - - -def run_sortmerna(cmd): - """Run the SortMeRNA program.""" - try: - stdout_arg = subprocess.PIPE - stderr_arg = subprocess.PIPE - child_process = subprocess.Popen(args=" ".join(cmd), shell=True, - stdin=None, stdout=stdout_arg, - stderr=stderr_arg) - stdout_str, stderr_str = child_process.communicate() - return_code = child_process.returncode - if return_code is not 0: - raise Exception(stderr_str) - except Exception, error: - raise Exception('Error while running SortMeRNA:\n' + - '\n'.join([str(error), stdout_str, stderr_str])) - - -if __name__ == "__main__": - main() diff -r a8ac09e937f3 -r b482293b2987 sortmerna.xml --- a/sortmerna.xml Mon Aug 03 08:18:26 2015 -0400 +++ b/sortmerna.xml Wed Aug 05 02:50:43 2015 -0400 @@ -1,7 +1,7 @@ - + Fast and accurate filtering of ribosomal RNAs in metatranscriptomic data - sortmerna + sortmerna &1|grep 'SortMeRNA version' ]]> - + - - - - - + + + + + + + + + + + + + + + + + + + - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + - - - - - - - - - + + + + + - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + - - outputs_selected and 'accept' in outputs_selected + + aligned_fastx['aligned_fastx_selector'] - outputs_selected and 'other' in outputs_selected + label="Rejected reads on ${on_string} (${input_reads.datatype.file_ext})"> + aligned_fastx['aligned_fastx_selector'] and aligned_fastx['other'] + + + aligned_sam['aligned_sam_selector'] - + + aligned_blast + + + + + log - - + - - - + + + - - + + + + + + + + + + + + + + diff -r a8ac09e937f3 -r b482293b2987 test-data/merged-paired-reads_output.fastq --- a/test-data/merged-paired-reads_output.fastq Mon Aug 03 08:18:26 2015 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,50 +0,0 @@ -@PHOSPHORE_0118:4:1101:5073:2473#ATCACG/1 -CTCGCGCCACTGGCTGGCTCTCCCCTCTCTCTCTCTCTCTCCCCTCTCTCTCCCCGAGTACACCCGCGCGGGC -+PHOSPHORE_0118:4:1101:5073:2473#ATCACG/1 -___cccccggggghhhhf_YbcgecZYa^Y^ceXccacXIXaafX`_]VMHVbdhSTFKU_Z]b]WFW_aZ__ -@PHOSPHORE_0118:4:1101:1143:2236#ATCACG/1 -CGCGCGCACCATGCCCGTGCTCGACTGCTCGCCGCCCGCGCCGCGCCACGCCTCGCCGTTCACCTTGACCTCCACGCGCAGGCC -+PHOSPHORE_0118:4:1101:1143:2236#ATCACG/1 -__beeeecggcgfhhhhhhhhhhfhhhhhhhhhhhgececcccccccacaccccaccccac`bbbcbbbcbbc[`_a[aLT[aa -@PHOSPHORE_0118:4:1101:17780:2432#ATCACG/1 -GCGGTCTGGAACCACCTGATCCCATCCCGAACTCAGCCGTGAAACAGACCAGCGCCGATGGTAGTGTGGCTTCTGCCCGT -+PHOSPHORE_0118:4:1101:17780:2432#ATCACG/1 -bbbeeeeegggggiiiiiiihiiihiiiiiiihiidhhigfhihiihgiiiiiggeccaccbbb_b`baccab`]_ba[T -@PHOSPHORE_0118:4:1101:1397:2156#ATCACG/1 -ATGCGCTTGAGGTAGCCAAGGGCACGCAGTTCACACCGTATGGCGGCGATGCGAGGCGGCGGCAGGAGG -+PHOSPHORE_0118:4:1101:1397:2156#ATCACG/1 -bbbeeeeegggggiiiiiiiiihfhiheghihgihiiihiiiiihifc^acac^_cEHVZaccaac^__ -@PHOSPHORE_0118:4:1101:7859:3729#ATCACG/1 -GCATCGCCGCCCTCGCCCTCGCCCCCGCCGCCCTCGCCCT -+PHOSPHORE_0118:4:1101:7859:3729#ATCACG/1 -___c`cc`Yeee[ePY_`UY^_eGL_F_`dUU\F\^Q\a_ -@PHOSPHORE_0118:4:1101:1633:2146#ATCACG/1 -TGACCAATGTTCTTCATAACGTCCCCGATAGACATGCCAACATTCGGTCCGTTGCGAAGATAGTCATCGCCGTCATAGTC -+PHOSPHORE_0118:4:1101:1633:2146#ATCACG/1 -_bbecdecgggggihhiiiihghhhiihiichdfghiffghfghiifiifhhhiiggaaaddddbdd`bccaQ\a`bc_b -@PHOSPHORE_0118:4:1101:7007:4197#ATCACG/1 -CCCCACTCACCACGCACGGCACCGCACGGCACTCTCACGGTTTCTCTCTCTCTCTCCCATCCTTCTCACTCCTTTCTCTCTCTCCCTATCTCTCTCTCTC -+PHOSPHORE_0118:4:1101:7007:4197#ATCACG/1 -bbbeeeeeggggghiiiiihfhifhhiihiagfhiiiiig^dceeeedcddddcccccRZ``bbcbcbbR_bb`b`bb`bbb`bbcbR]b_]]`bb_b_b -@PHOSPHORE_0118:4:1101:1719:2154#ATCACG/1 -GCATGCCAAACGTGGAGGCGGTGGGCGCGACGCCGGTGGGGAACTTGGACTTGTACATGGCCGGCCGCTGCGTCACCATCGA -+PHOSPHORE_0118:4:1101:1719:2154#ATCACG/1 -bbbeeeeeggggegifhiiii^fgh\edgffeeccc\acccWV_ccccbccccbbdccccccccaaacaZ]acaa^bccccb -@PHOSPHORE_0118:4:1101:7113:4916#ATCACG/1 -AGCCTGACGATAACCTACTTTCACACTGGTTGCAGCACTATCATCGGCGCAAAGTCGTTTCACGGTCCTGTTCGGGATGGGAAGGGGTGGGACCGACTTGC -+PHOSPHORE_0118:4:1101:7113:4916#ATCACG/1 -bbbeeeeegggggiiiiiiiiiihiiiiigiiiiiiiiiiiiiiiiiiiihiiggggeeeeddccccccccccccccccccbbccccHXaacccc]acacb -@PHOSPHORE_0118:4:1101:1730:2193#ATCACG/1 -GTTCGATTAGTCTTTCGCCCCTATACCCAAATTTGACGATCGATTTGCACGTCAGAATCGCTAC -+PHOSPHORE_0118:4:1101:1730:2193#ATCACG/1 -baaeeeeeggggghhifhiiihfgbghhhiihihfhihfhihdfghhiiifhiihdbfgddged -@PHOSPHORE_0118:4:1101:11597:5204#ATCACG/1 -GCCGAAGAAGGACGTGGGAATCTGCGATAAGCCTGGTGGAGTCGATAACCGGACGTTGAGACCAGGATTTCCGAATGGGGAAACCCCGCACGACGTGTCGT -+PHOSPHORE_0118:4:1101:11597:5204#ATCACG/1 -bbbeeeeegggggifgiiiiiiiiiihiihiiiiiicfhgheghiiihihggcecccccccccccccccccccccccccccaacccaccc__[_a[a^[_Q - -@PHOSPHORE_0118:4:1101:13427:6518#ATCACG/1 -GCGCTCGAGGAGGACGTGGCGCCCGAGGCTGCCGAGGCTG -+PHOSPHORE_0118:4:1101:13427:6518#ATCACG/1 -Z__cccc`ecg^_fe^[^^^Y_ffc]_`[eZ_c_W\]]b` - diff -r a8ac09e937f3 -r b482293b2987 test-data/read_small.fasta --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/read_small.fasta Wed Aug 05 02:50:43 2015 -0400 @@ -0,0 +1,2 @@ +>read1 +GCCATTGGAAACTGGGAGACTTGAGTGCAGAAGAGGAGAGTGGAATACCATGTGTAGCGGTGAAATGCGTAGATATATGGAGGAACACC diff -r a8ac09e937f3 -r b482293b2987 test-data/read_small.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/read_small.fastq Wed Aug 05 02:50:43 2015 -0400 @@ -0,0 +1,4 @@ +@read1 +GCCATTGGAAACTGGGAGACTTGAGTGCAGAAGAGGAGAGTGGAATACCATGTGTAGCGGTGAAATGCGTAGATATATGGAGGAACACC ++read1 +IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII diff -r a8ac09e937f3 -r b482293b2987 test-data/ref_small.fasta --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/ref_small.fasta Wed Aug 05 02:50:43 2015 -0400 @@ -0,0 +1,2 @@ +>EncFa169 count=1; cluster_weight=27830; cluster=EncFa169; cluster_score=1.000000; cluster_center=True; +AGAGTTTGATCCTGGCTCAGGACGAACGCTGGCGGCGTGCCTAATACATGCAAGTCGAACGCTTCTTTCCTCCCGAGTGCTTGCACTCAATTGGAAAGAGGAGTGGCGGACGGGTGAGTAACACGTGGGTAACCTACCCATCAGAGGGGGATAACACTTGGAAACAGGTGCTAATACCGCATAACAGTTTATGCCGCATGGCATAAGAGTGAAAGGCGCTTTCGGGTGTCGCTGATGGATGGACCCGCGGTGCATTAGCTAGTTGGTGAGGTAACGGCTCACCAAGGCCACGATGCATAGCCGACCTGAGAGGGTGATCGGCCACACTGGGACTGAGACACGGCCCAGACTCCTACGGGAGGCAGCAGTAGGGAATCTTCGGCAATGGACGAAAGTCTGACCGAGCAACGCCGCGTGAGTGAAGAAGGTTTTCGGATCGTAAAACTCTGTTGTTAGAGAAGAACAAGGACGTTAGTAACTGAACGTCCCCTGACGGTATCTAACCAGAAAGCCACGGCTAACTACGTGCCAGCAGCCGCGGTAATACGTAGGTGGCAAGCGTTGTCCGGATTTATTGGGCGTAAAGCGAGCGCAGGCGGTTTCTTAAGTCTGATGTGAAAGCCCCCGGCTCAACCGGGGAGGGTCATTGGAAACTGGGAGACTTGAGTGCAGAAGAGGAGAGTGGAATTCCATGTGTAGCGGTGAAATGCGTAGATATATGGAGGAACACCAGTGGCGAAGGCGGCTCTCTGGTCTGTAACTGACGCTGAGGCTCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGAGTGCTAAGTGTTGGAGGGTTTCCGCCCTTCAGTGCTGCAGCAAACGCATTAAGCACTCCGCCTGGGGAGTACGACCGCAAGGTTGAAACTCAAAGGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGAAGCAACGCGAAGAACCTTACCAGGTCTTGACATCCTTTGACCACTCTAGAGATAGAGCTTTCCCTTCGGGGACAAAGTGACAGGTGGTGCATGGTTGTCGTCAGCTCGTGTCGTGAGATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATTGTTAGTTGCCATCATTTAGTTGGGCACTCTAGCGAGACTGCCGGTGACAAACCGGAGGAAGGTGGGGATGACGTCAAATCATCATGCCCCTTATGACCTGGGCTACACACGTGCTACAATGGGAAGTACAACGAGTCGCTAGACCGCGAGGTCATGCAAATCTCTTAAAGCTTCTCTCAGTTCGGATTGCAGGCTGCAACTCGCCTGCATGAAGCCGGAATCGCTAGTAATCGCGGATCAGCACGCCGCGGTGAATACGTTCCCGGGCCTTGTACACACCGCCCGTCACACCACGAGAGTTTGTAACACCCGAAGTCGGTGAGGTAACCTTTTTGGAGCCAGCCGCCTAAGGTGGGATAGATGATTGGGGTGAAGTCGTAACAAGGTAGCCGTATCGGAAGGTGCGGCTGGATCACCT diff -r a8ac09e937f3 -r b482293b2987 test-data/sortmerna_wrapper_accept1.fastq --- a/test-data/sortmerna_wrapper_accept1.fastq Mon Aug 03 08:18:26 2015 -0400 +++ b/test-data/sortmerna_wrapper_accept1.fastq Wed Aug 05 02:50:43 2015 -0400 @@ -1,28 +1,4 @@ -@PHOSPHORE_0118:4:1101:5073:2473#ATCACG/1 -CTCGCGCCACTGGCTGGCTCTCCCCTCTCTCTCTCTCTCTCCCCTCTCTCTCCCCGAGTACACCCGCGCGGGC -+PHOSPHORE_0118:4:1101:5073:2473#ATCACG/1 -___cccccggggghhhhf_YbcgecZYa^Y^ceXccacXIXaafX`_]VMHVbdhSTFKU_Z]b]WFW_aZ__ -@PHOSPHORE_0118:4:1101:17780:2432#ATCACG/1 -GCGGTCTGGAACCACCTGATCCCATCCCGAACTCAGCCGTGAAACAGACCAGCGCCGATGGTAGTGTGGCTTCTGCCCGT -+PHOSPHORE_0118:4:1101:17780:2432#ATCACG/1 -bbbeeeeegggggiiiiiiihiiihiiiiiiihiidhhigfhihiihgiiiiiggeccaccbbb_b`baccab`]_ba[T -@PHOSPHORE_0118:4:1101:7859:3729#ATCACG/1 -GCATCGCCGCCCTCGCCCTCGCCCCCGCCGCCCTCGCCCT -+PHOSPHORE_0118:4:1101:7859:3729#ATCACG/1 -___c`cc`Yeee[ePY_`UY^_eGL_F_`dUU\F\^Q\a_ -@PHOSPHORE_0118:4:1101:7007:4197#ATCACG/1 -CCCCACTCACCACGCACGGCACCGCACGGCACTCTCACGGTTTCTCTCTCTCTCTCCCATCCTTCTCACTCCTTTCTCTCTCTCCCTATCTCTCTCTCTC -+PHOSPHORE_0118:4:1101:7007:4197#ATCACG/1 -bbbeeeeeggggghiiiiihfhifhhiihiagfhiiiiig^dceeeedcddddcccccRZ``bbcbcbbR_bb`b`bb`bbb`bbcbR]b_]]`bb_b_b -@PHOSPHORE_0118:4:1101:7113:4916#ATCACG/1 -AGCCTGACGATAACCTACTTTCACACTGGTTGCAGCACTATCATCGGCGCAAAGTCGTTTCACGGTCCTGTTCGGGATGGGAAGGGGTGGGACCGACTTGC -+PHOSPHORE_0118:4:1101:7113:4916#ATCACG/1 -bbbeeeeegggggiiiiiiiiiihiiiiigiiiiiiiiiiiiiiiiiiiihiiggggeeeeddccccccccccccccccccbbccccHXaacccc]acacb -@PHOSPHORE_0118:4:1101:11597:5204#ATCACG/1 -GCCGAAGAAGGACGTGGGAATCTGCGATAAGCCTGGTGGAGTCGATAACCGGACGTTGAGACCAGGATTTCCGAATGGGGAAACCCCGCACGACGTGTCGT -+PHOSPHORE_0118:4:1101:11597:5204#ATCACG/1 -bbbeeeeegggggifgiiiiiiiiiihiihiiiiiicfhgheghiiihihggcecccccccccccccccccccccccccccaacccaccc__[_a[a^[_Q -@PHOSPHORE_0118:4:1101:13427:6518#ATCACG/1 -GCGCTCGAGGAGGACGTGGCGCCCGAGGCTGCCGAGGCTG -+PHOSPHORE_0118:4:1101:13427:6518#ATCACG/1 -Z__cccc`ecg^_fe^[^^^Y_ffc]_`[eZ_c_W\]]b` +@read1 +GCCATTGGAAACTGGGAGACTTGAGTGCAGAAGAGGAGAGTGGAATACCATGTGTAGCGGTGAAATGCGTAGATATATGGAGGAACACC ++read1 +IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII diff -r a8ac09e937f3 -r b482293b2987 test-data/sortmerna_wrapper_accept2.fasta --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/sortmerna_wrapper_accept2.fasta Wed Aug 05 02:50:43 2015 -0400 @@ -0,0 +1,2 @@ +>read1 +GCCATTGGAAACTGGGAGACTTGAGTGCAGAAGAGGAGAGTGGAATACCATGTGTAGCGGTGAAATGCGTAGATATATGGAGGAACACC diff -r a8ac09e937f3 -r b482293b2987 test-data/sortmerna_wrapper_in1.fastq --- a/test-data/sortmerna_wrapper_in1.fastq Mon Aug 03 08:18:26 2015 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,48 +0,0 @@ -@PHOSPHORE_0118:4:1101:1143:2236#ATCACG/1 -CGCGCGCACCATGCCCGTGCTCGACTGCTCGCCGCCCGCGCCGCGCCACGCCTCGCCGTTCACCTTGACCTCCACGCGCAGGCC -+PHOSPHORE_0118:4:1101:1143:2236#ATCACG/1 -__beeeecggcgfhhhhhhhhhhfhhhhhhhhhhhgececcccccccacaccccaccccac`bbbcbbbcbbc[`_a[aLT[aa -@PHOSPHORE_0118:4:1101:1397:2156#ATCACG/1 -ATGCGCTTGAGGTAGCCAAGGGCACGCAGTTCACACCGTATGGCGGCGATGCGAGGCGGCGGCAGGAGG -+PHOSPHORE_0118:4:1101:1397:2156#ATCACG/1 -bbbeeeeegggggiiiiiiiiihfhiheghihgihiiihiiiiihifc^acac^_cEHVZaccaac^__ -@PHOSPHORE_0118:4:1101:5073:2473#ATCACG/1 -CTCGCGCCACTGGCTGGCTCTCCCCTCTCTCTCTCTCTCTCCCCTCTCTCTCCCCGAGTACACCCGCGCGGGC -+PHOSPHORE_0118:4:1101:5073:2473#ATCACG/1 -___cccccggggghhhhf_YbcgecZYa^Y^ceXccacXIXaafX`_]VMHVbdhSTFKU_Z]b]WFW_aZ__ -@PHOSPHORE_0118:4:1101:17780:2432#ATCACG/1 -GCGGTCTGGAACCACCTGATCCCATCCCGAACTCAGCCGTGAAACAGACCAGCGCCGATGGTAGTGTGGCTTCTGCCCGT -+PHOSPHORE_0118:4:1101:17780:2432#ATCACG/1 -bbbeeeeegggggiiiiiiihiiihiiiiiiihiidhhigfhihiihgiiiiiggeccaccbbb_b`baccab`]_ba[T -@PHOSPHORE_0118:4:1101:7859:3729#ATCACG/1 -GCATCGCCGCCCTCGCCCTCGCCCCCGCCGCCCTCGCCCT -+PHOSPHORE_0118:4:1101:7859:3729#ATCACG/1 -___c`cc`Yeee[ePY_`UY^_eGL_F_`dUU\F\^Q\a_ -@PHOSPHORE_0118:4:1101:1633:2146#ATCACG/1 -TGACCAATGTTCTTCATAACGTCCCCGATAGACATGCCAACATTCGGTCCGTTGCGAAGATAGTCATCGCCGTCATAGTC -+PHOSPHORE_0118:4:1101:1633:2146#ATCACG/1 -_bbecdecgggggihhiiiihghhhiihiichdfghiffghfghiifiifhhhiiggaaaddddbdd`bccaQ\a`bc_b -@PHOSPHORE_0118:4:1101:1719:2154#ATCACG/1 -GCATGCCAAACGTGGAGGCGGTGGGCGCGACGCCGGTGGGGAACTTGGACTTGTACATGGCCGGCCGCTGCGTCACCATCGA -+PHOSPHORE_0118:4:1101:1719:2154#ATCACG/1 -bbbeeeeeggggegifhiiii^fgh\edgffeeccc\acccWV_ccccbccccbbdccccccccaaacaZ]acaa^bccccb -@PHOSPHORE_0118:4:1101:7007:4197#ATCACG/1 -CCCCACTCACCACGCACGGCACCGCACGGCACTCTCACGGTTTCTCTCTCTCTCTCCCATCCTTCTCACTCCTTTCTCTCTCTCCCTATCTCTCTCTCTC -+PHOSPHORE_0118:4:1101:7007:4197#ATCACG/1 -bbbeeeeeggggghiiiiihfhifhhiihiagfhiiiiig^dceeeedcddddcccccRZ``bbcbcbbR_bb`b`bb`bbb`bbcbR]b_]]`bb_b_b -@PHOSPHORE_0118:4:1101:7113:4916#ATCACG/1 -AGCCTGACGATAACCTACTTTCACACTGGTTGCAGCACTATCATCGGCGCAAAGTCGTTTCACGGTCCTGTTCGGGATGGGAAGGGGTGGGACCGACTTGC -+PHOSPHORE_0118:4:1101:7113:4916#ATCACG/1 -bbbeeeeegggggiiiiiiiiiihiiiiigiiiiiiiiiiiiiiiiiiiihiiggggeeeeddccccccccccccccccccbbccccHXaacccc]acacb -@PHOSPHORE_0118:4:1101:11597:5204#ATCACG/1 -GCCGAAGAAGGACGTGGGAATCTGCGATAAGCCTGGTGGAGTCGATAACCGGACGTTGAGACCAGGATTTCCGAATGGGGAAACCCCGCACGACGTGTCGT -+PHOSPHORE_0118:4:1101:11597:5204#ATCACG/1 -bbbeeeeegggggifgiiiiiiiiiihiihiiiiiicfhgheghiiihihggcecccccccccccccccccccccccccccaacccaccc__[_a[a^[_Q -@PHOSPHORE_0118:4:1101:1730:2193#ATCACG/1 -GTTCGATTAGTCTTTCGCCCCTATACCCAAATTTGACGATCGATTTGCACGTCAGAATCGCTAC -+PHOSPHORE_0118:4:1101:1730:2193#ATCACG/1 -baaeeeeeggggghhifhiiihfgbghhhiihihfhihfhihdfghhiiifhiihdbfgddged -@PHOSPHORE_0118:4:1101:13427:6518#ATCACG/1 -GCGCTCGAGGAGGACGTGGCGCCCGAGGCTGCCGAGGCTG -+PHOSPHORE_0118:4:1101:13427:6518#ATCACG/1 -Z__cccc`ecg^_fe^[^^^Y_ffc]_`[eZ_c_W\]]b` diff -r a8ac09e937f3 -r b482293b2987 test-data/sortmerna_wrapper_other1.fastq --- a/test-data/sortmerna_wrapper_other1.fastq Mon Aug 03 08:18:26 2015 -0400 +++ b/test-data/sortmerna_wrapper_other1.fastq Wed Aug 05 02:50:43 2015 -0400 @@ -1,20 +0,0 @@ -@PHOSPHORE_0118:4:1101:1143:2236#ATCACG/1 -CGCGCGCACCATGCCCGTGCTCGACTGCTCGCCGCCCGCGCCGCGCCACGCCTCGCCGTTCACCTTGACCTCCACGCGCAGGCC -+PHOSPHORE_0118:4:1101:1143:2236#ATCACG/1 -__beeeecggcgfhhhhhhhhhhfhhhhhhhhhhhgececcccccccacaccccaccccac`bbbcbbbcbbc[`_a[aLT[aa -@PHOSPHORE_0118:4:1101:1397:2156#ATCACG/1 -ATGCGCTTGAGGTAGCCAAGGGCACGCAGTTCACACCGTATGGCGGCGATGCGAGGCGGCGGCAGGAGG -+PHOSPHORE_0118:4:1101:1397:2156#ATCACG/1 -bbbeeeeegggggiiiiiiiiihfhiheghihgihiiihiiiiihifc^acac^_cEHVZaccaac^__ -@PHOSPHORE_0118:4:1101:1633:2146#ATCACG/1 -TGACCAATGTTCTTCATAACGTCCCCGATAGACATGCCAACATTCGGTCCGTTGCGAAGATAGTCATCGCCGTCATAGTC -+PHOSPHORE_0118:4:1101:1633:2146#ATCACG/1 -_bbecdecgggggihhiiiihghhhiihiichdfghiffghfghiifiifhhhiiggaaaddddbdd`bccaQ\a`bc_b -@PHOSPHORE_0118:4:1101:1719:2154#ATCACG/1 -GCATGCCAAACGTGGAGGCGGTGGGCGCGACGCCGGTGGGGAACTTGGACTTGTACATGGCCGGCCGCTGCGTCACCATCGA -+PHOSPHORE_0118:4:1101:1719:2154#ATCACG/1 -bbbeeeeeggggegifhiiii^fgh\edgffeeccc\acccWV_ccccbccccbbdccccccccaaacaZ]acaa^bccccb -@PHOSPHORE_0118:4:1101:1730:2193#ATCACG/1 -GTTCGATTAGTCTTTCGCCCCTATACCCAAATTTGACGATCGATTTGCACGTCAGAATCGCTAC -+PHOSPHORE_0118:4:1101:1730:2193#ATCACG/1 -baaeeeeeggggghhifhiiihfgbghhhiihihfhihfhihdfghhiiifhiihdbfgddged diff -r a8ac09e937f3 -r b482293b2987 test-data/sortmerna_wrapper_other2.fasta diff -r a8ac09e937f3 -r b482293b2987 test-data/sortmerna_wrapper_sam1.sam --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/sortmerna_wrapper_sam1.sam Wed Aug 05 02:50:43 2015 -0400 @@ -0,0 +1,3 @@ +@HD VN:1.0 SO:unsorted +@PG ID:sortmerna VN:1.0 CL:sortmerna --ref /tmp/tmpY80cK0/files/000/dataset_2.dat,dataset_2 --reads /tmp/tmpY80cK0/files/000/dataset_1.dat --aligned aligned --fastx --other other_file.dat dat -a 1 +read1 0 EncFa169 645 255 2S87M * 0 0 GCCATTGGAAACTGGGAGACTTGAGTGCAGAAGAGGAGAGTGGAATACCATGTGTAGCGGTGAAATGCGTAGATATATGGAGGAACACC IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII AS:i:169 NM:i:1 diff -r a8ac09e937f3 -r b482293b2987 test-data/sortmerna_wrapper_sam2.sam --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/sortmerna_wrapper_sam2.sam Wed Aug 05 02:50:43 2015 -0400 @@ -0,0 +1,3 @@ +@HD VN:1.0 SO:unsorted +@PG ID:sortmerna VN:1.0 CL:sortmerna --ref /tmp/tmpY80cK0/files/000/dataset_7.dat,dataset_7 --reads /tmp/tmpY80cK0/files/000/dataset_6.dat --aligned aligned --fastx --other other_file.dat dat -a 1 +read1 0 EncFa169 645 255 2S87M * 0 0 GCCATTGGAAACTGGGAGACTTGAGTGCAGAAGAGGAGAGTGGAATACCATGTGTAGCGGTGAAATGCGTAGATATATGGAGGAACACC * AS:i:169 NM:i:1 diff -r a8ac09e937f3 -r b482293b2987 tool-data/rRNA_databases.loc --- a/tool-data/rRNA_databases.loc Mon Aug 03 08:18:26 2015 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,8 +0,0 @@ -rfam-5.8s Database Rfam 5.8s $SORTMERNADIR/rRNA_databases/rfam-5.8s-database-id98.fasta -rfam-5s Database Rfam 5s $SORTMERNADIR/rRNA_databases/rfam-5s-database-id98.fasta -silva-arc-16s Database Silva-Arc 16s $SORTMERNADIR/rRNA_databases/silva-arc-16s-database-id95.fasta -silva-arc-23s Database Silva-Arc 23s $SORTMERNADIR/rRNA_databases/silva-arc-23s-database-id98.fasta -silva-bac-16s Database Silva-Bac 16s $SORTMERNADIR/rRNA_databases/silva-bac-16s-database-id85.fasta -silva-bac-23s Database Silva-Bac 23s $SORTMERNADIR/rRNA_databases/silva-bac-23s-database-id98.fasta -silva-euk-18s Databse Silva-Euk 18s $SORTMERNADIR/rRNA_databases/silva-euk-18s-database-id95.fasta -silva-euk-28s Database Silva-Euk 28s $SORTMERNADIR/rRNA_databases/silva-euk-28s-database-id98.fasta diff -r a8ac09e937f3 -r b482293b2987 tool-data/rRNA_databases.loc.sample --- a/tool-data/rRNA_databases.loc.sample Mon Aug 03 08:18:26 2015 -0400 +++ b/tool-data/rRNA_databases.loc.sample Wed Aug 05 02:50:43 2015 -0400 @@ -1,26 +1,30 @@ #This is a sample file distributed with Galaxy that is used to define a -#list of public ribosomal databases, using three columns tab separated -#(longer whitespace are TAB characters): +#list of public ribosomal databases for SortMeRNA, using the following format +#(white space characters are TAB characters): # -# +# # -#It is important that the actual database name does not have a space in it, -#and that the first tab that appears in the line is right before the path. +#So, for example, if your database is rfam-5.8s-id98 and the path to your FASTA +#file is /data/rRNA_databases/rfam-5.8s-id98.fasta, then the rRNA_databases.loc +#entry would look like this: # -#So, for example, if your database is rfam-5.8s and the path to your base name -#is /data/rRNA_databases/rfam-5.8s, then the rRNA_databases.loc entry would look like this: +#rfam-5.8s-id98 Rfam 5.8S eukarya /data/rRNA_databases/rfam-5.8s-id98.fasta # -#rfam-5.8s Rfam 5.8S eukarya /data/rRNA_databases/rfam-5.8s +#For each rRNA database, you need to create the index files using the +#indexdb_rna program provided by SortMeRNA. You need to specify as index +#basename the path of the FASTA file without extension. For example, for the +#previous database the command is: # -#Since SortMeRNA comes bundled with eight ribosomal databases, which are ready -#for use after the tool installation, this sample file is in fact an actual file -#to save the user the trouble of setting it. +# indexdb_rna --ref /data/rRNA_databases/rfam-5.8s-id98.fasta,/data/rRNA_databases/rfam-5.8s-id98 # -rfam-5.8s Rfam 5.8S eukarya $SORTMERNADIR/rRNA_databases/rfam-5.8s-database-id98.fasta -rfam-5s Rfam 5S archaea/bacteria $SORTMERNADIR/rRNA_databases/rfam-5s-database-id98.fasta -silva-arc-16s SILVA 16S archaea $SORTMERNADIR/rRNA_databases/silva-arc-16s-database-id95.fasta -silva-arc-23s SILVA 16S bacteria $SORTMERNADIR/rRNA_databases/silva-arc-23s-database-id98.fasta -silva-bac-16s SILVA 16S bacteria $SORTMERNADIR/rRNA_databases/silva-bac-16s-database-id85.fasta -silva-bac-23s SILVA 23S bacteria $SORTMERNADIR/rRNA_databases/silva-bac-23s-database-id98.fasta -silva-euk-18s SILVA 18S eukarya $SORTMERNADIR/rRNA_databases/silva-euk-18s-database-id95.fasta -silva-euk-28s SILVA 28S eukarya $SORTMERNADIR/rRNA_databases/silva-euk-28s-database-id98.fasta +#Since SortMeRNA comes bundled with eight ribosomal databases, you can use them +#by creating the actual index files as explained above and uncommenting the +#following lines. +#rfam-5.8s-id98 Rfam 5.8S eukarya $SORTMERNADIR/rRNA_databases/rfam-5.8s-database-id98.fasta +#rfam-5s-id98 Rfam 5S archaea/bacteria $SORTMERNADIR/rRNA_databases/rfam-5s-database-id98.fasta +#silva-arc-16s-id95 SILVA v.119 16S archaea $SORTMERNADIR/rRNA_databases/silva-arc-16s-id95.fasta +#silva-arc-23s-id98 SILVA v.119 23S archaea $SORTMERNADIR/rRNA_databases/silva-arc-23s-id98.fasta +#silva-bac-16s-id90 SILVA v.119 16S bacteria $SORTMERNADIR/rRNA_databases/silva-bac-16s-id90.fasta +#silva-bac-23s-id98 SILVA v.119 23S bacteria $SORTMERNADIR/rRNA_databases/silva-bac-23s-id98.fasta +#silva-euk-18s-id95 SILVA v.119 18S eukarya $SORTMERNADIR/rRNA_databases/silva-euk-18s-id95.fasta +#silva-euk-28s-id98 SILVA v.119 28S eukarya $SORTMERNADIR/rRNA_databases/silva-euk-28s-id98.fasta diff -r a8ac09e937f3 -r b482293b2987 tool_dependencies.xml --- a/tool_dependencies.xml Mon Aug 03 08:18:26 2015 -0400 +++ b/tool_dependencies.xml Wed Aug 05 02:50:43 2015 -0400 @@ -1,9 +1,9 @@ - + - http://bioinfo.lifl.fr/RNA/sortmerna/code/sortmerna-1.9.tar.gz + https://github.com/biocore/sortmerna/archive/2.0.tar.gz $INSTALL_DIR