couldn't run CNVnator

[menenuh]

Hello,
Does anyone have experience with CNVnator?
It requires the ROOT package, and after two unsuccessful attempts I managed to install it finally. But now I'm having problems running it.
Am I the only one trying/failing?

[southan]

Before you use cnvnator, you can go to (or not) "root/bin/" and type this line:
source thisroot.sh
It will automatically set your ROOTSYS.

Hello,
Does anyone have experience with CNVnator?
It requires the ROOT package, and after two unsuccessful attempts I managed to install it finally. But now I'm having problems running it.
Am I the only one trying/failing?

[milesgr]

How did you get root to install? After like 15 minutes, I got the following:


/usr/bin/ld: /usr/local/lib/libfftw3.a(map-r2r-kind.o): relocation R_X86_64_32S against `a local symbol' can not be used when making a shared object; recompile with -fPIC
/usr/local/lib/libfftw3.a: could not read symbols: Bad value
collect2: ld returned 1 exit status
gmake: *** [lib/libFFTW.so] Error 1


I downloaded the source...then did
export ROOTSYS=<path>/root

I then ran these commands from the root directory...
./configure
gmake

gmake resulted in the above error. Any thoughts?

[milesgr]

Never mind, I figured it out. If anyone needs help, please feel free to post here.

[CNVboy]

yeah, i have problem with the root package....how to change the path for the library? to edit the .bash_profile?

[southan]

I use this line:
source /path/root/bin/thisroot.sh

[CNVboy]

thx southan. I've been OK with both root and CNVnator compiling.
But now the question is, the command is like:
./cnvnator -root NA12878.root -chrom chr1 chr2 chr3 -tree NA12878_ali.bam

what's this NA12878.root? Is it some output of root package?

the documentation of CNVnator is really briefly-written and could be confusing to beginners like me.

[ElMichael]

thx southan. I've been OK with both root and CNVnator compiling.
But now the question is, the command is like:
./cnvnator -root NA12878.root -chrom chr1 chr2 chr3 -tree NA12878_ali.bam

what's this NA12878.root? Is it some output of root package?

It's just an output file that would be generated by the execution of the first CNVnator command (where extracted data from the input BAM file is written?) and, as far as I understand, would be used as a part of the input during the next CNVnator steps.

[CNVboy]

thx Michael.

Now I've gone through the whole steps of CNVnator and seems quite good. The only problem now is at the final step about visualizing the CNV. I enter a genomic region at the prompt but got NO output shown as below:
./cnvnator -root NA18507.root -chrom chr1 -view 100
>chr1:246974801-247024800
>

It has already called CNV, but I just cannot visualize nor genotype.

btw, I'm wondering where is the output stored? Is it within some files?

thx

[CNVboy]

Just now I tried the whole genome (no longer specific to chr1): ./cnvnator -root NA18507. root -tree ***.bam

Then it pop out:
Out of bound coordinate 199501844 for chromosome 'chr3'.
Out of bound coordinate 100338932 for chromosome 'chr15'.
Out of bound coordinate 100338932 for chromosome 'chr15'.
Out of bound coordinate 100338931 for chromosome 'chr15'.
Out of bound coordinate 100338932 for chromosome 'chr15'.
Out of bound coordinate 78774759 for chromosome 'chr17'.
Out of bound coordinate 199501844 for chromosome 'chr3'.

Then it goes well:
........
........
Filling and saving tree for 'chr5' ...
Done with the tree.
Filling and saving tree for 'chr6' ...
Done with the tree.
Filling and saving tree for 'chr7' ...
Done with the tree.
Filling and saving tree for 'chr8' ...
Done with the tree.
Filling and saving tree for 'chr9' ...
Done with the tree.
Filling and saving tree for 'chr10' ...
Done with the tree.
.........
.........
I don't know what does this 'out of bound coordinate" mean....Anyone has any idea? thx

[southan]

thx Michael.

Now I've gone through the whole steps of CNVnator and seems quite good. The only problem now is at the final step about visualizing the CNV. I enter a genomic region at the prompt but got NO output shown as below:
./cnvnator -root NA18507.root -chrom chr1 -view 100
>chr1:246974801-247024800
>

It has already called CNV, but I just cannot visualize nor genotype.

btw, I'm wondering where is the output stored? Is it within some files?

thx

cnvnator -root NA18507.root -view 100

That's Ok.
Are you using ssh? If so, you should use ssh -X.

[CNVboy]

yeah, I'm using ssh.

I google ssh option, and -X seems to mean to enable X11 forwarding....But I'm still confused about how this has anything to do with my problem of visualizing/genotyping/storation....

thx

[gprakhar]

Hello CNVboy,

If you want to view anything other than text while using ssh, you will have to use the "-X" option.

Regards,
Prakhar Gaur

[CNVboy]

thx guys.

I tried using ssh -X, but still doesn't work.

[xug@pem610-006 src]$ ./cnvnator -root NA18507.root -view 100
>chr12 11396601 11436500
>

I can see nothing except the new prompt.

[gprakhar]

Hello CNVboy,

Can you please give a list of all the commands you are using for running CNVnator from starting,
and some details about the input files?

Regards,
Prakhar Gaur