Header Leaderboard Ad

Collapse

couldn't run CNVnator

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.

                            [[email protected] 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
                                A Brief Overview and Common Challenges in Single-cell Sequencing Analysis
                                by seqadmin


                                ​​​​​​The introduction of single-cell sequencing has advanced the ability to study cell-to-cell heterogeneity. Its use has improved our understanding of somatic mutations1, cell lineages2, cellular diversity and regulation3, and development in multicellular organisms4. Single-cell sequencing encompasses hundreds of techniques with different approaches to studying the genomes, transcriptomes, epigenomes, and other omics of individual cells. The analysis of single-cell sequencing data i...

                                01-24-2023, 01:19 PM
                              • seqadmin
                                Introduction to Single-Cell Sequencing
                                by seqadmin
                                Single-cell sequencing is a technique used to investigate the genome, transcriptome, epigenome, and other omics of individual cells using high-throughput sequencing. This technology has provided many scientific breakthroughs and continues to be applied across many fields, including microbiology, oncology, immunology, neurobiology, precision medicine, and stem cell research.

                                The advancement of single-cell sequencing began in 2009 when Tang et al. investigated the single-cell transcriptomes
                                ...
                                01-09-2023, 03:10 PM
                              • seqadmin
                                AVITI from Element Biosciences: Latest Sequencing Technologies—Part 6
                                by seqadmin
                                Element Biosciences made its sequencing market debut this year when it released AVITI, its first sequencer. The AVITI System uses avidity sequencing, a novel sequencing chemistry that delivers higher quality data, decreases cycle times, and requires lower reagent concentrations. This new instrument reportedly features lower operating and start-up costs while maintaining quality sequencing.

                                Read type and length
                                AVITI is a short-read benchtop sequencer that also offers an innovative...
                                12-29-2022, 10:44 AM

                              ad_right_rmr

                              Collapse
                              Working...
                              X