Question: Remove columns from the VCF file using vcftools
0
gravatar for mayara_salvian
19 months ago by
mayara_salvian0 wrote:

I have a vcf file and I would like to remove four columns (CHROM,POS,ID,REF). I used the following command to extract these columns of my vcf file, but it's not working.

vcf-subset -c CHROM,POS,ID,REF my.vcf > out.vcf

What else should I add to that command?

vcftools • 1.6k views
ADD COMMENTlink modified 19 months ago by genomax76k • written 19 months ago by mayara_salvian0
1

if you want to remove headers and first 4 columns:

 $ awk -v OFS="\t" '!/##/ {$1=$2=$3=$4="";print}' test.vcf|sed 's/^\s\+//g'

If you want to retain headers, but remove first 4 columns:

$ awk -v OFS="\t" '!/##/ {$1=$2=$3=$4=""}1' test.vcf |sed 's/^\s\+//g'

with sed:

$ sed -e '/##/! s/^\([^\t]*\t\)\{4\}//g' test.vcf | grep -v "##"
ADD REPLYlink modified 19 months ago • written 19 months ago by cpad011212k
2
gravatar for Pierre Lindenbaum
19 months ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum125k wrote:
cut  --complement -f 1-4 my.vcf | grep -v "^##"> out.vcf
ADD COMMENTlink modified 19 months ago • written 19 months ago by Pierre Lindenbaum125k
1

shouldn't it be (as OP wants to remove the columns, not retain)?

cut --complement -f1-4  my.vcf | grep -v "^##"> out.vcf
ADD REPLYlink modified 19 months ago • written 19 months ago by cpad011212k

yes, you're right ! I didn't notice. Thanks.

ADD REPLYlink written 19 months ago by Pierre Lindenbaum125k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 924 users visited in the last hour