SNP data filtering from VCF file
0
0
Entering edit mode
7.3 years ago
sukesh1411 ▴ 30

Hi

I am sorry although there are so many questions on this query i could not find the relavant answer. I used vcftools to extract only snps from vcf file using below command

command: vcftools --vcf variantsfilterQDPMQ.vcf --remove-indels --recode --recode-INFO-all --out onlysnps

I got the vcf file which has the information about each nucleotide with snp as well as without snp. I want vcf file only with snp information alone.

I wrote a perl program, its taking lot of time since i got a big vcf file which is 10gb..

Thank you all

SNP • 3.0k views
ADD COMMENT
0
Entering edit mode

You have already asked some questions here:

please update/comment or validate (green mark on the left) the answers.

ADD REPLY
0
Entering edit mode

Hi, what variant caller do you use? Some of them have an option to call INDELS and SNP separately - for example VarScan..

ADD REPLY
0
Entering edit mode

Hi,

I took vcf file from rice 3000 genome project database. Can i use varscan for exporting only snps from vcf file?

Thanks

ADD REPLY
1
Entering edit mode

Hi, no Varscan is variant caller - it generates VCF from BAM file (mpileup). Did you tried this Pierre solution?

ADD REPLY
0
Entering edit mode

Hi

Thank you Paul... I finished my job using awk script.. thanks to Pierre :)

ADD REPLY

Login before adding your answer.

Traffic: 1175 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6