Skip to content

impute-gene-expression
impute-gene-expression
Commit 148ef2a8 authored by Kevin Kunzmann's avatar Kevin Kunzmann
Browse files
Options

add MAF

parent b7409884
Show whitespace changes
Inline Side-by-side
Showing with 1 addition and 1 deletion
Snakefile
View file @ 148ef2a8
...@@ -72,7 +72,7 @@ rule vcf_to_dosages: ...@@ -72,7 +72,7 @@ rule vcf_to_dosages:
echo "decompress and filter out: 'MAF[0]<{config[min_MAF]} | INFO<{config[min_INFO]} | TYPE!="snp" | N_ALT!=1' ..." echo "decompress and filter out: 'MAF[0]<{config[min_MAF]} | INFO<{config[min_INFO]} | TYPE!="snp" | N_ALT!=1' ..."
bcftools filter -e 'MAF[0]<{config[min_MAF]} | INFO<{config[min_INFO]} | TYPE!="snp" | N_ALT!=1' {input.vcf_gz_file} > $prefix/chr{wildcards.i}_.vcf bcftools filter -e 'MAF[0]<{config[min_MAF]} | INFO<{config[min_INFO]} | TYPE!="snp" | N_ALT!=1' {input.vcf_gz_file} > $prefix/chr{wildcards.i}_.vcf
echo "add computing MAFs ..." echo "add computing MAFs ..."
bcftools $prefix/chr{wildcards.i}.vcf_ > $prefix/chr{wildcards.i}.vcf bcftools +fill-tags $prefix/chr{wildcards.i}.vcf_ > $prefix/chr{wildcards.i}.vcf
rm $prefix/chr{wildcards.i}_.vcf rm $prefix/chr{wildcards.i}_.vcf
echo 'querying dosages ...' echo 'querying dosages ...'
bcftools query -f '%CHROM %ID %POS %REF %ALT %INFO/MAF [%DS ]\n' $prefix/chr{wildcards.i}.vcf > $prefix/chr{wildcards.i}.dosage.txt bcftools query -f '%CHROM %ID %POS %REF %ALT %INFO/MAF [%DS ]\n' $prefix/chr{wildcards.i}.vcf > $prefix/chr{wildcards.i}.dosage.txt
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment