Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
CNVnator need help jchoo Bioinformatics 4 09-23-2016 04:18 AM
CNVnator concepts eleven Bioinformatics 0 04-01-2013 11:25 PM
A problem in CNVnator yuhao Bioinformatics 1 01-30-2013 12:07 AM
Help! a problem in CNVnator yuhao Bioinformatics 0 07-12-2012 07:17 AM
Tophat options to report unaligned reads and controlling Bowtie options Siva Bioinformatics 0 10-15-2010 08:38 PM

Thread Tools
Old 06-19-2019, 07:25 AM   #1
Junior Member
Location: Michigan, US

Join Date: Mar 2017
Posts: 4
Default Correct options for CNVnator


I'm trying to run CNVnator on a cohort of ~300 BAM files.

For testing purposes I'm trying to run it on just 3 BAM files first.
Here are the commands I'm running:

export HS37d5=/data/fastas/hs37d5

cnvnator -root chr22.root -genome GRCh37 -chrom 22 -tree \
	/nfs/wgs30x/bam.files/samp1.bam \
	/nfs/wgs30x/bam.files/samp2.bam \

cnvnator -root chr22.root -genome GRCh37 -chrom 22 -his 1000 -d $HS37D5

cnvnator -root chr22.root -genome GRCh37 -chrom 22 -stat 1000

cnvnator -root chr22.root -genome GRCh37 -chrom 22 -partition 1000

cnvnator -root chr22.root -genome GRCh37 -chrom 22 -call 1000 > chr22.cnvnator.out

Now all of this runs and I don't get an error (as far as I can tell).
The last command I try is:

cnvnator -pe /nfs/wgs30x/bam.files/samp1.bam \
   /nfs/wgs30x/bam.files/samp2.bam \
   /nfs/wgs30x/bam.files/samp3.bam \
   -qual 20 -over 0.8 -root chr22.root -f chr22.support_reads;
This command never finishes. I left it running for 36 hours and it didn't produce anything.

I believe I'm screwing up the cnvnator -pe command but I can't figure out what the correct command should be.

Using the produces a VCF file but it doesn't contain any read count information.

Does anyone have any suggestions as to what I'm doing wrong?
In case it helps I'm using CNVnator v0.4 (latest download from github).
The system is a RHEL 7.2 box fully patched.

Thank you in advance for any and all help.
dfermin315 is offline   Reply With Quote


Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 08:33 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO