Mercurial > repos > portiahollyoak > breakdancer
diff create_config_file.py @ 0:694cbef83948 draft
planemo upload for repository https://github.com/portiahollyoak/Tools commit 9700cff725be112c4a813f914eca946e158fbebd-dirty
author | portiahollyoak |
---|---|
date | Fri, 13 May 2016 11:47:32 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/create_config_file.py Fri May 13 11:47:32 2016 -0400 @@ -0,0 +1,21 @@ +import argparse + + +description = ("This script will create a configuration file for samples to be run in Breakdancer." + "If pooled analysis desired, only one config file needed for all samples." + "Otherwise, individual analysis requires individual config files for each sample.") + +parser = argparse.ArgumentParser() +parser.add_argument("--input_file", nargs="*", help="One or more alignment files") +parser.add_argument("--mean", nargs="+", help="Mean insert size") +parser.add_argument("--std_dev", nargs="+", help="The standard deviation of insert size") +parser.add_argument("--read_length", nargs="+", help="Average read length") +parser.add_argument("--sample_name", nargs="+", help="Sample name") +parser.add_argument("--output_config_file", help="Name of the output config file") +args = parser.parse_args() + +template = "map:{input_file}\tmean:{mean}\tstd:{std_dev}\treadlen:{read_length}\tsample:{sample_name}\texe:samtools view\n" +with open(args.output_config_file, "w") as output: + for input_file, mean, std_dev, read_length, sample_name in zip(args.input_file, args.mean, args.std_dev, args.read_length, args.sample_name): + config_line = template.format(input_file=input_file, mean=mean, std_dev=std_dev, read_length=read_length, sample_name=sample_name) + output.write(config_line) \ No newline at end of file