File tree Expand file tree Collapse file tree 2 files changed +14
-12
lines changed Expand file tree Collapse file tree 2 files changed +14
-12
lines changed Original file line number Diff line number Diff line change @@ -2,9 +2,10 @@ rule order_markers:
2
2
input :
3
3
datacall = "2_Filtering/data.filtered.lepmap3.gz" ,
4
4
filt_map = "LOD.master"
5
- output : "4_OrderMarkers/ordered.{lg_range}"
5
+ output :
6
+ lg = "4_OrderMarkers/ordered.{lg_range}" ,
7
+ runlog = temp ("4_OrderMarkers/logs/ordered.{lg_range}.running" )
6
8
log :
7
- runlog = temp ("4_OrderMarkers/logs/ordered.{lg_range}.running" ),
8
9
run = "4_OrderMarkers/logs/ordered.{lg_range}.log" ,
9
10
recomb = "4_OrderMarkers/recombination/ordered.{lg_range}.recombinations"
10
11
message : "Ordering linkage group {params.chrom} with {params.iterations} iterations"
@@ -15,10 +16,10 @@ rule order_markers:
15
16
threads : 2
16
17
shell :
17
18
"""
18
- zcat {input.datacall} | java -cp software/LepMap3 OrderMarkers2 map={input.filt_map} {params.extra} data=- numThreads={threads} numMergeIterations={params.iterations} chromosome={params.chrom} &> {log .runlog}
19
- sed -n '/\*\*\* LG \=/,$p' {log .runlog} > {output}
20
- grep "recombin" {log .runlog} > {log.recomb}
21
- awk '/#java/{{flag=1}} flag; /logL/{{flag=0}}' {log .runlog} > {log.run}
19
+ zcat {input.datacall} | java -cp software/LepMap3 OrderMarkers2 map={input.filt_map} {params.extra} data=- numThreads={threads} numMergeIterations={params.iterations} chromosome={params.chrom} &> {output .runlog}
20
+ sed -n '/\*\*\* LG \=/,$p' {output .runlog} > {output.lg }
21
+ grep "recombin" {output .runlog} > {log.recomb}
22
+ awk '/#java/{{flag=1}} flag; /logL/{{flag=0}}' {output .runlog} > {log.run}
22
23
"""
23
24
24
25
rule recomb_summary :
Original file line number Diff line number Diff line change @@ -3,9 +3,10 @@ rule reorder_markers:
3
3
datacall = "2_Filtering/data.filtered.lepmap3.gz" ,
4
4
filt_map = "LOD.master" ,
5
5
lg_order = "5_Trim/ordered.{lg_range}.trimmed"
6
- output : "6_OrderMarkers/ordered.{lg_range}"
6
+ output :
7
+ lg = "6_OrderMarkers/ordered.{lg_range}" ,
8
+ runlog = temp ("6_OrderMarkers/logs/ordered.{lg_range}.running" )
7
9
log :
8
- runlog = temp ("6_OrderMarkers/logs/ordered.{lg_range}.running" ),
9
10
run = "6_OrderMarkers/logs/ordered.{lg_range}.log" ,
10
11
recomb = "6_OrderMarkers/recombination/ordered.{lg_range}.recombination"
11
12
message : "Reordering linkage group {params.lg} with {params.iterations} iterations"
@@ -16,10 +17,10 @@ rule reorder_markers:
16
17
threads : 2
17
18
shell :
18
19
"""
19
- zcat {input.datacall} | java -cp software/LepMap3 OrderMarkers2 {params.extra} map={input.filt_map} data=- numThreads={threads} evaluateOrder={input.lg_order} numMergeIterations={params.iterations} &> {log .runlog}
20
- sed -n '/\*\*\* LG \=/,$p' {log .runlog} > {output}
21
- grep "recombin" {log .runlog} > {log.recomb}
22
- awk '/#java/{{flag=1}} flag; /logL/{{flag=0}}' {log .runlog} > {log.run}
20
+ zcat {input.datacall} | java -cp software/LepMap3 OrderMarkers2 {params.extra} map={input.filt_map} data=- numThreads={threads} evaluateOrder={input.lg_order} numMergeIterations={params.iterations} &> {output .runlog}
21
+ sed -n '/\*\*\* LG \=/,$p' {output .runlog} > {output.lg }
22
+ grep "recombin" {output .runlog} > {log.recomb}
23
+ awk '/#java/{{flag=1}} flag; /logL/{{flag=0}}' {output .runlog} > {log.run}
23
24
"""
24
25
25
26
rule reorder_summary :
You can’t perform that action at this time.
0 commit comments