SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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

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

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

Greetings,

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:

Code:
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 \
	/nfs/wgs30x/bam.files/samp3.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:

Code:
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 cnvnator2VCF.pl 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
Reply

Tags
cnvnator

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 05:50 PM.


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