VCF filtering by depth
3
1
Entering edit mode
5.0 years ago

Hi, I am trying to filtering a vcf file by --minDP however the log file indicates that it is not filtering out any SNP, any clue about this problem? I tried different values for DP but is not working any of them (I do have DP information in my vcf file)

/data/storage/software/vcftools/bin/vcftools --vcf /home/jperez/croc/res_cornell/vcf_croc_addinCONTIG.vcf --recode-INFO-all --minDP 100 —-out minDP_100.vcf

Thanks

filtering by depth • 4.8k views
ADD COMMENT
0
Entering edit mode

you should try --min-meanDP instead of --minDP if you didn't solve that problem. I'm also wondering why vcftools write some argument making user a little bit confusing.

Like --min-meanDP & --minDP or --minQ and --minGQ

ADD REPLY
2
Entering edit mode
5.0 years ago
agata88 ▴ 810
I have a script that can do what you need. 

#! python

import sys
import os
import vcf

vcf_reader = vcf.Reader(open("file.vcf", 'r'))
vcf_writer = vcf.Writer(open('out.vcf', 'w'), vcf_reader)
for record in vcf_reader:
    if record.INFO['DP']>100:
        vcf_writer.write_record(record)

Hope it helps,
Best,
Agata
ADD COMMENT
0
Entering edit mode

How to filter AF in vcf file with that way? There is a KeyError.

ADD REPLY
0
Entering edit mode
5.0 years ago
igor 12k

VCFtools allows filtering. See the "--minDP" argument.

https://vcftools.github.io/man_latest.html

ADD COMMENT
0
Entering edit mode
5.0 years ago

Thanks Agata, I have mac and try to run the script and I got this message:

Traceback (most recent call last): File "DP.py", line 3, in <module> import vcf ImportError: No module named vcf

ADD COMMENT
0
Entering edit mode

You need to install PyVCF module for example from here https://pypi.python.org/pypi/PyVCF Make sure that it is compatible with your python version! Best, Agata

ADD REPLY

Login before adding your answer.

Traffic: 1069 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