Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • couldn't run CNVnator

    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?

  • #2
    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.


    Originally posted by menenuh View Post
    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?

    Comment


    • #3
      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?

      Comment


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

        Comment


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

          Comment


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

            Comment


            • #7
              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.

              Comment


              • #8
                Originally posted by CNVboy View Post
                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.

                Comment


                • #9
                  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

                  Comment


                  • #10
                    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

                    Comment


                    • #11
                      Originally posted by CNVboy View Post
                      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.

                      Comment


                      • #12
                        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

                        Comment


                        • #13
                          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

                          Comment


                          • #14
                            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.

                            Comment


                            • #15
                              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

                              Comment

                              Latest Articles

                              Collapse

                              • seqadmin
                                Strategies for Sequencing Challenging Samples
                                by seqadmin


                                Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                                03-22-2024, 06:39 AM
                              • seqadmin
                                Techniques and Challenges in Conservation Genomics
                                by seqadmin



                                The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                                Avian Conservation
                                Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                                03-08-2024, 10:41 AM

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by seqadmin, Yesterday, 06:37 PM
                              0 responses
                              10 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, Yesterday, 06:07 PM
                              0 responses
                              9 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-22-2024, 10:03 AM
                              0 responses
                              49 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-21-2024, 07:32 AM
                              0 responses
                              67 views
                              0 likes
                              Last Post seqadmin  
                              Working...
                              X