Filtering the vcf file
0
0
Entering edit mode
5 weeks ago
trezini • 0

I'm wondering how snpgap and indelgap work. Let's say snpgap 1 indelgap 10, and the second question, do I understand correctly that max-idepth takes the average coverage of a variant in a sample, and if it is higher than the threshold that I set, it will not call it? If I have one sample and a threshold of, for example, 7000, then all variants with coverage greater than 7000 will not be called? And if I have two samples, in one of which the coverage of a variant is 6000 and in the other of the same variant 10000, then it will not be called in both cases?

vcf • 400 views
ADD COMMENT
0
Entering edit mode

I hope this may help:
what is the properties of filtering the vcf files

For the depth part, let us assume we used the following command bcftools view --max-depth 500 input.vcf, any variant where the depth for an individual exceeds 500 will be filtered out. This can be useful to avoid sites with unusually high coverage, which may indicate problematic regions, such as repetitive sequences or alignment artifacts.

ADD REPLY
0
Entering edit mode

I understand how max-depth works, but I asked about max-idepth, indelgap and snpgap

ADD REPLY

Login before adding your answer.

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