annotate varamplicnv-5bafb1c69d03/VarAmpliCNV_Count.xml @ 0:febc6023d37b draft

Uploaded
author geert-vandeweyer
date Fri, 25 Sep 2020 08:29:36 +0000
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
1 <tool id="VarAmpliCNV_Count" name="Count Aligned Amplicons" version="0.1.0" >
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
2 <requirements>
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
3 <container type="docker">cmgantwerpen/varamplicnv:1.0.0</container>
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
4 </requirements>
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
5 <command detect_errors="exit_code"><![CDATA[
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
6 ##prepare input and indices
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
7 ln -s '$bam' infile &&
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
8 #if $bam.is_of_type('bam'):
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
9 #if str( $bam.metadata.bam_index ) != "None":
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
10 ln -s '${bam.metadata.bam_index}' infile.bai &&
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
11 #else:
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
12 samtools index infile infile.bai &&
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
13 #end if
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
14 #elif $bam.is_of_type('cram'):
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
15 #if str( $bam.metadata.cram_index ) != "None":
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
16 ln -s '${bam.metadata.cram_index}' infile.crai &&
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
17 #else:
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
18 samtools index infile infile.crai &&
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
19 #end if
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
20 #end if
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
21
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
22 ## Run Counter.
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
23 parseBAM.py '$bed' 'infile' '$counts' '$unmapped' '$stats'
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
24
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
25 ]]></command>
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
26 <inputs>
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
27 <param name="bed" type="data" format="BED" optional="false" label="Unique Amplicons (BED)" help="This file is created by the VarAmpliCNV 'Annotate Amplicons' tool." />
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
28 <param name="bam" type="data" format="bam" optional="false" label="BAM file" help="Sequencing data file" />
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
29 </inputs>
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
30 <outputs>
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
31 <data name="counts" format="tabular" label="VarAmpliCNV on ${on_string}: Counts" />
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
32 <data name="unmapped" format="bam" label="VarAmpliCNV on ${on_string}: UnMapped Reads" />
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
33 <data name="stats" format="txt" label="VarAmpliCNV on ${on_string}: Statistics" />
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
34 </outputs>
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
35 <help>
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
36 **VarAmpliCNV : Counting**
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
37
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
38 BAM files are parsed for readpairs exactly matching specified amplicons, based on start and end position.
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
39
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
40 **Parameters are :**
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
41
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
42 * Amplicon Design File (BED) : The de-duplicated amplicon list, generated by "varAmpliCNV Annotate".
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
43 * Sample Data (BAM) : The sample read data, provided as a single BAM file, or a collection of BAM files.
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
44
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
45 **Output files :**
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
46
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
47 * Count file (txt): the amplicon-read table.
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
48 * Unmapped (BAM) : Reads not matching amplicons. This bam file can be used to investigate issues.
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
49 * Stats (txt): Overview of the matching performance. Use it to investigate issues.
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
50
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
51 In case a collection of BAM files is provided, the output files will be grouped in collections as well.
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
52
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
53 </help>
febc6023d37b Uploaded
geert-vandeweyer
parents:
diff changeset
54 </tool>