filtering large vcf. GQ and DP
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?
• 1,025 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
kazakkaty •
•
written
by
kazakkaty •
1 answer
use bctools stats followed by plot-vcfstats https://samtools.github.io/bcftools/bcftools.html#plot-vcfstats
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
Putting colors of different samples on PCA plot
written by prity6459 •
-
filtering lowly expressed reads
written by Shaimaa Gamal •
-
filter barplot large VCF
written by kazakkaty •
-
RNASeq - How to compare the expression of a single gene across batches
written by oakhamwolf •
-
Error in extract.gt(vcf) - ID column contains non-unique names
written by gimeno •
-
How to label custom multiple gene lists object into volcano plot?
written by choijamtsm
-
Should I remove samples after normalization of miRNA seq read counts ?
written by Björn •
-
Ma plot labelling
written by op263 •
-
Signature genes from RNA-seq time series (Developmental)
written by Dataminer
-
Function for barplot
written by dktarathym •
