|
jpayne@0
|
1 <tool id="call_sites" name="2. Call Sites" version="1.0.1" profile="16.10">
|
|
jpayne@0
|
2 <description>of high-confidence SNPs, or lookup cached sitecall</description>
|
|
jpayne@0
|
3 <requirements>
|
|
jpayne@61
|
4 <!-- <container type="docker">quay.io/biocontainers/mulled-v2-e518cb713b80020daa244e2ce559c37edb42dcde:5cb90165033bd74a97d0a807f439cf4b1aaabb89-0</container> -->
|
|
jpayne@61
|
5 <requirement type="package" version="2.4.5">bowtie2</requirement>
|
|
jpayne@50
|
6 <requirement type="package" version="1.9">samtools</requirement>
|
|
jpayne@58
|
7 <requirement type="package" version="1.0.8">bzip2</requirement>
|
|
jpayne@58
|
8 <requirement type="package" version="2.27.4">picard</requirement>
|
|
jpayne@58
|
9 <requirement type="package" version="2.4.4">varscan</requirement>
|
|
jpayne@61
|
10 <!-- <requirement type="package" version="1.9.134">boto3</requirement> -->
|
|
jpayne@58
|
11 <!-- <requirement type="package" version="3.6.8">python</requirement> -->
|
|
jpayne@0
|
12 </requirements>
|
|
jpayne@0
|
13 <command detect_errors="exit_code"><![CDATA[
|
|
jpayne@40
|
14 export LD_LIBRARY_PATH="\$CONDA_DEFAULT_ENV/lib" &&
|
|
jpayne@9
|
15 python $__tool_directory__/snp-cache.py snp_sitecalls
|
|
jpayne@0
|
16 "\$(md5sum $input $sample | cut -c -32 | md5sum | cut -c -32)"
|
|
jpayne@0
|
17 -c "
|
|
jpayne@0
|
18 cp $input ./reference.fasta &&
|
|
jpayne@0
|
19 samtools faidx ./reference.fasta &&
|
|
jpayne@0
|
20 samtools view -bS -F 4 $sample |
|
|
jpayne@0
|
21 samtools sort -o ./sorted.bam - &&
|
|
jpayne@0
|
22 picard MarkDuplicates INPUT=./sorted.bam
|
|
jpayne@0
|
23 OUTPUT=/dev/stdout
|
|
jpayne@0
|
24 METRICS_FILE=./metrics |
|
|
jpayne@0
|
25 samtools mpileup -f ./reference.fasta - -o ./pileup
|
|
jpayne@0
|
26 &&
|
|
jpayne@0
|
27 varscan mpileup2snp ./pileup
|
|
jpayne@0
|
28 --min_var-freq 0.90
|
|
jpayne@0
|
29 --output-vcf 1
|
|
jpayne@0
|
30 > ./calls
|
|
jpayne@0
|
31 && tar -zcf /dev/stdout ./metrics ./pileup ./calls
|
|
jpayne@0
|
32 "
|
|
jpayne@0
|
33 -o ./archive.tar.gz
|
|
jpayne@0
|
34 -l $cache_log
|
|
jpayne@0
|
35 && cat $cache_log
|
|
jpayne@0
|
36 && tar -zxvf ./archive.tar.gz
|
|
jpayne@0
|
37 && cp ./metrics $metrics
|
|
jpayne@0
|
38 && cp ./pileup $pileup
|
|
jpayne@0
|
39 && cp ./calls $calls
|
|
jpayne@0
|
40 ]]></command>
|
|
jpayne@0
|
41 <inputs>
|
|
jpayne@0
|
42 <param name="input" label="FASTA Reference from your history" type="data" format="fasta" />
|
|
jpayne@0
|
43 <!-- <conditional name="reads">
|
|
jpayne@0
|
44 <param name="reads_select" type="select" label="Paired-end collection, or two datasets from your history">
|
|
jpayne@0
|
45 <option value="collection">Paired collection from your history</option>
|
|
jpayne@0
|
46 <option value="history">Two FASTQ datasets from your history</option>
|
|
jpayne@0
|
47 </param>
|
|
jpayne@0
|
48 <when value="collection">
|
|
jpayne@0
|
49 <param name="coll" type="data_collection" format="fastq,fastqsanger" collection_type="paired">
|
|
jpayne@0
|
50 </when>
|
|
jpayne@0
|
51 <when value="history">
|
|
jpayne@0
|
52 <param type="data" name="forward" format="fastq,fastqsanger" />
|
|
jpayne@0
|
53 <param type="data" name="reverse" format="fastq,fastqsanger" />
|
|
jpayne@0
|
54 </when>
|
|
jpayne@0
|
55 </conditional> -->
|
|
jpayne@0
|
56 <param name="sample" label="Read alignment to reference" type="data" format="sam" />
|
|
jpayne@0
|
57 </inputs>
|
|
jpayne@0
|
58 <outputs>
|
|
jpayne@0
|
59 <data name="calls" label="${sample.name.split(' ')[0]} unfiltered SNPs" format="vcf" />
|
|
jpayne@0
|
60 <data name="pileup" label="${sample.name.split(' ')[0]} pileup" format="pileup" hidden="true" />
|
|
jpayne@0
|
61 <data name="metrics" label="Metrics from Picard" format="txt" hidden="true" />
|
|
jpayne@0
|
62 <data label="S3 Cache log" name="cache_log" format="txt" hidden="true" />
|
|
jpayne@0
|
63 </outputs>
|
|
jpayne@0
|
64 <tests>
|
|
jpayne@0
|
65 <test>
|
|
jpayne@0
|
66 <param name="input" value="reference/lambda_virus.fasta" />
|
|
jpayne@0
|
67 <param name="sample" value="samples/sample2/reads.sam" />
|
|
jpayne@0
|
68 <output name="calls" value="samples/sample2/var.flt.vcf" />
|
|
jpayne@0
|
69 <output name="pileup" value="samples/sample2/reads.all.pileup" />
|
|
jpayne@0
|
70 </test>
|
|
jpayne@0
|
71 </tests>
|
|
jpayne@0
|
72 <help><![CDATA[
|
|
jpayne@0
|
73 <a href="http://snp-pipeline.readthedocs.io/en/latest/index.html">http://snp-pipeline.readthedocs.io/en/latest/index.html</a>
|
|
jpayne@0
|
74 ]]></help>
|
|
jpayne@0
|
75 <citations>
|
|
jpayne@0
|
76 <citation type="doi">10.7717/peerj-cs.20</citation>
|
|
jpayne@0
|
77 </citations>
|
|
jpayne@0
|
78 </tool> |
