# HG changeset patch
# User mini
# Date 1417080718 -3600
# Node ID 1c8dcda28be76968016b175adc32ac76cf669fb0
# Parent f785839fb39544a28b31bcd8ee39cff3d40de7ee
version 27/11/2014, corrected extra arguments bug
diff -r f785839fb395 -r 1c8dcda28be7 macro_configuration.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macro_configuration.xml Thu Nov 27 10:31:58 2014 +0100
@@ -0,0 +1,23 @@
+
+ macro_customConfiguration.xml
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
\ No newline at end of file
diff -r f785839fb395 -r 1c8dcda28be7 macro_customConfiguration.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macro_customConfiguration.xml Thu Nov 27 10:31:58 2014 +0100
@@ -0,0 +1,27 @@
+
+ macro_extraStrelkaArguments.xml
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
\ No newline at end of file
diff -r f785839fb395 -r 1c8dcda28be7 macro_extraStrelkaArguments.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macro_extraStrelkaArguments.xml Thu Nov 27 10:31:58 2014 +0100
@@ -0,0 +1,30 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff -r f785839fb395 -r 1c8dcda28be7 macro_help.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macro_help.xml Thu Nov 27 10:31:58 2014 +0100
@@ -0,0 +1,56 @@
+
+
+
+ Strelka, a method for somatic SNV and small indel detection from sequencing data of matched tumor-normal samples. You can see more information at : https://sites.google.com/site/strelkasomaticvariantcaller.
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff -r f785839fb395 -r 1c8dcda28be7 macro_inputs.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macro_inputs.xml Thu Nov 27 10:31:58 2014 +0100
@@ -0,0 +1,33 @@
+
+ macro_configuration.xml
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
\ No newline at end of file
diff -r f785839fb395 -r 1c8dcda28be7 macro_main.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macro_main.xml Thu Nov 27 10:31:58 2014 +0100
@@ -0,0 +1,76 @@
+
+
+ Strelka
+
+ SCRIPT_PATH_STRELKA
+ samtools
+ vcftools
+
+ strelka_wrapper.py --tumorBam $tumorBam --normalBam $normalBam
+ #silent sys.stderr.write("!!!! Cheetah Template Variables !!!!\n")
+ #for k,v in $searchList[2].items()
+ #silent sys.stderr.write(" %s = %s\n" % (str(k), str(v) ))
+ #end for
+ #silent sys.stderr.write("!!!! end-of-list !!!!\n")
+ #import string
+
+ --scriptPath \$SCRIPT_PATH_STRELKA
+ #if $genomeSource.refGenomeSource == "history":
+ --refFile "${genomeSource.ownFile}"
+ #else:
+ --refFile "${genomeSource.index.fields.path}"
+ #end if
+
+ #if $configuration.configuration_switch == 'Default for Bwa':
+ --configFile strelka_config_bwa_default.ini
+ #else if $configuration.configuration_switch == 'Default for Isaac':
+ --configFile strelka_config_isaac_default.ini
+ #else if $configuration.configuration_switch == 'Default for Eland':
+ --configFile strelka_config_eland_default.ini
+ #else if $configuration.configuration_switch == 'Path':
+ --configFile $configuration.configFile
+ #else:
+ --configFile Custom
+ --depthFilterMultiple $configuration.depthFilterMultiple
+ --snvMaxFilteredBasecallFrac $configuration.snvMaxFilteredBasecallFrac
+ --snvMaxSpanningDeletionFrac $configuration.snvMaxSpanningDeletionFrac
+ --indelMaxRefRepeat $configuration.indelMaxRefRepeat
+ --indelMaxWindowFilteredBasecallFrac $configuration.indelMaxWindowFilteredBasecallFrac
+ --indelMaxIntHpolLength $configuration.indelMaxIntHpolLength
+ --ssnvPrior $configuration.ssnvPrior
+ --sindelPrior $configuration.sindelPrior
+ --ssnvNoise $configuration.ssnvNoise
+ --sindelNoise $configuration.sindelNoise
+ --ssnvNoiseStrandBiasFrac $configuration.ssnvNoiseStrandBiasFrac
+ --minTier1Mapq $configuration.minTier1Mapq
+ --minTier2Mapq $configuration.minTier2Mapq
+ --ssnvQuality_LowerBound $configuration.ssnvQuality_LowerBound
+ --sindelQuality_LowerBound $configuration.sindelQuality_LowerBound
+ --isWriteRealignedBam $configuration.isWriteRealignedBam
+ --binSize $configuration.binSize
+ --isSkipDepthFilters $configuration.isSkipDepthFilters
+ --maxInputDepth $configuration.maxInputDepth
+ #if $configuration.extra_arguments.extra_arguments_switch == 'Yes':
+ --extraStrelkaArguments yes
+ #if $configuration.extra_arguments.a:
+ -a
+ #end if
+ #if $configuration.extra_arguments.b:
+ -b $configuration.extra_arguments.barg
+ #end if
+ #if $configuration.extra_arguments.c:
+ -c $configuration.extra_arguments.carg
+ #end if
+ #if $configuration.extra_arguments.d:
+ -d $configuration.extra_arguments.darg
+ #end if
+ #if $configuration.extra_arguments.e:
+ -e $configuration.extra_arguments.earg
+ #end if
+ #end if
+ #end if
+
+
+
+
+
diff -r f785839fb395 -r 1c8dcda28be7 macro_outputs.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macro_outputs.xml Thu Nov 27 10:31:58 2014 +0100
@@ -0,0 +1,27 @@
+
+
+
+
+ output1_file_check == True
+
+
+ output2_file_check == True
+
+
+ output3_file_check == True
+
+
+ output4_file_check == True
+
+
+ conf_file_check == True
+
+
+ output5_file_check == True
+
+
+ output6_file_check == True
+
+
+
+
diff -r f785839fb395 -r 1c8dcda28be7 macro_tests.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macro_tests.xml Thu Nov 27 10:31:58 2014 +0100
@@ -0,0 +1,16 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff -r f785839fb395 -r 1c8dcda28be7 strelka.xml
--- a/strelka.xml Wed Oct 15 15:47:03 2014 +0200
+++ b/strelka.xml Thu Nov 27 10:31:58 2014 +0100
@@ -1,172 +1,24 @@
-
-
- Strelka
-
- SCRIPT_PATH_STRELKA
- samtools
- vcftools
-
- strelka_wrapper.py --tumorBam $tumorBam --normalBam $normalBam
- #if $genomeSource.refGenomeSource == "history":
- --refFile "${genomeSource.ownFile}"
- #else:
- --refFile "${genomeSource.index.fields.path}"
- #end if
+
- #if $configuration.configuration_switch == 'Default for Bwa':
- --configFile strelka_config_bwa_default.ini
- #else if $configuration.configuration_switch == 'Default for Isaac':
- --configFile strelka_config_isaac_default.ini
- #else if $configuration.configuration_switch == 'Default for Eland':
- --configFile strelka_config_eland_default.ini
- #else if $configuration.configuration_switch == 'Path':
- --configFile $configuration.configFile
- #else:
- --configFile Custom
- --depthFilterMultiple $configuration.depthFilterMultiple
- --snvMaxFilteredBasecallFrac $configuration.snvMaxFilteredBasecallFrac
- --snvMaxSpanningDeletionFrac $configuration.snvMaxSpanningDeletionFrac
- --indelMaxRefRepeat $configuration.indelMaxRefRepeat
- --indelMaxWindowFilteredBasecallFrac $configuration.indelMaxWindowFilteredBasecallFrac
- --indelMaxIntHpolLength $configuration.indelMaxIntHpolLength
- --ssnvPrior $configuration.ssnvPrior
- --sindelPrior $configuration.sindelPrior
- --ssnvNoise $configuration.ssnvNoise
- --sindelNoise $configuration.sindelNoise
- --ssnvNoiseStrandBiasFrac $configuration.ssnvNoiseStrandBiasFrac
- --minTier1Mapq $configuration.minTier1Mapq
- --minTier2Mapq $configuration.minTier2Mapq
- --ssnvQuality_LowerBound $configuration.ssnvQuality_LowerBound
- --sindelQuality_LowerBound $configuration.sindelQuality_LowerBound
- --isWriteRealignedBam $configuration.isWriteRealignedBam
- --binSize $configuration.binSize
- --isSkipDepthFilters $configuration.isSkipDepthFilters
- --maxInputDepth $configuration.maxInputDepth
- #if $configuration.extra_arguments.extra_arguments_switch == 'Yes':
- --extraStrelkaArguments $configuration.extra_arguments.extraStrelkaArguments
- #end if
- #end if
- --scriptPath \$SCRIPT_PATH_STRELKA
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- output1_file_check == True
-
-
- output2_file_check == True
-
-
- output3_file_check == True
-
-
- output4_file_check == True
-
-
- conf_file_check == True
-
-
- output5_file_check == True
-
-
- output6_file_check == True
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-Strelka, a method for somatic SNV and small indel detection from sequencing data of matched tumor-normal samples.
-You can see more information at : https://sites.google.com/site/strelkasomaticvariantcaller.
-
-
+
+
+
+
+ macro_main.xml
+ macro_inputs.xml
+
+ macro_outputs.xml
+ macro_tests.xml
+ macro_help.xml
+
+
+
+
+
+
+
+
diff -r f785839fb395 -r 1c8dcda28be7 strelka_wrapper.py
--- a/strelka_wrapper.py Wed Oct 15 15:47:03 2014 +0200
+++ b/strelka_wrapper.py Thu Nov 27 10:31:58 2014 +0100
@@ -20,9 +20,23 @@
def _create_config(args, config_path):
conf_file = open(config_path, "w")
conf_file.write("[user]\n")
- for option in args:
- if not option in ["tumorBam", "normalBam", "refFile", "configFile", "scriptPath"] and args[option]!=None:
- conf_file.write("%s=%s\n" % (option, args[option]))
+ args2 = vars(args)
+ for option in args2:
+ if not option in ["tumorBam", "normalBam", "refFile", "configFile", "scriptPath", "a", "b", "c", "d", "e", "extraStrelkaArguments"] and args2[option]!=None:
+ conf_file.write("%s=%s\n" % (option, args2[option]))
+ if args.extraStrelkaArguments == "yes":
+ conf_file.write("extraStrelkaArguments=")
+ if args.a:
+ conf_file.write("--ignore-conflicting-read-names ")
+ if args.b != None:
+ conf_file.write("-used-allele-count-min-qscore %s " % (args.b))
+ if args.c != None:
+ conf_file.write("--candidate-indel-input-vcf %s " % (args.c))
+ if args.d != None:
+ conf_file.write("--force-output-vcf %s " % (args.d))
+ if args.e != None:
+ conf_file.write("-min-small-candidate-indel-read-frac %s " % (args.e))
+ conf_file.write("\n")
conf_file.close()
def my_Popen(cmd, prefix_for_stderr_name, tmp_dir, msg_error):
@@ -60,10 +74,10 @@
#Manage options
print(os.environ['PATH'])
parser = argparse.ArgumentParser()
- parser.add_argument( '-t', '--tumorBam', help='path to tumor bam file', required = False )
- parser.add_argument( '-n', '--normalBam', help='', required = False )
- parser.add_argument( '-r', '--refFile', help='', required = False )
- parser.add_argument( '-c', '--configFile', help='', required = False )
+ parser.add_argument( '--tumorBam', help='path to tumor bam file', required = False )
+ parser.add_argument( '--normalBam', help='', required = False )
+ parser.add_argument( '--refFile', help='', required = False )
+ parser.add_argument( '--configFile', help='', required = False )
parser.add_argument( '--depthFilterMultiple', help='', required = False )
parser.add_argument( '--snvMaxFilteredBasecallFrac', help='', required = False )
parser.add_argument( '--snvMaxSpanningDeletionFrac', help='', required = False )
@@ -85,6 +99,11 @@
parser.add_argument( '--isSkipDepthFilters', help='', required = False )
parser.add_argument( '--maxInputDepth', help='', required = False )
parser.add_argument( '--scriptPath', help='', required = False )
+ parser.add_argument( '-a', action="store_true", help='', required = False )
+ parser.add_argument( '-b', help='', required = False )
+ parser.add_argument( '-c', help='', required = False )
+ parser.add_argument( '-d', help='', required = False )
+ parser.add_argument( '-e', help='', required = False )
args = parser.parse_args()
root_dir= args.scriptPath
@@ -114,7 +133,7 @@
#creating config file if needed
if args.configFile == "Custom":
- _create_config(vars(args), config_ini)
+ _create_config(args, config_ini)
elif args.configFile in ["strelka_config_bwa_default.ini", "strelka_config_isaac_default.ini", "strelka_config_eland_default.ini"]:
cmdbash="cp %s %s" % (root_dir + "/lib/" + args.configFile, config_ini)
my_Popen(cmdbash, "copy_default_file_err", tmp_dir, "Error during the copy of default config file, maybe it was removed")
diff -r f785839fb395 -r 1c8dcda28be7 strelkarecovery.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/strelkarecovery.xml Thu Nov 27 10:31:58 2014 +0100
@@ -0,0 +1,208 @@
+
+
+ Strelka
+
+ SCRIPT_PATH_STRELKA
+ samtools
+ vcftools
+
+ strelka_wrapper.py --tumorBam $tumorBam --normalBam $normalBam
+ --scriptPath \$SCRIPT_PATH_STRELKA
+ #if $genomeSource.refGenomeSource == "history":
+ --refFile "${genomeSource.ownFile}"
+ #else:
+ --refFile "${genomeSource.index.fields.path}"
+ #end if
+
+ #if $configuration.configuration_switch == 'Default for Bwa':
+ --configFile strelka_config_bwa_default.ini
+ #else if $configuration.configuration_switch == 'Default for Isaac':
+ --configFile strelka_config_isaac_default.ini
+ #else if $configuration.configuration_switch == 'Default for Eland':
+ --configFile strelka_config_eland_default.ini
+ #else if $configuration.configuration_switch == 'Path':
+ --configFile $configuration.configFile
+ #else:
+ --configFile Custom
+ --depthFilterMultiple $configuration.depthFilterMultiple
+ --snvMaxFilteredBasecallFrac $configuration.snvMaxFilteredBasecallFrac
+ --snvMaxSpanningDeletionFrac $configuration.snvMaxSpanningDeletionFrac
+ --indelMaxRefRepeat $configuration.indelMaxRefRepeat
+ --indelMaxWindowFilteredBasecallFrac $configuration.indelMaxWindowFilteredBasecallFrac
+ --indelMaxIntHpolLength $configuration.indelMaxIntHpolLength
+ --ssnvPrior $configuration.ssnvPrior
+ --sindelPrior $configuration.sindelPrior
+ --ssnvNoise $configuration.ssnvNoise
+ --sindelNoise $configuration.sindelNoise
+ --ssnvNoiseStrandBiasFrac $configuration.ssnvNoiseStrandBiasFrac
+ --minTier1Mapq $configuration.minTier1Mapq
+ --minTier2Mapq $configuration.minTier2Mapq
+ --ssnvQuality_LowerBound $configuration.ssnvQuality_LowerBound
+ --sindelQuality_LowerBound $configuration.sindelQuality_LowerBound
+ --isWriteRealignedBam $configuration.isWriteRealignedBam
+ --binSize $configuration.binSize
+ --isSkipDepthFilters $configuration.isSkipDepthFilters
+ --maxInputDepth $configuration.maxInputDepth
+ #if $configuration.extra_arguments.extra_arguments_switch == 'Yes':
+ --extraStrelkaArguments yes
+ #if $configuration.extra_arguments.a:
+ -a
+ #end if
+ #if 'b' in $configuration.extra_arguments.extraArgumentsChoice.extra_choice:
+ -b $configuration.extra_arguments.extraArgumentsChoice.b
+ #end if
+ #if 'c' in $configuration.extra_arguments.extraArgumentsChoice.extra_choice:
+ -c $configuration.extra_arguments.extraArgumentsChoice.c
+ #end if
+ #if 'd' in $configuration.extra_arguments.extraArgumentsChoice.extra_choice:
+ -d $configuration.extra_arguments.extraArgumentsChoice.d
+ #end if
+ #if 'e' in $configuration.extra_arguments.extraArgumentsChoice.extra_choice:
+ -e $configuration.extra_arguments.extraArgumentsChoice.e
+ #end if
+ #end if
+ #end if
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+ output1_file_check == True
+
+
+ output2_file_check == True
+
+
+ output3_file_check == True
+
+
+ output4_file_check == True
+
+
+ conf_file_check == True
+
+
+ output5_file_check == True
+
+
+ output6_file_check == True
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+Strelka, a method for somatic SNV and small indel detection from sequencing data of matched tumor-normal samples.
+You can see more information at : https://sites.google.com/site/strelkasomaticvariantcaller.
+
+
+