# HG changeset patch # User bgruening # Date 1761734807 0 # Node ID 7643953b19f975b003f83a35d078c1582384c65f # Parent a896b8b7d26455bd66927d0baf74ffdb2d47efc7 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/salmon commit 1971e35a365a91ae36c990b374a552a46196e888 diff -r a896b8b7d264 -r 7643953b19f9 macros.xml --- a/macros.xml Mon Mar 18 09:35:17 2024 +0000 +++ b/macros.xml Wed Oct 29 10:46:47 2025 +0000 @@ -1,8 +1,8 @@ 1.10.1 - 2 + 3 q7 - 20.01 + 24.1 salmon @@ -19,8 +19,11 @@ scipy + + + - + @@ -28,7 +31,7 @@ - + @@ -68,7 +71,7 @@ - + @@ -78,14 +81,11 @@
- - + - - +
@@ -95,34 +95,20 @@ - - - + - - - + + - - - - - - - - - - - - @@ -131,40 +117,40 @@ - - - - - - - - - - - - + @@ -177,152 +163,135 @@ - -
- + - + - - - - - - - - - - -
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
- - - - - - - - - - - - - - 10.1186/s13059-019-1670-y @@ -333,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 ]]> @@ -362,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} @@ -431,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: @@ -497,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: @@ -590,33 +524,7 @@ -o ./output ]]> - - + - + + + + + + + + + + + + diff -r a896b8b7d264 -r 7643953b19f9 test-data/full_data_structure.txt --- a/test-data/full_data_structure.txt Mon Mar 18 09:35:17 2024 +0000 +++ b/test-data/full_data_structure.txt Wed Oct 29 10:46:47 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 a896b8b7d264 -r 7643953b19f9 umiout.sh --- a/umiout.sh Mon Mar 18 09:35:17 2024 +0000 +++ b/umiout.sh Wed Oct 29 10:46:47 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