Skip to content

Commit c16f68e

Browse files
committed
outputs from GLIMPSE2_phase and GLIMPSE2_ligate are automatically indexed
1 parent e369d02 commit c16f68e

File tree

1 file changed

+0
-3
lines changed

1 file changed

+0
-3
lines changed

gwaspy/imputation/glimpse2_impute.py

-3
Original file line numberDiff line numberDiff line change
@@ -236,8 +236,6 @@ def glimpse2_impute_phase(
236236
--threads 4 \
237237
--output {j.imputed_phased['bcf']}""")
238238

239-
j.command(f"""bcftools index -f {j.imputed_phased['bcf']} --output {j.imputed_phased['bcf.csi']} --threads {ncpu}""")
240-
241239
return j
242240

243241
# 7. Ligate imputed chunks
@@ -272,7 +270,6 @@ def ligate_chunks(
272270
j.command(f'echo "{imp_chunk_bcf_names}" > chunk_list.txt')
273271

274272
j.command(f"""GLIMPSE2_ligate --input chunk_list.txt --output {j.ligated['bcf']} --threads {ncpu}""")
275-
j.command(f"""bcftools index {j.ligated['bcf']} --output {j.ligated['bcf.csi']} --threads {ncpu}""")
276273

277274
b.write_output(j.ligated,
278275
f'{out_dir}/glimpse/{output_vcf_name}_{chrom}.phased.imputed')

0 commit comments

Comments
 (0)