SeqSQC can be used for sample QC of Whole Genome Sequencing Data

I am using the R librabry SeqSQC from Bioconductor to do sample QC of around 50 Whole Genome Sequencing samples.

The vignette says Through incorporation a benchmark data assembled from the 1000 Genomes Project, it can accommodate NGS study of small sample size and low number of variants.

https://www.bioconductor.org/packages/devel/bioc/vignettes/SeqSQC/inst/doc/vignette.html

I could generate plots for MissingRate, SexCheck, Inbreeding, IBD, PCA, . Not sure if WGS samples with high coverage are giving the right results.

seqsqc wgs

• 583 views • link

written 2.8 years ago by ttom

0 answers

No answers yet.

More posts like this

Eigencor plot after limma::removeBatchEffect
written 1.5 years ago by f2ff7a46 •

I'm trying to create an eigencor plot of clinical variables against principle components using the PCAtools guidelines. [https://www.bioconductor.org/packages/devel/bioc/vignettes/PCAtools/inst/doc/PCAtools.html#correlate-the-principal-components-back-to-the-clinical-data][1] I want to first eliminate batch effects …
How to access TCGA samples that were treated with a specific drug?
written 2.3 years ago by Qroid •

I'm new to working with TCGA, but I'd like to look at RNAseq expression of a certain gene between responders and non-responders of a given …
Pre-filtering in DESeq2 vignette
written 2.9 years ago by dylannicoembros •

I am doing a pre-filtering of my counts following deseq2 vignette: > https://bioconductor.org/packages/devel/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#pre-filtering In the section Pre-filtering they select the significant counts (>= 10), but …
Tutorials for forging BSGenome packages?
written 7.3 years ago by am3

Lots of Bioconductor packages require a BSGenome data package to work, but they're only available for a limited number of organisms, so as far as …
Is it possible to get .bed file of sequenced regions after NGS experiment?
written 8.3 years ago by Adam •

Hello, As it is in post title, Is it possible to get .bed file of sequenced regions after NGS experiment? I mean that you have …
Is there a name for these charts?
written 8.4 years ago by russhh

There's some plots on the fgsea [vignette][1] that I'd like to refer to by name. It's the chart returned by the plotGseaTable function that I'm …
What is the meaning of Bumphunter output?
written 8.5 years ago by boyu93 •

I tried using ChAMP to analyze differential Methylation Regions (DMR).[ChAMP][1] The method I used was based on Bumphunter. But its output made me confused. Especially …
How do I find out regions of methylation enrichment using MEDIPS bioconductor package?
written 9.8 years ago by lakhujanivijay

I have a single MEDIP-seq sample and I am only interested in finding out regions of methylation. How do I go about it? Since, I …
Replicates in ExpressionSet
written 10.0 years ago by jacobsen.jeremy •

I'm trying to benchmark TDAracne for gene network reconstruction on a simulated dataset. TDAracne takes in data as an ExpressionSet. This is my first time …
Tool: Drug repositioning based on co-expressed gene set enrichment analysis, cogena
written 11.3 years ago by Zhilong Jia

[Cogena][1] is a framework for co-expressed gene set enrichment analysis. It can be used for pathway/GO analysis or drug repositioning, etc. as smaller scale, but …