|
0
|
1 <tool id="call_sites" name="CFSAN SNP pipeline: call sites" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="@PROFILE@">
|
|
|
2 <description>of high-confidence SNPs</description>
|
|
|
3 <macros>
|
|
|
4 <import>macros.xml</import>
|
|
|
5 </macros>
|
|
|
6 <expand macro="requirements"/>
|
|
|
7 <command detect_errors="exit_code"><![CDATA[
|
|
|
8 cp '$reference' 'reference.fasta' &&
|
|
|
9 samtools faidx 'reference.fasta' &&
|
|
|
10 picard MarkDuplicates INPUT='$input' OUTPUT=/dev/stdout METRICS_FILE='$metrics' | samtools mpileup -f ./reference.fasta - -o '$pileup' &&
|
|
|
11 varscan mpileup2snp '$pileup' --min_var-freq 0.90 --output-vcf 1 > '$calls'
|
|
|
12 ]]></command>
|
|
|
13 <inputs>
|
|
|
14 <param name="input" label="Read alignment to reference" type="data" format="bam"/>
|
|
|
15 <param name="reference" label="FASTA Reference from your history" type="data" format="fasta"/>
|
|
|
16 </inputs>
|
|
|
17 <outputs>
|
|
|
18 <data name="calls" label="${input.name.split(' ')[0]} unfiltered SNPs" format="vcf"/>
|
|
|
19 <data name="pileup" label="${input.name.split(' ')[0]} pileup" format="pileup"/>
|
|
|
20 <data name="metrics" label="Metrics from Picard" format="txt"/>
|
|
|
21 <data label="S3 Cache log" name="cache_log" format="txt"/>
|
|
|
22 </outputs>
|
|
|
23 <tests>
|
|
|
24 <test>
|
|
|
25 <param name="reference" value="lambda_virus.fasta"/>
|
|
|
26 <param name="input" value="reads.bam"/>
|
|
|
27 <output name="calls" ftype="vcf">
|
|
|
28 <assert_contents>
|
|
|
29 <has_size value="9976" delta="100"/>
|
|
|
30 </assert_contents>
|
|
|
31 </output>
|
|
|
32 <output name="pileup" ftype="pileup">
|
|
|
33 <assert_contents>
|
|
|
34 <has_size value="4259620" delta="100"/>
|
|
|
35 </assert_contents>
|
|
|
36 </output>
|
|
|
37 </test>
|
|
|
38 </tests>
|
|
|
39 <help><![CDATA[
|
|
|
40 <a href="http://snp-pipeline.readthedocs.io/en/latest/index.html">http://snp-pipeline.readthedocs.io/en/latest/index.html</a>
|
|
|
41 ]]></help>
|
|
|
42 <expand macro="citations"/>
|
|
|
43 </tool>
|
|
|
44
|
