This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extracting Novel SNPs from VEP output

I run command based VEP for annotation. I want to extract only NOVEL SNPSs according to rs ID in first column. How I can do this.

1_17538_C/T     1:17538 T       ENSG00000278267 ENST00000619216 Transcript      upstream_gene_variant   -
1_17538_C/T     1:17538 T       ENSG00000223972 ENST00000456328 Transcript      downstream_gene_variant -
1_17538_C/T     1:17538 T       ENSG00000227232 ENST00000488147 Transcript
1_17538_C/T     1:17538 T       ENSG00000223972 ENST00000450305 Transcript      downstream_gene_variant -
rs201057270     1:17614 A       ENSG00000278267 ENST00000619216 Transcript      upstream_gene_variant   -       - rs201057270     1:17614 A       ENSG00000223972 ENST00000456328 Transcript      downstream_gene_variant -     rs201057270     1:17614 A       ENSG00000227232 ENST00000488147 Transcript
rs201057270     1:17614 A       ENSG00000223972 ENST00000450305 Transcript      downstream_gene_variant -
1_19342_G/A     1:19342 A       ENSG00000278267 ENST00000619216 Transcript      upstream_gene_variant   -
1_19342_G/A     1:19342 A       ENSG00000223972 ENST00000456328 Transcript      downstream_gene_variant -
1_19342_G/A     1:19342 A       ENSG00000227232 ENST00000488147 Transcript
snp
• 2,533 views • link
updated by Ram Ram • written by Ishtiaq Ahmad Khan Ishtiaq Ahmad Khan

1 answer 

grep -v '^rs'
• 0 views • link
updated by Ram Ram • written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this