Mercurial > repos > hepcat72 > lumpyexpress
diff lumpy_wrapper.tcsh @ 2:bfd745cf84c2 draft
planemo upload for repository https://github.com/hepcat72/robs_galaxy_tools/tree/master/tools/lumpyexpress commit 3e9afda4e4ad444b21dc30b866074ed6fa527bfa
author | hepcat72 |
---|---|
date | Thu, 25 Oct 2018 13:53:31 -0400 |
parents | e5150e64206a |
children | ea7edefc72c3 |
line wrap: on
line diff
--- a/lumpy_wrapper.tcsh Tue Oct 16 17:34:40 2018 -0400 +++ b/lumpy_wrapper.tcsh Thu Oct 25 13:53:31 2018 -0400 @@ -20,6 +20,22 @@ setenv DSCDS '' setenv SPLTOPTS '' +## +## The following is a work-around to a condition that causes a segfault in lumpy from lumpyexpress due to a new min_element threshold +## + +setenv LUMPYEXPRESSSCPT lumpyexpress +setenv LUMPYEXPRESSLOC `which lumpyexpress` +setenv LUMPYDIR `dirname $LUMPYEXPRESSLOC` +setenv LUMPYCONFIG lumpyexpress.config +setenv DISTROSCPT pairend_distro.py +setenv DISTROLOC `which $DISTROSCPT` +perl -e 'while(<STDIN>){s/min_elements = 10+/min_elements = 1/;print}' < $DISTROLOC > ./$DISTROSCPT +perl -e 'while(<STDIN>){s/PAIREND_DISTRO=.*/PAIREND_DISTRO=$ARGV[0]/;print}' ./$DISTROSCPT < $LUMPYDIR/$LUMPYCONFIG > ./$LUMPYCONFIG +perl -e 'while(<STDIN>){s/PAIREND_DISTRO=.*/PAIREND_DISTRO=$ARGV[0]/;print}' ./$DISTROSCPT < $LUMPYEXPRESSLOC > ./$LUMPYEXPRESSSCPT +chmod 555 ./$DISTROSCPT +chmod 555 ./$LUMPYEXPRESSSCPT + foreach b ( $MYARGV ) echo @@ -122,8 +138,8 @@ if ( $PAIREDENDRUN ) then echo - echo "Running: lumpyexpress -B $BAMS -S $SPLTS -D $DSCDS -o $OUTVCF" - lumpyexpress -B $BAMS -S $SPLTS -D $DSCDS -o $OUTVCF + echo "Running: ./$LUMPYEXPRESSSCPT -B $BAMS -S $SPLTS -D $DSCDS -o $OUTVCF" + ./$LUMPYEXPRESSSCPT -B $BAMS -S $SPLTS -D $DSCDS -o $OUTVCF if ( $status ) then echo "lumpyexpress failed"