Skip to content

Commit

Permalink
minor syntax changes
Browse files Browse the repository at this point in the history
  • Loading branch information
LindoNkambule committed Feb 22, 2024
1 parent 0d0a7ea commit df1a30a
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions gwaspy/imputation/glimpse2_impute.py
Original file line number Diff line number Diff line change
Expand Up @@ -87,7 +87,7 @@ def extract_ref_sites(
j.command(f"""bcftools view -G -Ou -o {j.ref_sites['bcf']} {ref_bcf['bcf']} """)
j.command(f"""bcftools index {j.ref_sites['bcf']} --output {j.ref_sites['bcf.csi']} --threads {ncpu}""")

j.command(f"""bcftools query -f'%CHROM\t%POS\t%REF,%ALT\n' {j.ref_sites['bcf']} | bgzip -c > {j.ref_sites['tsv.gz']}""")
j.command(f'''bcftools query -f'%CHROM\t%POS\t%REF,%ALT\n' {j.ref_sites['bcf']} | bgzip -c > {j.ref_sites['tsv.gz']}''')
j.command(f"""tabix -s1 -b2 -e2 {j.ref_sites['tsv.gz']}""")

return j
Expand Down Expand Up @@ -121,7 +121,7 @@ def compute_gls(
)

# compute GLs
j.command(f"""bcftools mpileup -f {ref_genome_fasta['fasta']} -I -E -a 'FORMAT/DP' -T {ref_sites['bcf']} -r {chrom} {sample_bam['bam']} -Ou | bcftools call -Aim -C alleles -T {ref_sites['tsv.gz']} -Ou -o {j.gl_vcf['bcf']}""")
j.command(f'''bcftools mpileup -f {ref_genome_fasta['fasta']} -I -E -a 'FORMAT/DP' -T {ref_sites['bcf']} -r {chrom} {sample_bam['bam']} -Ou | bcftools call -Aim -C alleles -T {ref_sites['tsv.gz']} -Ou -o {j.gl_vcf['bcf']}''')

j.command(f"""bcftools index {j.gl_vcf['bcf']} --output {j.gl_vcf['bcf.csi']} --threads {ncpu}""")

Expand Down

0 comments on commit df1a30a

Please sign in to comment.