Contra control file
2
0
Entering edit mode
2.3 years ago
BAGeno ▴ 180

Hi,

I am want to find copy number variations in my exome sequencing data. For this purpose i would like to use contra. I do not have control bam so can I use bam file of 1000 genome as my control?

contra cnv • 546 views
ADD COMMENT
1
Entering edit mode

What makes you think that a completely unrelated sequencing experiment (different experiment exome vs WGS, different kit, different sequencing regime etc.) would serve as an adequate control?

ADD REPLY
0
Entering edit mode
2.3 years ago
kanika.151 ▴ 120

To find CNVs you need control samples from the same sequencing type. For example, if you have a case file sequenced by Sure Select V4 then you need control from Sure Select V4. You should rather use tools which would detect CNV without control bams. I think cn.mops detects CNVs without control bam files. But, please check before. Also, There might be control samples available on NCBI which would have undergone identical sequencing experiment to yours. Try find those samples for your controls. (Do this only if you could not find any other tool for detection of CNV).

ADD COMMENT
0
Entering edit mode
2.2 years ago
halina • 0

Hi BAGeno, Contra has a baseline.py script that can create a control from your sample bams. I'm running into some errors with it at the moment but the documentation is here: http://contra-cnv.sourceforge.net/ at the end of the page.

ADD COMMENT

Login before adding your answer.

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