SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Correct options for CNVnator (http://seqanswers.com/forums/showthread.php?t=89829)

dfermin315 06-19-2019 06:25 AM

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.


All times are GMT -8. The time now is 05:18 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.