Hi all, i have analyzed 35 normal WGS sample for cnvs using cnvNator. Now i want to know common cnv region among these files so that those can be used as control panel.
Is the a tool or method to obtain these common region among all files at once?
thank you
•
written
by
popayekid55
GAIA will find recurrent copy number regions from your input data and assign a p-value to each region to help with filtering them. I believe you have the required information to run GAIA. The starting data is a row-binded list of all regions, with an extra column that indicates the sample from which the region derived. You decide your own cut-off points for gain (1) and loss (0) based on the segment mean.
A practical example for cancer is given here: C: How to extract the list of genes from TCGA CNV data
Kevin
i did not understand completely. I am looking for common (overlapping) cnv coordinates among these 35 files.
GAIA will find the common regions and assign a p-value based on how recurrent (frequent) they are in your dataset. The idea is that the more recurrent ones are more important.
If you literally just want to see the overlapping BED regions, even if it occurs in just 2 samples, then use the BEDTools solutions that were suggested. However, what would you do in the situation were one region is gain (amplified) in one sample but loss (deleted) in another? - does it make sense to merge these in light of what is your downstream analysis plan?
Log in to answer this question.
We are not necessarily familiar with the output format of cnvNator, so it would be best if you could elaborate on which files you have.
output will be converted into bed file format like below
chromosome start and end of cnv, type of cnv and a score
bedtools multiinter will help you
Parse the chr and start end form CNVnator results then overlap the output file using
Bedtools multiIntersectBedorBedops.