Advertisement
Guest User

Untitled

a guest
Aug 18th, 2019
111
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.31 KB | None | 0 0
  1. # VCF file = test.filt.vcf.gz
  2. # replace where necessary
  3.  
  4. (printf "chr\npos\nref\n"; bcftools query -l test.filt.vcf.gz) | datamash transpose > test.filt.mat.bin
  5. bcftools query -f '%CHROM\t%POS\t%REF[\t%GT]\n' test.filt.vcf.gz | tr '|' '/' | sed 's/\.\/\./N/g' | sed 's/0\/0/0/g' | sed 's/.\/./1/g' >> test.filt.mat.bin
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement