filtering large vcf. GQ and DP
1
0
Entering edit mode
6 weeks ago
kazakkaty • 0

hi everyone, i'm new in bioinformatics, i have a big VCF file, i get dp from there using dp = extract.gt (vcf, element = 'DP', as.numeric = TRUE) dpp = apply (dp, 2, sum) and I build a plot barplot (sort (dpp), las = 3, main = "DP", col = 1:12), I need to take averages from this graph and do a GQ analysis on them, how can I do that?

DP GQ vcf • 149 views
ADD COMMENT
0
Entering edit mode
6 weeks ago

use bctools stats followed by plot-vcfstats https://samtools.github.io/bcftools/bcftools.html#plot-vcfstats

ADD COMMENT

Login before adding your answer.

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