diff --git a/src/download_encode_eclip/merge_eclip.py b/src/download_encode_eclip/merge_eclip.py index 0e14931127702065cff3164712aad169586457d8..db0df8d8a231b2f7d2fa4cf6b33b6bd8044c3d6f 100644 --- a/src/download_encode_eclip/merge_eclip.py +++ b/src/download_encode_eclip/merge_eclip.py @@ -143,8 +143,8 @@ def merge_files(dfiles: Dict[str, List[Path]], target: str) -> None: list_file = " ".join(map(str, dfiles[k])) cmd = f"zcat {list_file} | sort -k1,1 -k2,2n > {tmp_out}" cmd2 = f"bedtools merge -i {tmp_out} -s -c 4,5,6 " \ - f"-o collapse,mean,distinct | cut -f1-3,5-7 | gzip -c " \ - f"> {res_out}" + f"-o collapse,mean,distinct | cut -f1-3,5-7 | sed 's/^chr//g' " \ + f"| gzip -c > {res_out}" sp.check_output(cmd, shell=True) sp.check_output(cmd2, shell=True) tmp_out.unlink()