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