Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • CREST structural variation software error

    Hi all, I'm getting an error using CREST's example pipeline for structural variation detection

    Code:
    $CREST.pl -f tumor.bam.cover -d tumor.bam -g germline.bam --ref_genome \
            hg18.fa -t /genome/hg18.2bit -r chr1
    Use of uninitialized value $tid in array element at ./CREST.pl line 338.
    It's not obvious to me from the code how $tid failed to initialize, and my search for similar errors came up empty. Before I try to isolate the problem I thought I would raise the error here to see if anyone else has seen it before.

    Thanks!

  • #2
    Hm, perhaps the authors don't hang out here. I'll send a message directly. If I get a resolution I'll post here.

    Meanwhile, any CREST success stories?

    Comment


    • #3
      I have run it and got a few of those errors. It still completed though and I've compared the results with what I found with BreakDancer - didn't find anything new, BreakDancer picked up around 80% of CREST results and the others had been picked up by Freec. Best thing was the list was much shorter than BreakDancer, of the 30ish predictions only 2 look a bit dubious in the mappings

      Comment


      • #4
        Oh! the error I saw only occurred when I attempted to use the -r option. Without an -r parameter everything worked fine, cool.

        But! It only works for paired-end input in my hands. You're supposed to be able to use CREST in --nopaired mode, but when I attempt --nopaired mode extractSClip.pl gives me a torrent of errors like:
        Code:
        Use of uninitialized value $sdna in substr at /usr/local/lib64/perl5/Bio/DB/Bam/AlignWrapper.pm line 243.
        extractSClip does provide some output, in spite of the errors. But, when extractSClip's output is supplied to CREST.pl, CREST.pl returns a null file for xxx.predSV.txt and another mighty stream to STDERR.

        So, it seems my issue remains. Has anyone successfully used CREST for unpaired SV analysis?
        Last edited by bubfranks; 08-24-2011, 12:38 PM. Reason: punctuation

        Comment


        • #5
          hm, so it turned out that I misinterpreted the instructions. CREST instructions clearly state you need to use the same genome file in the CREST pipeline as you used in the initial mapping, but I took that to mean it was OK as long as the assemblies were the same. After using the exact same hg18.fasta file in BWA and CREST then the program started to work for me.

          Comment


          • #6
            I'm trying as well but getting a different error than you, regarding the BLAT server.

            /Crest$ CREST.pl -f tumor.bam.cover -d tumor.bam -g germline.bam --ref_genome /home/aquinom/hg18.fa -t /home/aquinom/hg18.2bit --2bitdir /home/aquinom
            Replacement list is longer than search list at /usr/local/share/perl/5.12.4/Bio/Range.pm line 251.
            4 125893227 + 5
            getaddrinfo() error on hostName=sjblat: No address associated with hostname

            Connection timed out
            Sorry, the BLAT/iPCR server seems to be down. Please try again later.
            10 66301865 - 9
            getaddrinfo() error on hostName=sjblat: No address associated with hostname

            Connection timed out
            Sorry, the BLAT/iPCR server seems to be down. Please try again later.
            10 66301858 + 4
            getaddrinfo() error on hostName=sjblat: No address associated with hostname

            Connection timed out


            I have gfServer running with the command: gfServer start localhost 3503 ~/hg18.2bit
            and tested that it works using gfClient. Anyone know what I need to change to get this working?

            Comment


            • #7
              Fixed issue, had to change the hardcoded values in the perl script.

              Comment


              • #8
                Hi,
                I am trying to use the new version of CREST program, CREST1-0-1. But I'm having some difficulties, I have tried to do some changes and managed to eliminate some of the errors.

                Let me know if anyone is having the same issues..

                1. on line 539 "my ($rs, $re) = ($m>0 ? $a->[$m-1][0] : $a->[0][0], $m<scalar(@{$a})-1 ? $a->[$m+1][0] : $a->[-1][0];" , the program was giving error for $m, as it was not declared
                below is the subroutine

                sub bin_search {
                my ($a, $p) = @_;
                my ($s, $e) = (0, scalar(@{$a})-1);
                # my $m = int( ($s + $e)/2);
                while(1) {
                return $s if($a->[$s][0] >= $p);
                return $e if($a->[$e][0] <= $p);
                # return $m if($a->[$m][0] == $p || ($a->[$m-1][0] < $p && $a->[$m+1][0] > $p));
                my ($rs, $re) = ($m>0 ? $a->[$m-1][0] : $a->[0][0], $m<scalar(@{$a})-1 ? $a->[$m+1][0] : $a->[-1][0];
                return $m if($a->[$m][0] == $p || ($rs < $p && $re > $p));

                if($a->[$m][0] > $p) {
                $e = $m;
                }
                else {
                $s = $m;
                }
                $m = int( ($s+$e)/2 );
                }
                }
                I had then uncommented the line 534 "my $m = int( ($s + $e)/2);"

                2. After this I got another error in the line 539 "my ($rs, $re) = ($m>0 ? $a->[$m-1][0] : $a->[0][0], $m<scalar(@{$a})-1 ? $a->[$m+1][0] : $a->[-1][0];" , syntax error at ];"
                I then added a ) at the end and changed it to my ($rs, $re) = ($m>0 ? $a->[$m-1][0] : $a->[0][0], $m<scalar(@{$a})-1 ? $a->[$m+1][0] : $a->[-1][0];" , syntax error at ]);

                3. The program now seems to work, but I'm getting another error which I'm not able to understand
                Can't open /tmp/150008.1.all.q/DmcqZ93HGj/Q4KaHbp14u.fa.cap.ace:No such file or directory at /home/ssinha/sw/CREST1-0-1/CREST.pl line 595

                Reading through the code I didn't understand where to change for the tmp folder

                Comment


                • #9
                  Hi Shruti,

                  I'm also facing the same error. Did you solve this error, please let me know.

                  can't open /tmp/aPNssgMFN0/QARvoJZlvh.fa.cap.contigs.clip.fa.psl.sorted:No such file or directory at ./CREST.pl line 606
                  Thanks,

                  Comment


                  • #10
                    Yes I managed to run it.. this was to do with some permission problem.
                    so I changed the path of the output to a location where I have full permission.

                    Comment


                    • #11
                      I'm also getting a similar error message when running CREST.pl:

                      Can't open /tmp/7iYXT0zlo_/e9TZ74LVc9.fa.cap.ace:No such file or directory at /crest/1.0/bin/CREST.pl line 593.

                      Update: It is running now that our gfServer is running and as long as I use the exact same path to start BLAT server on the .2bit file as I then use to return the prompt and start CREST.pl, similar to what other users report.
                      Last edited by dakin; 07-08-2013, 10:39 AM. Reason: working now

                      Comment

                      Latest Articles

                      Collapse

                      • seqadmin
                        The Impact of AI in Genomic Medicine
                        by seqadmin



                        Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
                        Yesterday, 02:07 PM
                      • seqadmin
                        Multiomics Techniques Advancing Disease Research
                        by seqadmin


                        New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

                        A major leap in the field has
                        ...
                        02-08-2024, 06:33 AM

                      ad_right_rmr

                      Collapse

                      News

                      Collapse

                      Topics Statistics Last Post
                      Started by seqadmin, 02-23-2024, 04:11 PM
                      0 responses
                      54 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 02-21-2024, 08:52 AM
                      0 responses
                      62 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 02-20-2024, 08:57 AM
                      0 responses
                      51 views
                      0 likes
                      Last Post seqadmin  
                      Started by seqadmin, 02-14-2024, 09:19 AM
                      0 responses
                      65 views
                      0 likes
                      Last Post seqadmin  
                      Working...
                      X