Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • crh
    Member
    • Dec 2009
    • 46

    cuffdiff Error: cannot open reference GTF file –L for reading

    Hi

    I ran tophat with replicate samples for a control & treated sample, then ran cufflinks, followed by cuffmerge. I provided a GTF file to cuffmerge.

    When I attempt to run cuffdiff I get the following error:
    You are using Cufflinks v1.3.0, which is the most recent release.
    Error: cannot open reference GTF file –L for reading


    my command line is:
    cuffdiff -o diff_out -b genome.fa –L C,T -u merged_asm/merged.gtf ./C1/accepted_hits.bam,./C2/accepted_hits.bam ./T1/accepted_hits.bam,./T2/accepted_hits.bam

    Any suggestions?

    thanks

    Charles
  • mpiro
    Member
    • Nov 2009
    • 10

    #2
    the dash sign behind the L seems distorted .. I assume you copied the command from a web page .. please delete the -L and retype it.

    Comment

    • crh
      Member
      • Dec 2009
      • 46

      #3
      Correct you are - retyping seems to have fixed the error - thanks!

      Charles

      Comment

      • Charitra
        Member
        • Feb 2013
        • 57

        #4
        I have the same problem but after removing -L doesnt work.
        my commands are:
        cuffdiff -o diff_out -b genome.fa -p 8 calu,h358 -u merged.gtf, \ accepted_hits.bam accepted_hits1.bam
        or
        cuffdiff -o diff_out -b genome.fa -p 8 calu,h358 -u merged.gtf, \ ./accepted_hits.bam, accepted_hits1.bam

        Error: cannot open reference GTF file merged.gtf, for reading

        My both .bam files are in the same working directory but still the error ?
        ?????

        Comment

        • jeales
          Member
          • Oct 2012
          • 13

          #5
          You do need the "-L" before your labels
          But sometimes webpages, text editors and Word can insert non-standard hyphen-like characters like this "–", these then get copied and pasted in command lines and scripts and cufdiff etc have no idea what they are
          Just type out the command fully yourself and use the proper hyphen button on your keyboard

          Comment

          Latest Articles

          Collapse

          • GATTACAT
            Reply to Nine Things a Sample Prep Scientist Thinks About Before Sequencing
            by GATTACAT
            Love this - good data definitely starts from good input, and poor input can only give relatively poor data. I particularly like the mention of Nanodrop/absorbance based methods for quantification. It's such a toss up if you'll get an accurate reading or what amounts to a randomly generated number, and a lot of library/sequencing related issues can be traced back to poor quant.
            07-01-2026, 11:43 AM
          • SEQadmin2
            Nine Things a Sample Prep Scientist Thinks About Before Sequencing
            by SEQadmin2


            I’m not a sequencing expert. I’m a purification scientist who uses NGS to evaluate workflows my group develops. With this perspective, we think about the sample first and the NGS workflow second. The sequencer is an exceptionally honest reporter, but it can only report on what you give it, so whether you get clean, interpretable data from an NGS workflow is largely determined before you begin.

            Here are nine questions we think about, in roughly the order they matter, before...
            06-18-2026, 07:11 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, 07-02-2026, 11:08 AM
          0 responses
          14 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 06-30-2026, 05:37 AM
          0 responses
          15 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 06-26-2026, 11:10 AM
          0 responses
          20 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 06-17-2026, 06:09 AM
          0 responses
          54 views
          0 reactions
          Last Post SEQadmin2  
          Working...