Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • tophat/cufflinks bam vs. RPKM

    Is there any legitimate reason a tophat bam file would have reads for a gene whose RPKM was not reported in the output from tophat or cufflinks?

  • #2
    Originally posted by mgogol View Post
    Is there any legitimate reason a tophat bam file would have reads for a gene whose RPKM was not reported in the output from tophat or cufflinks?
    You are finding reads in tophat bam file that are not reported in output from tophat (which is .bam format)? I don't really understand your question-can you give an example.
    I thought tophat uses FPKM units and not RKPM.

    Comment


    • #3
      If I look at the list of genes whose FPKMs are reported by cufflinks, a certain group of genes is not found in that list at all.

      If I load a bam file generated from the tophat sam file into IGV, and browse around looking at those genes, they have reads. So they have reads, but FPKMs are not reported.

      I'm just wondering if this an error or if there is some legitimate explanation?

      Comment


      • #4
        if the reads are too few, they may be filtered by the F/--min-isoform-fraction parameter where TopHat filters out junctions supported by too few alignments.
        May be you can set it to 0.

        Also, I don't intuitively understand -a/--min-anchor-length <int> parameter, but it may also be limiting your reads that are counted for the FPKM.

        Comment


        • #5
          Thanks for the ideas.

          Comment


          • #6
            sure, if you find out whats going on, let us know. If the reads you saw on IGV don't span junctions, then its probably a parameter in bowtie that is throwing it- off-low quality or reads mapping to multiple places in the genome or something else. Its a good thing you caught it though

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Strategies for Sequencing Challenging Samples
              by seqadmin


              Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
              03-22-2024, 06:39 AM
            • seqadmin
              Techniques and Challenges in Conservation Genomics
              by seqadmin



              The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

              Avian Conservation
              Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
              03-08-2024, 10:41 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, Yesterday, 06:37 PM
            0 responses
            12 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, Yesterday, 06:07 PM
            0 responses
            10 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-22-2024, 10:03 AM
            0 responses
            52 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-21-2024, 07:32 AM
            0 responses
            68 views
            0 likes
            Last Post seqadmin  
            Working...
            X