# HG changeset patch # User bgruening # Date 1761734739 0 # Node ID 58b9ef13fcca48912096383cba2959e2723dce1c # Parent 3c3cf3ff1a914b90a169c809bf7645236512374a planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/salmon commit 1971e35a365a91ae36c990b374a552a46196e888 diff -r 3c3cf3ff1a91 -r 58b9ef13fcca macros.xml --- a/macros.xml Fri Jun 23 17:01:36 2023 +0000 +++ b/macros.xml Wed Oct 29 10:45:39 2025 +0000 @@ -1,8 +1,13 @@ 1.10.1 - 0 + 3 q7 - 20.01 + 24.1 + + + salmon + + salmon @@ -14,8 +19,11 @@ scipy + + + - + @@ -23,7 +31,7 @@ - + @@ -63,7 +71,7 @@ - + @@ -73,14 +81,11 @@
- - + - - +
@@ -90,75 +95,62 @@ - - - + - - - + + - - - - - - - - - - - - - + + + + - - - - - - - - - - - - + @@ -171,152 +163,135 @@ - -
- + - + - - - - - - - - - - -
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- - - - - - - - - - - - - - 10.1186/s13059-019-1670-y @@ -327,26 +302,26 @@ " | cut -d " " -f 1 > 'decoys.txt' && + #if $TRANSCRIPTSOURCE_CONDITIONAL.TranscriptSource != "indexed": + #if $TRANSCRIPTSOURCE_CONDITIONAL.s_index.genome + cat '${TRANSCRIPTSOURCE_CONDITIONAL.s_index.genome}' | grep "^>" | cut -d " " -f 1 > 'decoys.txt' && sed -i.bak -e 's/>//g' 'decoys.txt' && - cat '${quant_type.refTranscriptSource.s_index.fasta}' '${quant_type.refTranscriptSource.s_index.genome}' > 'input_index.fasta' && + cat '${TRANSCRIPTSOURCE_CONDITIONAL.s_index.fasta}' '${TRANSCRIPTSOURCE_CONDITIONAL.s_index.genome}' > 'input_index.fasta' && #else - ln -s '${quant_type.refTranscriptSource.s_index.fasta}' 'input_index.fasta' && + ln -s '${TRANSCRIPTSOURCE_CONDITIONAL.s_index.fasta}' 'input_index.fasta' && #end if - salmon index -i ./index - --kmerLen '${quant_type.refTranscriptSource.s_index.kmer}' + salmon --no-version-check index -i ./index + --kmerLen '${TRANSCRIPTSOURCE_CONDITIONAL.s_index.kmer}' --gencode - --threads "\${GALAXY_SLOTS:-4}" + --threads "\$SLOTS" --transcripts 'input_index.fasta' - #if $quant_type.refTranscriptSource.s_index.genome + #if $TRANSCRIPTSOURCE_CONDITIONAL.s_index.genome --decoy 'decoys.txt' #end if && #set $index_path = './index' #else - #set $index_path = $quant_type.refTranscriptSource.index.fields.path + #set $index_path = $TRANSCRIPTSOURCE_CONDITIONAL.index.fields.path #end if ]]> @@ -356,58 +331,33 @@ #if $quant_type.input.single_or_paired.input_singles.ext == 'fasta': #set $ext = 'fasta' #else: - #if $quant_type.input.single_or_paired.input_singles.is_of_type("fastq.gz", "fastqsanger.gz"): + #if $quant_type.input.single_or_paired.input_singles.is_of_type("fastqsanger.gz"): #set compressed = 'GZ' - #else if $quant_type.input.single_or_paired.input_singles.is_of_type("fastq.bz2", "fastqsanger.bz2"): + #else if $quant_type.input.single_or_paired.input_singles.is_of_type("fastqsanger.bz2"): #set compressed = 'BZ2' #end if #set $ext = 'fastq' #end if ln -s $quant_type.input.single_or_paired.input_singles ./single.$ext && - #else if $quant_type.input.single_or_paired.single_or_paired_opts == 'paired': - #if $quant_type.input.single_or_paired.input_mate1.ext == 'fasta': - #set $ext = 'fasta' - #else: - #if $quant_type.input.single_or_paired.input_mate1.is_of_type("fastq.gz", "fastqsanger.gz"): - #set compressed = 'GZ' - #else if $quant_type.input.single_or_paired.input_mate1.is_of_type("fastq.bz2", "fastqsanger.bz2"): - #set compressed = 'BZ2' - #end if - #set $ext = 'fastq' - #end if - ln -s $quant_type.input.single_or_paired.input_mate1 ./mate1.$ext && - ln -s $quant_type.input.single_or_paired.input_mate2 ./mate2.$ext && #else if $quant_type.input.single_or_paired.single_or_paired_opts == 'paired_collection': #if $quant_type.input.single_or_paired.input_1.forward.ext == 'fasta': #set $ext = 'fasta' #else: - #if $quant_type.input.single_or_paired.input_1.forward.is_of_type("fastq.gz", "fastqsanger.gz"): + #if $quant_type.input.single_or_paired.input_1.forward.is_of_type("fastqsanger.gz"): #set compressed = 'GZ' - #else if $quant_type.input.single_or_paired.input_1.forward.is_of_type("fastq.bz2", "fastqsanger.bz2"): + #else if $quant_type.input.single_or_paired.input_1.forward.is_of_type("fastqsanger.bz2"): #set compressed = 'BZ2' #end if #set $ext = 'fastq' #end if ln -s '${quant_type.input.single_or_paired.input_1.forward}' ./mate1.$ext && ln -s '${quant_type.input.single_or_paired.input_1.reverse}' ./mate2.$ext && - #else if '$quant_type.input.single_or_paired.single_or_paired_opts' == 'paired_interleaved': - #if $quant_type.input.single_or_paired.input_1.ext == 'fasta': - #set $ext = 'fasta' - #else: - #if $quant_type.input.single_or_paired.input_1.is_of_type("fastq.gz", "fastqsanger.gz"): - #set compressed = 'GZ' - #else if $quant_type.input.single_or_paired.input_1.is_of_type("fastq.bz2", "fastqsanger.bz2"): - #set compressed = 'BZ2' - #end if - #set $ext = 'fastq' - #end if - ln -s '$quant_type.input.single_or_paired.input_1' ./mate1.$ext && #end if #if $geneMap: ln -s '$geneMap' ./geneMap.${geneMap.ext} && #end if - salmon quant + salmon --no-version-check quant --index '$index_path' #if $quant_type.input.single_or_paired.single_or_paired_opts == 'single': --libType ${quant_type.input.single_or_paired.libtype.strandedness} @@ -425,28 +375,18 @@ #else ${quant_type.input.single_or_paired.libtype.orientation}${quant_type.input.single_or_paired.libtype.strandedness} #end if - #if $quant_type.input.single_or_paired.single_or_paired_opts == 'paired_interleaved': - #if $compressed == 'BZ2': - --mates1 <(bzcat < ./mate1.$ext | seqtk seq -1) - --mates2 <(bzcat < ./mate1.$ext | seqtk seq -2) - #else: - --mates1 <(seqtk seq -1 ./mate1.$ext) - --mates2 <(seqtk seq -2 ./mate1.$ext) - #end if + #if $compressed == 'GZ': + --mates1 <(zcat < ./mate1.$ext) + --mates2 <(zcat < ./mate2.$ext) + #else if $compressed == 'BZ2': + --mates1 <(bzcat < ./mate1.$ext) + --mates2 <(bzcat < ./mate2.$ext) #else: - #if $compressed == 'GZ': - --mates1 <(zcat < ./mate1.$ext) - --mates2 <(zcat < ./mate2.$ext) - #else if $compressed == 'BZ2': - --mates1 <(bzcat < ./mate1.$ext) - --mates2 <(bzcat < ./mate2.$ext) - #else: - --mates1 ./mate1.$ext - --mates2 ./mate2.$ext - #end if + --mates1 ./mate1.$ext + --mates2 ./mate2.$ext #end if #end if - --threads "\${GALAXY_SLOTS:-4}" + --threads "\$SLOTS" ${quant_type.discardOrphansQuasi} ${quant_type.validmap.validateMappings} #if $quant_type.validmap.validateMappings: @@ -491,12 +431,12 @@ #if $geneMap: ln -s "$geneMap" ./geneMap.${geneMap.ext} && #end if - salmon quant - -t '${quant_type.transcript}' - -l '${quant_type.libtype.strandedness}' - -a '${quant_type.afile}' + salmon --no-version-check quant + --targets '${quant_type.transcript}' + --libType '${quant_type.libtype.strandedness}' + --alignments '${quant_type.afile}' $quant_type.ont - --threads "\${GALAXY_SLOTS:-4}" + --threads "\$SLOTS" ${quant_type.discardOrphans} ${quant_type.noErrorModel} #if $quant_type.numErrorBins: @@ -584,33 +524,7 @@ -o ./output ]]> - - macros.xml
+ + - + + + + + + + + + + + + diff -r 3c3cf3ff1a91 -r 58b9ef13fcca test-data/full_data_structure.txt --- a/test-data/full_data_structure.txt Fri Jun 23 17:01:36 2023 +0000 +++ b/test-data/full_data_structure.txt Wed Oct 29 10:45:39 2025 +0000 @@ -544,7 +544,7 @@ 8 198 199 200 201 202 203 204 205 1 1 32 1 TACGGGCTGC 1 1 170 1 1 21 1 GAACATTTGT 1 3 29 30 31 2 2 34 1 CGGCGTCCGG 1 1 1 AGCGGCTGGT 1 -11 5 6 7 8 9 10 11 12 13 14 16 1 1 93 1 TTGACACCAT 1 +1 16 1 1 93 1 TTGACACCAT 1 1 146 1 1 46 1 AACACAGATA 1 4 287 289 291 293 2 2 140 1 CCGTACGGGG 1 35 1 GGATTCAGTA 1 3 273 274 275 1 1 86 1 CCGACTCGTT 1 diff -r 3c3cf3ff1a91 -r 58b9ef13fcca umiout.sh --- a/umiout.sh Fri Jun 23 17:01:36 2023 +0000 +++ b/umiout.sh Wed Oct 29 10:45:39 2025 +0000 @@ -1,11 +1,14 @@ #!/usr/bin/bash -mkdir fixed; +set -e + +mkdir -p fixed; for file in ./umiout/*; -do prefix="${file%.dot.gz}"; -prefix=${prefix/.\/umiout\//}; -gunzip $file; -sed "s/umiout\/$prefix.dot.gz/$prefix/" umiout/$prefix.dot > fixed/$prefix.dot; -dot -Tpdf fixed/$prefix.dot -o fixed/$prefix.pdf; +do +prefix="${file%.dot.gz}"; +prefix="${prefix/.\/umiout\//}"; +gunzip "$file"; +sed "s/umiout\/$prefix.dot.gz/$prefix/" "umiout/$prefix.dot" > "fixed/$prefix.dot"; +dot -Tpdf "fixed/$prefix.dot" -o "fixed/$prefix.pdf"; done ls fixed