Hello,
I want to detect CNVs from BAM files (from a batch and one BAM file independently), I read a couple of posts about it, I read about recent tools: biorxiv review and Nature review. Read-depth is the most common approach and I wondered why not normalise by read-depth then use SAMTools to plot the coverage. Comparing with a control, you can infer the CNVs from a BAM file. Sound logical however I feel like I'm trying to reinvent the wheel since there are tools built for this purpose. Do you have any thoughts on it?
I know that it is possible to infer the CNVs from multiple samples but is it correct to infer CNV's from one BAM file by comparing the read-depth between a patient and control (after normalization)?
•
written
by
Sammy •
1 answer
Maybe cnv-kit could be the better choice.
I usually do something even more simple for 2 samples - GC-content normalization and manual selection of regions with copy-numbers of approx 3 or 1. But yeap it is a powerful magic.
In general, CNV callers which are based on read depth and do not use a large cohort for normalization are suitable only for the detection of large CNVs (>10 exonic regions involved). Smaller CNVs become unreliable.
Log in to answer this question.
-
Herald: The Biostar Herald for Tuesday, December 19, 2023
written by Biostar
-
Herald: The Biostar Herald for Tuesday, September 26, 2023
written by Biostar
-
NGS - CNVnator
written by Nameless •
-
Help with CNV calling using ExomeDepth
written by luffy •
-
Tool: ClinCNV: CNV detection from short reads
written by German.M.Demidov
-
RNASeq read coverage in protein space of MSA?
written by CephBirk •
-
Best way of Merging Multiple Lanes generated by HiSeq 2000, index BAM or VCF?
written by venuraherath •
-
CONTRA for CNV detection. troubleshooting
written by Kizuna
-
Standalone Gc Correction Tools
written by Plantae
-
The Need Of Repeatmask For Read-Depth Based Cnv Detection Algorithm
written by Bioscientist
Are you trying to detect Germline CNVs? Just use cnvkit or ExomeDepth for this task!
Yes, germline, I have 2 BAMs from related individuals. Is it possible to calculate the LOH as well only from these two files?
you can call LOH even in one individual - try calling Runs-of-Homozygosity instead (plink can do this, also RohHunter https://github.com/imgag/ngs-bits/blob/master/doc/tools/RohHunter/index.md )
is it WGS, panel or WES?
you have only 2 samples, right?
WES, yes, I do have a patient and a healthy individual.