comparison 5_call_consensus.xml @ 0:eefdd97a6749

planemo upload commit b'7f6183b769772449fbcee903686b8d5ec5b7439f\n'-dirty
author jpayne
date Wed, 24 Jan 2018 14:18:21 -0500
parents
children 62191abf5426
comparison
equal deleted inserted replaced
-1:000000000000 0:eefdd97a6749
1 <tool id="call_consensus" name="5. Call consensus SNPs" version="1.0.1" profile="16.10">
2 <description>for a sample, or lookup cached consensus SNP calls</description>
3 <requirements>
4 <requirement type="package" version="1.0.1">snp-pipeline</requirement>
5 <requirement type="package" version="1.5.11">boto3</requirement>
6 </requirements>
7 <command detect_errors="exit_code"><![CDATA[
8 #set $identifier = $call_pileup.name.split(' ')[0]
9 $__tool_directory__/snp-cache.py snp_consenus_calls
10 "\$(md5sum $snps $calls $call_pileup | cut -c -32 | md5sum | cut -c -32)"
11 -c "
12 mkdir ./$identifier/
13 && cp $call_pileup ./$identifier/reads.all.pileup
14 && cfsan_snp_pipeline call_consensus -v 4 -l $snps --vcfFileName "consensus.vcf" ./$identifier/reads.all.pileup 1>&2 && tar zcv ./consensus.vcf ./consensus.fasta
15 "
16 -l $cache_log
17 | tar zvx ./
18 && cp ./consensus.fasta $calls_consensus_fa
19 && cp ./consensus.vcf $calls_consensus_vcf
20 && cat $cache_log
21 ]]></command>
22 <inputs>
23 <param type="data" name="snps" format="txt" />
24 <param type="data" name="calls" format="vcf" />
25 <param type="data" name="call_pileup" format="pileup" />
26 </inputs>
27 <outputs>
28 <data label="${calls.name.split(' ')[0]} consensus SNPs (VCF)" name="calls_consensus_vcf" format="vcf" />
29 <data label="${calls.name.split(' ')[0]} consensus SNPs (FASTA)" name="calls_consensus_fa" format="fasta" hidden="true"/>
30 <data label="S3 Cache log" name="cache_log" format="txt" hidden="true" />
31 </outputs>
32 <tests>
33 <test>
34 <param name="snps" value="snplist.txt" />
35 <param name="calls" value="samples/sample1/var.flt.vcf" />
36 <param name="call_pileup" value="samples/sample1/reads.all.pileup" />
37 <output name="calls_consensus_vcf" value="samples/sample1/consensus.vcf" lines_diff="6" />
38 <output name="calls_consensus_fa" value="samples/sample1/consensus.fasta" lines_diff="2" />
39 </test>
40 </tests>
41 <help><![CDATA[
42 <a href="http://snp-pipeline.readthedocs.io/en/latest/index.html">http://snp-pipeline.readthedocs.io/en/latest/index.html</a>
43 ]]></help>
44 <citations>
45 <citation type="doi">10.7717/peerj-cs.20</citation>
46 </citations>
47 </tool>