Mercurial > repos > drosofff > lumpy

comparison lumpy.xml @ 10:3856d59b4f2a draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/lumpy commit a7bd16b6c6dbfe8b41de2a51c29b197ca660c5c6
author drosofff
date Mon, 12 Dec 2016 08:48:47 -0500
parents ee2ee0b9de52
children 4ae2b4c68862
comparison
equal deleted inserted replaced
9:ee2ee0b9de52 10:3856d59b4f2a
23 samtools view -h input.bam | python $__tool_directory__/extractSplitReads_BwaMem.py -i stdin | samtools view -Sb - > "input.splitters.unsorted.bam" && 23 samtools view -h input.bam | python $__tool_directory__/extractSplitReads_BwaMem.py -i stdin | samtools view -Sb - > "input.splitters.unsorted.bam" &&
24 samtools sort input.discordants.unsorted.bam > input.discordants.bam && 24 samtools sort input.discordants.unsorted.bam > input.discordants.bam &&
25 samtools sort input.splitters.unsorted.bam > input.splitters.bam && 25 samtools sort input.splitters.unsorted.bam > input.splitters.bam &&
26 samtools view input.bam 26 samtools view input.bam
27 |tail -n +1 27 |tail -n +1
28 |python $__tool_directory__/pairend_distro.py -r $$analysis_type.readLength -X 4 -N $seq_method.additional_params.samplingValue -o input.lib.histo > meandev.txt && 28 |python $__tool_directory__/pairend_distro.py -r $analysis_type.readLength -X 4 -N $seq_method.additional_params.samplingValue -o input.lib.histo > meandev.txt &&
29 mean=\$(cat meandev.txt | sed s/mean:// | sed -r s/stdev:.+//) && 29 mean=\$(cat meandev.txt | sed s/mean:// | sed -r s/stdev:.+//) &&
30 stdev=\$(cat meandev.txt | sed -r s/mean:.+stdev://) && 30 stdev=\$(cat meandev.txt | sed -r s/mean:.+stdev://) &&
31 lumpy -mw $seq_method.additional_params.mw -tt $seq_method.additional_params.tt 31 lumpy -mw $seq_method.additional_params.mw -tt $seq_method.additional_params.tt
32 #if $output_format == "BEDPE": 32 #if $output_format == "BEDPE":
33 -b 33 -b
62 samtools sort input.B.discordants.unsorted.bam > input.B.discordants.bam && 62 samtools sort input.B.discordants.unsorted.bam > input.B.discordants.bam &&
63 samtools sort input.A.splitters.unsorted.bam > input.A.splitters.bam && 63 samtools sort input.A.splitters.unsorted.bam > input.A.splitters.bam &&
64 samtools sort input.B.splitters.unsorted.bam > input.B.splitters.bam && 64 samtools sort input.B.splitters.unsorted.bam > input.B.splitters.bam &&
65 samtools view input.A.bam 65 samtools view input.A.bam
66 |tail -n +1 66 |tail -n +1
67 |python $__tool_directory__/pairend_distro.py -r $$analysis_type.readLength -X 4 -N $seq_method.additional_params.samplingValue -o input.A.lib.histo > meandevA.txt && 67 |python $__tool_directory__/pairend_distro.py -r $analysis_type.readLength -X 4 -N $seq_method.additional_params.samplingValue -o input.A.lib.histo > meandevA.txt &&
68 samtools view input.B.bam 68 samtools view input.B.bam
69 |tail -n +1 69 |tail -n +1
70 |python $__tool_directory__/pairend_distro.py -r $$analysis_type.readLengthB -X 4 -N $seq_method.additional_params.samplingValue -o input.B.lib.histo > meandevB.txt && 70 |python $__tool_directory__/pairend_distro.py -r $analysis_type.readLengthB -X 4 -N $seq_method.additional_params.samplingValue -o input.B.lib.histo > meandevB.txt &&
71 meanA=\$(cat meandevA.txt | sed s/mean:// | sed -r s/stdev:.+//) && 71 meanA=\$(cat meandevA.txt | sed s/mean:// | sed -r s/stdev:.+//) &&
72 meanB=\$(cat meandevB.txt | sed s/mean:// | sed -r s/stdev:.+//) && 72 meanB=\$(cat meandevB.txt | sed s/mean:// | sed -r s/stdev:.+//) &&
73 stdevA=\$(cat meandevA.txt | sed -r s/mean:.+stdev://) && 73 stdevA=\$(cat meandevA.txt | sed -r s/mean:.+stdev://) &&
74 stdevB=\$(cat meandevB.txt | sed -r s/mean:.+stdev://) && 74 stdevB=\$(cat meandevB.txt | sed -r s/mean:.+stdev://) &&
75 <<<<<<< HEAD 75 <<<<<<< HEAD