Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • prbndr
    Junior Member
    • Sep 2011
    • 5

    #1

    cufflinks rna-seq coverage

    i'm a little confused as to why cufflinks didn't output a coverage estimate in the gene level tracking file. does anyone know what the reason is for this, and how i can derive an estimate of gene level coverage from the cufflinks output files?
  • zhangju
    Member
    • May 2011
    • 18

    #2
    Cufflinks coverage no output

    I run Cufflinks with no coverage output. Have you got the answer yet? Let me know. Thanks, Justin

    Comment

    • zhangju
      Member
      • May 2011
      • 18

      #3
      Cufflinks coverage from Tophat

      Originally posted by zhangju View Post
      I run Cufflinks with no coverage output. Have you got the answer yet? Let me know. Thanks, Justin
      I think I figure it out. Use "--coverage-search" option at Tophat run. I am rerunning it now.

      Comment

      • rflrob
        Member
        • May 2010
        • 50

        #4
        Originally posted by zhangju View Post
        I think I figure it out. Use "--coverage-search" option at Tophat run. I am rerunning it now.
        Did this work for you? I'm also failing to get any gene-level coverage information (though it does work at the isoform level), and I usually run my tophat with --no-novel-juncs, so I'm a little hesitant to turn on --coverage-search, but if that's what I have to do...

        Comment

        • zhangju
          Member
          • May 2011
          • 18

          #5
          This is my command line

          "tophat -p ${NUM_PROCS} -o $OUTPUT --coverage-search -G $GTF --transcriptome-index=$TRANSINDEX $BINDEX $READ1 $READ2 $Orphan"

          It works for me.

          Originally posted by rflrob View Post
          Did this work for you? I'm also failing to get any gene-level coverage information (though it does work at the isoform level), and I usually run my tophat with --no-novel-juncs, so I'm a little hesitant to turn on --coverage-search, but if that's what I have to do...

          Comment

          • zhangju
            Member
            • May 2011
            • 18

            #6
            You would try a small subset data to test this command. I could not remember well. I check my most current command without "coverage-search", so "coverage-search" may not help. I am running Bacterial genome, so I am using isoform data as gene data. I did not set "no-novel-junc" on purpose, and it did not find new juncs by program.

            Originally posted by zhangju View Post
            This is my command line

            "tophat -p ${NUM_PROCS} -o $OUTPUT --coverage-search -G $GTF --transcriptome-index=$TRANSINDEX $BINDEX $READ1 $READ2 $Orphan"

            It works for me.

            Comment

            Latest Articles

            Collapse

            • SEQadmin2
              Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
              by SEQadmin2



              Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

              Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
              ...
              Today, 09:55 AM
            • SEQadmin2
              New Genomics Technologies Take Aim at Long-Standing Limits
              by SEQadmin2


              Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

              We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
              ...
              09-28-2026, 10:25 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, Today, 05:39 AM
            0 responses
            9 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, Yesterday, 09:56 AM
            0 responses
            11 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 10-09-2026, 11:14 AM
            0 responses
            25 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 09-29-2026, 09:51 AM
            0 responses
            77 views
            0 reactions
            Last Post SEQadmin2  
            Working...