Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • error message of cuffmerge

    Hi,

    When running cuffmerge, I only got very implicit error messages. It seemed that the gtf files were converted into sam files, but then stuck.

    I tried to troubleshot in several aspects: checking the name of chromosome(with chr or not); sorting the gtf files (chr name then coordinate); with larger memories. But none of them helped. I'm wondering if any of you might also meet this.

    Any suggestion will be helpful, thanks.


    =============command and message====================
    $ cuffmerge -g dmel-all-r6.09.gtf -s dmel_chr_r609.fa assemblies.txt

    [Sun Feb 21 17:29:11 2016] Beginning transcriptome assembly merge
    -------------------------------------------

    [Sun Feb 21 17:29:11 2016] Preparing output location ./merged_asm/
    [Sun Feb 21 17:29:14 2016] Converting GTF files to SAM
    [17:29:14] Loading reference annotation.
    [17:29:16] Loading reference annotation.
    [17:29:17] Loading reference annotation.
    [17:29:18] Loading reference annotation.
    [Sun Feb 21 17:29:19 2016] Quantitating transcripts
    You are using Cufflinks v2.2.1, which is the most recent release.
    Command line:
    cufflinks -o ./merged_asm/ -F 0.05 -g
    ../../../dmelGenomes/dmel-all-r6.09.gtf -q --overhang-tolerance 200
    --library-type=transfrags -A 0.0 --min-frags-per-transfr
    ag 0 --no-5-extend -p 1
    ./merged_asm/tmp/mergeSam_fileTqiEXC
    [17:29:19] Loading reference annotation.
    [17:29:23] Inspecting reads and determining fragment length distribution.
    [FAILED]
    Error: could not execute cufflinks

  • #2
    I assume you have run cufflinks already so the program must be installed and working on your machine. Simplest explanation may be that you do not have "cufflinks" program in your $PATH/current directory. You can amend your $PATH variable and then re-run.

    Comment


    • #3
      Thanks, GenoMax. But this is not the reason, I'm sure cufflinks is in the PATH. I also tried "cufflinks -o ./merged_asm/ -F 0.05 -g
      ../../../dmelGenomes/dmel-all-r6.09.gtf -q --overhang-tolerance 200
      --library-type=transfrags -A 0.0 --min-frags-per-transfr
      ag 0 --no-5-extend -p 1" separately, the message is more concise:

      =========
      [11:45:13] Loading reference annotation.
      [11:45:17] Inspecting reads and determining fragment length distribution.
      Segmentation fault

      Comment


      • #4
        It looks like cufflinks is seg faulting/crashing. Which could be because of many reasons. Have you used this copy of cufflinks before?

        Comment


        • #5
          The version of cufflinks is v2.2.1. I loaded it as a module in a computing cluster (and all commands of the cufflinks package were loaded), so I think the software is correctly located and functional.

          I mainly used cuffdiff to estimate differently expressed annotated genes before, and didn't use cuffmerge often.

          Comment


          • #6
            Not sure what to tell you in this case. Since someone else installed this program we should assume that they built the program from source or got the right binary for this OS.

            You are sure there is no other cufflinks in $PATH from an older version that may be getting preferentially executed.
            Last edited by GenoMax; 02-29-2016, 09:31 AM.

            Comment


            • #7
              I think this is the only version available. Thank you, GenoMax, all the same.

              Comment


              • #8
                Hi GenoMax, I finally found the problem is caused by the cufflinks (v2.2.1) in the cluster. And when using an earlier version (v2.2.0) by myself, it all went well.

                Thanks, and sorry for bothering.

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Best Practices for Single-Cell Sequencing Analysis
                  by seqadmin



                  While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
                  06-06-2024, 07:15 AM
                • seqadmin
                  Latest Developments in Precision Medicine
                  by seqadmin



                  Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

                  Somatic Genomics
                  “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
                  05-24-2024, 01:16 PM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 06-07-2024, 06:58 AM
                0 responses
                179 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-06-2024, 08:18 AM
                0 responses
                227 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-06-2024, 08:04 AM
                0 responses
                184 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-03-2024, 06:55 AM
                0 responses
                18 views
                0 likes
                Last Post seqadmin  
                Working...
                X