Skip to content

Commit

Permalink
fix 1025
Browse files Browse the repository at this point in the history
  • Loading branch information
[email protected] committed Sep 5, 2024
1 parent bfc7cb1 commit cda0065
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion snakePipes/shared/rules/ChIP_peak_calling.snakefile
Original file line number Diff line number Diff line change
Expand Up @@ -122,7 +122,7 @@ rule MACS2_peak_qc:
# compute peak genome coverage
peak_len=`awk '{{total+=$3-$2}}END{{print total}}' {params.peaks}`
genome_size=`awk '{{total+=$3-$2}}END{{print total}}' {params.genome_index}`
genome_size=`awk '{{total+=$3-$2}}END{{printf( "%14d",total) }}' {params.genome_index}`
genomecov=`bc -l <<< "$peak_len/$genome_size"`
# write peak-based QC metrics to output file
Expand Down

0 comments on commit cda0065

Please sign in to comment.