Skip to content
Snippets Groups Projects
Commit e9e7281c authored by Sascha Meiers's avatar Sascha Meiers
Browse files

dirty hack to remove chrX and chrY

parent 2398f57b
No related branches found
No related tags found
No related merge requests found
......@@ -192,7 +192,10 @@ rule determine_initial_strand_states:
sce_command = "Rscript " + config["sce_script"]
shell:
"""
{params.sce_command} {input} {output}
echo "[Note] This is a dirty hack to remove chrX and chrY."
tmp=$(mktemp)
{params.sce_command} {input} $tmp
grep -vP '^Y|^chrY' $tmp | grep -vP '^X|^chrX' > {output}
"""
# Strandphaser needs a different input format which contains the path names to
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment