bash - Extracting SNPs and INDELs from a VCF using VCF tools -


how extract these snps listed in excel file this vcf file using vcf tools?

i have tried using script, took regions, not actual variants.

#! /bin/bash  module load java module load vcftools  pvcf='yeast_combined.g.vcf.filtered.vcf' cvcf='yeast_combined.updated2.g.vcf.filtered.vcf'  #extracts snp peak regions  vcftools --vcf $pvcf --bed snp_regions.bed --recode --recode-info-all --out psnps vcftools --vcf $cvcf --bed snp_regions.bed --recode --recode-info-all --out csnps  #extracts indel region  vcftools --vcf $pvcf --bed ind_regions.bed --recode --recode-info-all --out pinds vcftools --vcf $cvcf --bed ind_regions.bed --recode --recode-info-all --out cinds 


Comments

Popular posts from this blog

Retrieving ETA (estimated time of arrival) with Google Distance Matrix API and public transit as transport mode -

javascript - Confirm a form & display message if form is valid with JQuery -

ionic framework - Meteor - Error: Failed to execute 'insertBefore' on 'Node' -