We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 789c12e commit 180b6a7Copy full SHA for 180b6a7
1 file changed
main.nf
@@ -2055,7 +2055,7 @@ process bedtools {
2055
script:
2056
"""
2057
## Create genome file from bam header
2058
- samtools view -H ${bam} | grep @SQ | sed 's/@SQ\tSN:\|LN://g' > genome.txt
+ samtools view -H ${bam} | grep '@SQ' | sed 's#@SQ\tSN:\\|LN:##g' > genome.txt
2059
2060
## Run bedtools
2061
bedtools coverage -nonamecheck -g genome.txt -sorted -a ${anno_file} -b ${bam} | pigz -p ${task.cpus - 1} > "${bam.baseName}".breadth.gz
0 commit comments