# HG changeset patch # User bgruening # Date 1761734829 0 # Node ID 4dbb23a9fcccb58d26ef72f9c0ce19a98e9c732e # Parent b0fe3e482cb74c96b09c113da185a5df9fc9db42 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/salmon commit 1971e35a365a91ae36c990b374a552a46196e888 diff -r b0fe3e482cb7 -r 4dbb23a9fccc alevin.xml --- a/alevin.xml Mon Oct 28 17:50:35 2024 +0000 +++ b/alevin.xml Wed Oct 29 10:47:09 2025 +0000 @@ -5,63 +5,55 @@ + 1 ? \${GALAXY_SLOTS:-4} - 1 : 1 )); + + #set TRANSCRIPTSOURCE_CONDITIONAL = $refTranscriptSource + @indexing@ + #if $pairstraight.readselect == 'paired': - #if $pairstraight.file1.is_of_type("fastq.gz"): - && cp '${pairstraight.file1}' ./mate1.fastq.gz - && gunzip ./mate1.fastq.gz - && cp '${pairstraight.file2}' ./mate2.fastq.gz - && gunzip ./mate2.fastq.gz - #else if $pairstraight.file1.is_of_type("fastq.bz2"): - && cp '${pairstraight.file1}' ./mate1.fastq.bz2 - && bzip2 -d ./mate1.fastq.bz2 - && cp '${pairstraight.file2}' ./mate2.fastq.bz2 - && bzip2 -d ./mate2.fastq.bz2 + #if $pairstraight.mates.forward.is_of_type("fastq.gz"): + gunzip -c '${pairstraight.mates.forward}' > ./mate1.fastq && + gunzip -c '${pairstraight.mates.reverse}' > ./mate2.fastq && + #else if $pairstraight.mates.forward.is_of_type("fastq.bz2"): + bunzip2 -c '${pairstraight.mates.forward}' > ./mate1.fastq && + bunzip2 -c '${pairstraight.mates.reverse}' > ./mate2.fastq && #else: - && ln -s '${pairstraight.file1}' ./mate1.fastq - && ln -s '${pairstraight.file2}' ./mate2.fastq + ln -s '${pairstraight.mates.forward}' ./mate1.fastq && + ln -s '${pairstraight.mates.reverse}' ./mate2.fastq && #end if #else: #if $pairstraight.unmatedreads.is_of_type("fastq.gz"): - && cp '${pairstraight.unmatedreads}' ./unmate.fastq.gz - && gunzip ./unmate.fastq.gz + gunzip -c '${pairstraight.unmatedreads}' > ./unmate.fastq && #else if $pairstraight.unmatedreads.is_of_type("fastq.bz2"): - && cp '${pairstraight.unmatedreads}' ./unmate.fastq.bz2 - && bzip2 -d unmate.fastq.bz2 + bunzip2 -c '${pairstraight.unmatedreads}' > ./unmate.fastq && #else: - && ln -s '${pairstraight.unmatedreads}' ./unmate.fastq + ln -s '${pairstraight.unmatedreads}' ./unmate.fastq && #end if #end if - && ln -s '${tgmap}' ./alevinmap.tsv - && salmon alevin -l + ln -s '${tgmap}' ./alevinmap.tsv && + salmon --no-version-check alevin + --libType #if $pairstraight.readselect == 'paired': #if $pairstraight.libtype.strandedness == 'A' A #else ${pairstraight.libtype.orientation}${pairstraight.libtype.strandedness} #end if - -i $index_path - -1 ./mate1.fastq - -2 ./mate2.fastq + --index $index_path + --mates1 ./mate1.fastq + --mates2 ./mate2.fastq #else: '${pairstraight.libtype.strandedness}' - -i $index_path - -r zcat ./unmate.fastq + --index $index_path + --unmatedReads ./unmate.fastq #end if -o ./output - -p "\${GALAXY_SLOTS:-4}" + -p "\$SLOTS" ${protocol_cond.protocol} #if $protocol_cond.protocol == '--indropV2' --w1 '${protocol_cond.w1}' @@ -113,7 +105,7 @@ #end if #if 'dumpUmiGraph' in $output_files: && python '$__tool_directory__/vpolo_convert.py' -u - && sh '$__tool_directory__/umiout.sh' + && bash '$__tool_directory__/umiout.sh' #end if ## those gzip file include binary datasets ## #if $optional.numCellBootstraps: @@ -133,27 +125,26 @@ - - + - + - - - - - - - - - - + + + + + + + + + + @@ -173,7 +164,7 @@ - + @@ -255,10 +246,16 @@ - - - - + + + + + + + + + + @@ -274,7 +271,7 @@ - + @@ -287,10 +284,16 @@ - - - - + + + + + + + + + + @@ -307,10 +310,11 @@ - + + @@ -353,10 +363,16 @@ - - - - + + + + + + + + + + @@ -370,7 +386,7 @@ - + diff -r b0fe3e482cb7 -r 4dbb23a9fccc macros.xml --- a/macros.xml Mon Oct 28 17:50:35 2024 +0000 +++ b/macros.xml Wed Oct 29 10:47:09 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 ]]> - - 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