Skip to content

Commit

Permalink
Update integrated_annotation.wdl
Browse files Browse the repository at this point in the history
  • Loading branch information
stellaning1120 authored Aug 18, 2023
1 parent 420b42c commit 3ffa949
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions AnnotateBed/integrated_annotation.wdl
Original file line number Diff line number Diff line change
Expand Up @@ -153,10 +153,10 @@ task samtools_coverage {
Int disk_size
}
command <<<
sed 's/^chr//' ~{bam_file} > rm_chr_tmp.bed
for region in `awk '{print $1":"$2"-"$3}' rm_chr_tmp.bed`
sed 's/^chr//' ~{bed_to_annotate} > rm_chr_tmp.bed
for region in `awk '{print $1":"$2"-"$3}' ${rm_chr_tmp.bed}`
do
samtools coverage rm_chr_tmp.bed -r ${region} >> tmp
samtools coverage ~{bam_file} -r ${region} >> tmp
done
awk '/startpos/&&c++>0 {next} 1' tmp > ~{sample_id}.coverage.txt
>>>
Expand Down

0 comments on commit 3ffa949

Please sign in to comment.