Skip to content

Instantly share code, notes, and snippets.

@jodyphelan
Last active January 7, 2020 19:37
Show Gist options
  • Save jodyphelan/f90cddccd56576171f89c1527056665b to your computer and use it in GitHub Desktop.
Save jodyphelan/f90cddccd56576171f89c1527056665b to your computer and use it in GitHub Desktop.
bcftools view -V indels merged.raw.vcf.gz | setGT.py | bcftools view -a | bcftools filter -e 'GT="het"' -S . | bcftools view -i 'F_PASS(GT!="mis")>0.9' | bcftools view -c 1 | bcftools norm -f ../reference/Mycobacterium_abscessus.7396_7_44.dna.toplevel.fa | bcftools view -Oz -o merged.filt.vcf.gz
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment