diff genap2_kb_python/kb_count.xml @ 0:e8d93f1429c2 draft

Uploaded

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/genap2_kb_python/kb_count.xml	Fri Apr 16 18:18:21 2021 +0000
@@ -0,0 +1,143 @@
+<tool id="kb_python" name="kb_python" version="@VERSION@+galaxy0">
+	<description>performs gene and feature quantification on single-cell sequencing data.</description>
+	<macros>
+        	<import>macros.xml</import>
+    </macros>
+	<requirements>
+		<requirement type="package" version="0.25.1">kb-python</requirement>
+	</requirements>
+	<stdio>
+		<exit_code range="1:" />
+		<exit_code range=":-1" />
+		<regex match="Error:" />
+		<regex match="Exception:" />
+		<regex match="Exception :" />
+	</stdio>
+	<command detect_errors="exit_code"><![CDATA[
+		mkdir ./index
+		&& mkdir ./kb_outs
+		#if $refTranscriptSource.TranscriptSource == "history":
+			&& ln -s '${refTranscriptSource.h_index.index_file}' './index/kb_ref.idx'
+			&& ln -s '${refTranscriptSource.h_index.t2g_file}' './index/t2g.txt'
+			#if $workflow == "lamanno" or $workflow == "nucleus":
+				&& ln -s '${refTranscriptSource.h_index.history_lamanno.cdna_t2c}' './index/cdna_t2c.txt'
+				&& ln -s '${refTranscriptSource.h_index.history_lamanno.intron_t2c}' './index/intron_t2c.txt'
+			#end if
+			#set $index_path = './index'
+		#else if $refTranscriptSource.TranscriptSource == "built":
+			&& kb ref -i ./index/kb_ref.idx
+				-g ./index/t2g.txt
+				-f1 ./index/cdna.fa	
+				--workflow $workflow
+				#if $workflow == "lamanno" or $workflow == "nucleus":
+					-f2 ./index/intron.fa
+					-c1 ./index/cdna_t2c.txt
+					-c2 ./index/intron_t2c.txt
+					--workflow $workflow
+				#end if
+				#if $workflow != "kite":
+					'${refTranscriptSource.s_index.genomic_fasta}'
+					'${refTranscriptSource.s_index.genomic_gtf}'
+				#else:
+					'${refTranscriptSource.s_index.kite_table}'
+				#end if
+				#set $index_path = './index'
+        #end if
+		&& kb count
+		-t \${GALAXY_SLOTS:-1}
+		-i $index_path/kb_ref.idx 
+		-g $index_path/t2g.txt
+		#if $workflow == "lamanno" or  $workflow == "nucleus":
+			-c1 $index_path/cdna_t2c.txt
+			-c2 $index_path/intron_t2c.txt
+		#end if
+		-x $technology
+        #if $whitelist:
+            --whitelist '${optional.whitelist}'
+        #end if
+		#if $optional.multimap:
+			--mm
+		#end if
+		--workflow $workflow
+		#if $extra_dtype != "none":
+			$extra_dtype
+		#end if
+		-o ./kb_outs 
+		--cellranger 
+		$FASTQ1
+		$FASTQ2
+		]]>
+	</command>
+	<inputs>
+		<param name="workflow" label="Select the workflow you want to use:" type="select" multiple="false" format="text" help="Type of workflow. Use `lamanno` for RNA velocity based on La Manno et al. 2018 logic.
+		Use `nucleus` for RNA velocity on single-nucleus RNA-seq reads. Use `kite` for feature barcoding. Use `kite:10xFB` for 10x Genomics Feature Barcoding technology.">
+			<option value="standard">standard</option>
+			<option value="lamanno">lamanno</option>
+			<option value="nucleus">nucleus</option>
+			<option value="kite">kite</option>
+		</param>
+		<expand macro="index"/>
+		<param name="technology" label="Select the scRNA-seq technology:" type="select" multiple="false" format="text" help="Choose the scRNA-seq technology used to generate the fastq data.">
+			<option value="10XV1">10XV1</option>
+			<option value="10XV2">10XV2</option>
+			<option value="10XV3">10XV3</option>
+			<option value="CELSEQ">CELSEQ</option>
+			<option value="CELSEQ2">CELSEQ2</option>
+			<option value="DROPSEQ">DROPSEQ</option>
+			<option value="INDROPSV1">INDROPSV1</option>
+			<option value="INDROPSV2">INDROPSV2</option>
+			<option value="INDROPSV3">INDROPSV3</option>
+			<option value="SCRUBSEQ">SCRUBSEQ</option>
+			<option value="SURECELL">SURECELL</option>
+		</param>
+		<param name="FASTQ1" label="Select the R1 fastq file:" type="data" format="fastqsanger.gz" multiple="false"/>
+		<param name="FASTQ2" label="Select the R2 fastq file:" type="data" format="fastqsanger.gz" multiple="false"/> 
+		<param name="extra_dtype" type="select" label="Do you want any additional output data type beside CellRanger?">
+			<option value="none" selected = "true">No</option>
+			<option value="--loom">Loom</option>
+			<option value="--h5ad">H5ad</option>
+		</param>
+		<section name="optional" title="Optional commands" expanded="false">
+			<param name="whitelist" type="data" format="tsv,tabular" optional="true" label="Whitelist file" help="Whitelisted barcodes to correct to. If not provided and bustools supports the technology, a pre-packaged whitelist is used. If not, the bustools
+			whitelist command is used."/>
+			<param name="multimap" type="boolean" optional="true" label="Include multi pseudoaligned reads?" help="Do you want to include reads that pseudoalign to multiple genes?"/>
+		</section>
+	</inputs>
+    <outputs>
+		<data name="barcodes" label="cellranger barcodes" format="txt" from_work_dir="kb_outs/counts_unfiltered/cellranger/barcodes.tsv"/>
+		<data name="genes" label="cellranger genes" format="txt" from_work_dir="kb_outs/counts_unfiltered/cellranger/genes.tsv"/>
+		<data name="matrix" label="cellranger matrix" format="mtx" from_work_dir="kb_outs/counts_unfiltered/cellranger/matrix.mtx" />
+		<data name="inspect" label="inspect_report" format="json" from_work_dir="kb_outs/inspect.json"/>
+		<data name="runinfo" label="run info" format="json" from_work_dir="kb_outs/run_info.json"/>
+		<data name="kbinfo" label="kb info" format="json" from_work_dir="kb_outs/kb_info.json"/>
+		<data name="kb_ref" label="kb_ref.idx" format="kallisto.idx" from_work_dir="index/kb_ref.idx"/>
+		<data name="t2g" label="t2g.txt" format="txt" from_work_dir="index/t2g.txt"/>
+		<data name="cdna_t2c" label="-c1 cdna_t2c.txt" format="txt" from_work_dir="index/cdna_t2c.txt">
+			<filter>workflow == "lamanno"</filter>
+		</data>
+		<data name="intron_t2c" label="-c2 intron_t2c.txt" format="txt" from_work_dir="index/intron_t2c.txt">
+			<filter>workflow == "lamanno"</filter>
+		</data>	
+		<data name="adata" label="adata.h5ad" format="h5ad" from_work_dir="kb_outs/counts_unfiltered/adata.h5ad">
+			<filter	>extra_dtype == "--h5ad"</filter>
+		</data>
+		<data name="adata" label="adata.loom" format="loom" from_work_dir="kb_outs/counts_unfiltered/adata.loom">
+			<filter>extra_dtype == "--loom"</filter>
+		</data>
+	</outputs>
+	<tests>
+		<expand macro="tests"/>
+	</tests>
+	<help><![CDATA[
+		This is a galaxy wrapper for kallisto-bustools (kb) for quantification of different types of single-cell sequencing data. The main kallisto-bustools homepage can be found under: 
+		'kallisto-bustools<https://www.kallistobus.tools/>'.
+
+		The kb count command runs the kallisto and bustools programs. It can be used for pre-processing of data from a variety of single-cell RNA-seq technologies, 
+		and for a number of different workflows (e.g. production of gene count matrices, RNA velocity analyses, etc.).
+  
+		The kb package is developed and maintained by the Pachterlab under the MIT License.
+	  ]]></help>
+	<citations>
+        <citation type="doi">https://doi.org/10.1101/673285</citation>
+    </citations>
+</tool>
\ No newline at end of file