Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Help with Cufflinks in command line

    Hello,
    I am new to programming/command line and am having trouble with the simplest of tasks I am sure. I am used to using Galaxy for analysis but I am trying to go towards programming/command line usage. I am trying to run cufflinks via command line but got the following error:

    You are using Cufflinks v2.0.2, which is the most recent release.
    Error: cannot open reference GTF file /--GTF for reading



    Here is my input text into the command line(could have errors but I am not sure):

    cufflinks -G/--GTF /Users/chriscard/Desktop/Galaxy2-\[UCSC_Main_on_Cow__all_mrna_\(genome\)\].gtf /Users/chriscard/Desktop/Galaxy14-\[Tophat_high_fert__accepted_hits\].bam o- /Users/chriscard/Desktop/Command_Line_Output

    I plan on using tophat eventually as the first step to command line data analysis but for now I just tried to pull the .bam file from a tophat run off of one of my successful galaxy runs as well as the .gtf file for the reference mRNA that I want cufflinks to use(this is successful via Galaxy). Is there something wrong with my input text into command line? Do I need a different mRNA .gtf file? Sorry if this seems like a very simple task to be having problems with, trying to dive head first into this command line and coding stuff via a couple books I have started reading.

    Thank You,
    Chris

  • #2
    Good on you for leaving Galaxy and going for it on the command line.

    The reason you're seeing this error is actually explained in the error, but if you're not used to command line work it won't seem that way.

    -G/--GTF

    is being interpreted as passing the 'argument' of /--GTF to the 'flag' -G

    When you've read the command line you have assumed that you need -G/--GTF because that's what the manual says (no shame in taking it literally) however -G is the same as --GTF and the "/" just means 'use one or the other'.

    You may still run into errors by just using one flag, you seem to be also trying to escape all your bracketed characters (preceding ) and ] with \) and I'm not sure why you feel you need to do this. I'd try avoiding having filenames with these kind of characters in. I'm wondering if you're doing this because it's how the command line appears when you're running in Galaxy?

    Comment


    • #3
      Chris,

      Unix commands in general will have the following format:

      Code:
      name_of_command [-option 1 -option 2 -option n]
      Some of the commands may not need any options at all.

      Most of the times if you type the command, with -h or --help option, then you will get a list of all the possible arguments a command will accept, make a note of the default parameters shown for various cufflinks options by running the following.
      Code:
      cufflinks --help
      Check out this guide that will show you how to run the Tophat/Cufflinks suite of commands written by the program's author: http://www.nature.com/nprot/journal/....2012.016.html

      Comment


      • #4
        I can get cufflinks to run with just this in the command line: cufflinks highferttophat.bam

        However when I try and add the reference genome option in which I am typing this into the command line:

        new-hostesktop chriscard$ cufflinks [--GTF cowallmrna.gtf] highferttophat.bam


        I get this message which I am not sure how to fix

        You are using Cufflinks v2.0.2, which is the most recent release.
        open: No such file or directory
        File [--GTF doesn't appear to be a valid BAM file, trying SAM...
        Error: cannot open alignment file [--GTF for reading

        This is the same gtf file I had been using successfully as a reference gtf file to use for cufflinks through Galaxy. It seems like it has a problem with the file but I was wondering if maybe my option typing in the command line has errors I am failing to recognize? Even though this is the reference file I am using for galaxy does it need to be a different format through command line? The file type in the manual is gtf which is what this is but I keep getting this error.


        Sorry again, I am very new to this. I appreciate all of the help I have gotten and any more anyone can offer.

        Thank You,
        Chris

        Comment


        • #5
          You don't need the square brackets, so it should look something like this:
          cufflinks --GTF cowallmrna.gtf highferttophat.bam

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Recent Advances in Sequencing Analysis Tools
            by seqadmin


            The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
            05-06-2024, 07:48 AM
          • seqadmin
            Essential Discoveries and Tools in Epitranscriptomics
            by seqadmin




            The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
            04-22-2024, 07:01 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 05-14-2024, 07:03 AM
          0 responses
          19 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 05-10-2024, 06:35 AM
          0 responses
          44 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 05-09-2024, 02:46 PM
          0 responses
          54 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 05-07-2024, 06:57 AM
          0 responses
          42 views
          0 likes
          Last Post seqadmin  
          Working...
          X