|
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>
|
|
1
|
18 <data name="calls" label="Call Sites (unfiltered SNPs)" format="vcf"/>
|
|
|
19 <data name="pileup" label="Call Sites (pileup)" format="pileup"/>
|
|
|
20 <data name="metrics" label="Call sites (Picard metrics)" format="txt"/>
|
|
0
|
21 </outputs>
|
|
|
22 <tests>
|
|
|
23 <test>
|
|
|
24 <param name="reference" value="lambda_virus.fasta"/>
|
|
|
25 <param name="input" value="reads.bam"/>
|
|
|
26 <output name="calls" ftype="vcf">
|
|
|
27 <assert_contents>
|
|
|
28 <has_size value="9976" delta="100"/>
|
|
1
|
29 </assert_contents>
|
|
0
|
30 </output>
|
|
|
31 <output name="pileup" ftype="pileup">
|
|
|
32 <assert_contents>
|
|
|
33 <has_size value="4259620" delta="100"/>
|
|
1
|
34 </assert_contents>
|
|
0
|
35 </output>
|
|
|
36 </test>
|
|
|
37 </tests>
|
|
|
38 <help><![CDATA[
|
|
|
39 <a href="http://snp-pipeline.readthedocs.io/en/latest/index.html">http://snp-pipeline.readthedocs.io/en/latest/index.html</a>
|
|
|
40 ]]></help>
|
|
|
41 <expand macro="citations"/>
|
|
|
42 </tool>
|
|
|
43
|
