|
46 | 46 | Somatic variant calling pipeline
|
47 | 47 | Pipeline version: $version
|
48 | 48 |
|
49 |
| -<<<<<<< HEAD |
50 |
| -$yellow Usage: perl $0 --srg --step --sre --rdir --ref --log --q --mincovt --mincovn --minvaf --maxindsize --exonic --wgs --smg |
51 |
| -======= |
52 | 49 | $yellow Usage: perl $0 --srg --sre --wgs --rdir --ref --log --q --mincovt --mincovn --minvaf --maxindsize --exonic --smg --groupname --users --step
|
53 |
| ->>>>>>> a4a79c27f83ba1ebce51a18020e77b941eb33d2b |
54 | 50 |
|
55 | 51 | $normal
|
56 | 52 |
|
|
60 | 56 | <groupname> = job group name
|
61 | 57 | <users> = user name for job group
|
62 | 58 | <srg> = bam having read group or not: 1, yes and 0, no (default 1)
|
63 |
| -<<<<<<< HEAD |
64 |
| -<sre> = re-run: 1, yes and 0, no (default 0) |
65 |
| -<wgs> = 1 if it is wgs data and otherwise it is 0; If you want to output the maf for all variants, set exonic to 0 |
66 |
| -======= |
67 | 59 | <sre> = re-run and overwrite previous results: 1, yes and 0, no (default 0)
|
68 | 60 | >>>>>>> a4a79c27f83ba1ebce51a18020e77b941eb33d2b
|
69 | 61 | <step> run this pipeline step by step. (user must provide)
|
|
123 | 115 | my $status = &GetOptions (
|
124 | 116 | "step=i" => \$step_number,
|
125 | 117 | "srg=i" => \$status_rg,
|
126 |
| -<<<<<<< HEAD |
127 | 118 | "sre=i" => \$status_rerun,
|
128 |
| - "wgs=i" => \$s_wgs, |
129 |
| -======= |
130 |
| - "wgs=i" => \$s_wgs, |
131 |
| - "sre=i" => \$status_rerun, |
132 |
| - "groupname=s" => \$group_name, |
| 119 | + "wgs=i" => \$s_wgs, |
| 120 | + "groupname=s" => \$group_name, |
133 | 121 | "users=s" => \$compute_username,
|
134 |
| ->>>>>>> a4a79c27f83ba1ebce51a18020e77b941eb33d2b |
135 |
| - "exonic=i" => \$status_exonic, |
| 122 | + "exonic=i" => \$status_exonic, |
136 | 123 | "rdir=s" => \$run_dir,
|
137 | 124 | "ref=s" => \$h38_REF,
|
138 | 125 | "smg=s" => \$db_smg,
|
@@ -550,16 +537,13 @@ sub bsub_strelka{
|
550 | 537 |
|
551 | 538 | my $sh_file=$job_files_dir."/".$current_job_file;
|
552 | 539 |
|
553 |
| -<<<<<<< HEAD |
554 | 540 | if($q_name eq "research-hpc")
|
555 | 541 | {
|
556 | 542 | $bsub_com = "bsub -q research-hpc -n 1 -R \"select[mem>30000] rusage[mem=30000]\" -M 30000000 -a \'docker(scao/dailybox)\' -o $lsf_out -e $lsf_err bash $sh_file\n"; }
|
557 | 543 | else {
|
558 | 544 | $bsub_com = "bsub -q $q_name -n 1 -R \"select[mem>30000] rusage[mem=30000]\" -M 30000000 -o $lsf_out -e $lsf_err bash $sh_file\n";
|
559 | 545 | }
|
560 |
| -======= |
561 | 546 | #$bsub_com = "LSF_DOCKER_PRESERVE_ENVIRONMENT=false bsub -q $q_name -n 1 -R \"select[mem>30000] rusage[mem=30000]\" -M 30000000 -a \'docker(scao/dailybox)\' -o $lsf_out -e $lsf_err bash $sh_file\n";
|
562 |
| ->>>>>>> a4a79c27f83ba1ebce51a18020e77b941eb33d2b |
563 | 547 |
|
564 | 548 | $bsub_com = "bsub -g /$compute_username/$group_name -q $q_name -n 1 -R \"select[mem>30000] rusage[mem=30000]\" -M 30000000 -a \'docker(scao/dailybox)\' -o $lsf_out -e $lsf_err bash $sh_file\n";
|
565 | 549 | print $bsub_com;
|
@@ -695,19 +679,16 @@ sub bsub_varscan{
|
695 | 679 |
|
696 | 680 | my $sh_file=$job_files_dir."/".$current_job_file;
|
697 | 681 |
|
698 |
| -<<<<<<< HEAD |
699 | 682 | if($q_name eq "research-hpc")
|
700 | 683 | {
|
701 | 684 | $bsub_com = "bsub -q research-hpc -n 1 -R \"select[mem>30000] rusage[mem=30000]\" -M 30000000 -a \'docker(scao/dailybox)\' -w \"$hold_job_file\" -o $lsf_out -e $lsf_err bash $sh_file\n"; }
|
702 | 685 | else { $bsub_com = "bsub -q $q_name -n 1 -R \"select[mem>30000] rusage[mem=30000]\" -M 30000000 -w \"$hold_job_file\" -o $lsf_out -e $lsf_err bash $sh_file\n"; }
|
703 | 686 | print $bsub_com;
|
704 | 687 | system ($bsub_com);
|
705 |
| -======= |
706 | 688 | $bsub_com = "LSF_DOCKER_PRESERVE_ENVIRONMENT=false bsub -g /$compute_username/$group_name -q $q_name -n 1 -R \"select[mem>30000] rusage[mem=30000]\" -M 30000000 -a \'docker(scao/dailybox)\' -o $lsf_out -e $lsf_err bash $sh_file\n";
|
707 | 689 |
|
708 | 690 | print $bsub_com;
|
709 | 691 | system ($bsub_com);
|
710 |
| ->>>>>>> a4a79c27f83ba1ebce51a18020e77b941eb33d2b |
711 | 692 |
|
712 | 693 | #$bsub_com = "bsub < $job_files_dir/$current_job_file\n";
|
713 | 694 | #system ( $bsub_com );
|
@@ -1287,21 +1268,17 @@ sub bsub_mutect{
|
1287 | 1268 |
|
1288 | 1269 | $bsub_com = "LSF_DOCKER_PRESERVE_ENVIRONMENT=false bsub -g /$compute_username/$group_name -q $q_name -n 4 -R \"select[mem>30000] rusage[mem=30000]\" -M 30000000 -a \'docker(scao/dailybox)\' -o $lsf_out -e $lsf_err bash $sh_file\n";
|
1289 | 1270 |
|
1290 |
| -<<<<<<< HEAD |
1291 | 1271 | #$bsub_com = "bsub < $job_files_dir/$current_job_file\n";
|
1292 |
| - my $sh_file=$job_files_dir."/".$current_job_file; |
| 1272 | + my $sh_file=$job_files_dir."/".$current_job_file; |
1293 | 1273 | if($q_name eq "research-hpc")
|
1294 | 1274 | {
|
1295 | 1275 | $bsub_com = "bsub -q research-hpc -n 4 -R \"select[mem>30000] rusage[mem=30000]\" -M 30000000 -a \'docker(scao/dailybox)\' -w \"$hold_job_file\" -o $lsf_out -e $lsf_err bash $sh_file\n"; }
|
1296 | 1276 | else
|
1297 | 1277 | {
|
1298 | 1278 | $bsub_com = "bsub -q $q_name -n 4 -R \"select[mem>30000] rusage[mem=30000]\" -M 30000000 -w \"$hold_job_file\" -o $lsf_out -e $lsf_err bash $sh_file\n";
|
1299 | 1279 | }
|
1300 |
| - print $bsub_com; |
1301 |
| -======= |
1302 | 1280 | print $bsub_com;
|
1303 | 1281 | system ($bsub_com);
|
1304 |
| ->>>>>>> a4a79c27f83ba1ebce51a18020e77b941eb33d2b |
1305 | 1282 |
|
1306 | 1283 | }
|
1307 | 1284 |
|
|
0 commit comments